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Chapter 1 

Outline 



Quantum chromodynamics (QCD) is the quantum field theory of the strong interac- 
tion, with quarks and gluons being its elementary degrees of freedom. Nevertheless 
only color-neutral hadrons as bound states of quarks and gluons appear as observable 
particles in detector facilities, and a detailed study of hadron properties is of fundamen- 
tal importance for an understanding of the quark-gluon dynamics of QCD. A wealth of 
information has been collected in the past decades' high-energy scattering experiments 
but many issues are still unresolved. Present and future experimental programs at 
JLAB, SLAC, COMPASS/CERN and FAIR/GSI are devoted to charmed meson and 
baryon spectroscopy, the search for exotic mesons, or the investigation of the nucleon's 
electromagnetic and spin structure in terms of form factors and generalized parton 
distributions. 

The interplay between QCD's elementary and observable degrees of freedom ad- 
dresses two phenomena whose origin is not yet fully understood: color confinement 
and dynamical chiral symmetry breaking. The latter is the mass generation mecha- 
nism that equips light quarks with large dynamical constituent-quark masses whereas 
it retains the light pseudoscalar mesons as would-be Goldstone bosons of spontaneously 
broken chiral symmetry. While QCD is well under control in the high-energy region 
whose weak coupling strength enables the application of perturbative methods, these 
phenomena characterize the large-distance or low-energy structure of QCD where the 
coupling is strong and demands a non-perturbative treatment. 

Many approaches which describe particular aspects of QCD have emerged in the 
past decades. For instance, quark models have pioneered our understanding of hadron 
structure and successfully described many of their properties, see e.g. [1-5]; effective 
field theories provide rigorous results in certain limits, e.g. [6-11]; experimentally ex- 
tracted generalized parton distributions (GPDs) combine electromagnetic form factors 
and parton distribution functions and establish a three-dimensional tomography of 
hadrons, see [12-14] for recent reviews. Each has its own strengths but also weak- 
nesses, e.g., parameter dependence, limited range of applicability, or the reliance upon 
a factorization scale. 
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Outline 



The long-term goal of rigorously solving QCD necessitates a quantum-field theo- 
retical, non-perturbative description starting directly from the QCD Lagrangian as a 
prerequisite. Lattice QCD is constantly pushing forward towards the physical quark 
mass in its investigation of hadron structure [15, 16]. On the other hand, the increase 
of computational power in the past decade has also played a vital role in the devel- 
opment of functional continuum methods: the Dyson-Schwinger equations (DSEs) of 
QCD provide a tool to map out the infrared structure of QCD where confinement 
and dynamical chiral symmetry breaking occur; non-perturbative bound-state equa- 
tions have been utilized to perform hadron spectroscopy and thereby determine the 
meson and baryon amplitudes that are needed to calculate hadronic observables such 
as electromagnetic form factors and decay properties. For reviews on Dyson-Schwinger 
equations and their application to hadron physics, see e.g. Refs. [17-21]. 

Functional methods are in some sense complementary to lattice QCD. While access 
to the full momentum region and quark-mass range is readily available without any 
need for extrapolations, one faces an infinite system of integral equations which needs 
to be truncated with regard to a practical numerical computation. Such truncations, 
while subject to specific constraints, induce a model dependence which is not easily 
quantified. Lattice calculations and their chiral extrapolations are not only useful 
for comparing results but can also serve as a guideline to construct truncations which 
capture the important physics. Much of the current effort is focused on further resolving 
the dynamics of quarks and gluons from their DSEs and elevate simplistic truncations 
towards a more complete picture of the involved physical content. An exploration of 
the full potential contained in such functional approaches has only begun, and the 
future holds many exciting perspectives: from a phenomenologist's point of view, some 
of the intermediate goals are access to charm and bottom physics, the investigation of 
excited hadrons, electromagnetic form factors in the timelike region, or the medium- 
and large-Q^ domain of spacelike form factors. 

The aim of this thesis is to report progress that has been achieved in the framework of 
meson and baryon bound-state equations and thereby contributes to an understanding 
of the dynamics of hadrons from their underlying constituents. After recapitulating 
basic concepts in QCD and the bound-state formalism in Chapter 2, we provide a 
systematic study of quark-antiquark (Chapter 3), three-quark (Chapter 4), diquark 
and quark-diquark bound states (Chapter 5). The common parenthesis is provided 
by a rainbow-ladder (RL) truncation, i.e. an iterated vector-vector gluon exchange in 
the qq kernel, which guarantees the correct implementation of chiral symmetry and its 
spontaneous breaking and is elaborated in Chapter 3. In the qq sector we report results 
for vr and p-meson observables; the three-body equation is solved for the nucleon, and 
the quark-diquark approach is applied to A'^ and A baryons. In Chapter 6 we present 
results for the nucleon's electromagnetic form factors in the quark-diquark framework. 
In Chapter 7 we summarize and conclude. 

Calculational details, as well as details on the structure of Green functions and 
bound-state amplitudes that appear throughout the text, are collected in Appendices A 
and B. We work in the Euclidean formulation (corresponding relations are given in 
App. B.l) and restrict ourselves to two-flavor QCD with isospin symmetry. 



Chapter 2 

Bound states in QCD 



The dynamical content of QCD as a local quantum field theory of quarks and gluons 
is described by its Lagrangian or, equivalently, its action. It is a difficult task to 
unfold the physical content of the QCD Lagrangian toward a quantum-field theoretical 
description of hadrons. At large energies and small distances (the 'ultraviolet' region), 
the interaction between quarks and gluons is weak and can be described by perturbative 
methods. At low energies and large distances (the 'infrared'), the coupling becomes 
strong and perturbation theory can no longer explain the most striking features of QCD 
in this regime: the dynamical breaking of chiral symmetry (D^SB), the generation of 
large constituent-quark masses from almost massless quarks, the formation of hadrons 
— mesons as quark-antiquark and baryons as three-quark bound states — , and the 
confinement of quarks and gluons inside hadrons. 

A rigorous starting point to describe the non-perturbative dynamics of quarks and 
gluons is the path-integral approach. Dyson-Schwinger equations (DSEs) are the quan- 
tum equations of motion. They relate QCD's Green functions — the basic propagators 
and vertices of the theory, e.g. the dressed quark and gluon propagators, and the 
quark-gluon vertex — to each other and lead to a self-consistent and infinitely coupled 
system of integral equations. 

Hadrons appear as free-particle poles in the respective n-particle Green functions. 
For instance, the quark-antiquark and three-quark Green functions exhibit meson and 
baryon poles, respectively. Hadron properties can be extracted upon solving bound- 
state equations which are valid at these poles and need the Green-function content of 
QCD as an input. In combination they provide a powerful tool to calculate experi- 
mentally accessible hadron observables, e.g. meson and baryon mass spectra, decay 
constants, scattering processes, and electromagnetic properties such as form factors, 
magnetic moments and charge radii. 

In the present chapter we briefly introduce the theoretical foundation and phe- 
nomenological aspects of QCD. Starting from its generating functional, we sketch the 
derivation of Dyson-Schwinger equations, discuss confinement and D^SB and their 
manifestation in the bound-state approach, and derive the general form of a bound- 
state equation which will be referred to in the following chapters. 
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Bound states in QCD 



2.1 Basic concepts in QCD 



QCD action. Based upon the principle of local gauge invariance, the QCD action 
which describes the interaction of quark and gluon fields "0) is written as 

Sqcb[A, V, ^ ] = jd^'x [i, {-p + m)i; + l F^, F^,] , (2.1) 

where = + igA^ is the covariant derivative and the gluon field-strength tensor 
reads -F^i, = df^A^, — d^A^ + ig [A^, A^]. The action is per construction invariant along 
the gauge orbit 

V^^ = C/V, A^ = UA^U^ -'-U{d^U^), Fl^, = UF^,U\ (2.2) 

where U{x) £ SU{3)c is a local gauge transformation and A^ = A'^t"', F^y = F^^jt"' 
are elements of the corresponding color algebra whose basis elements i", a = 1 ... 8 
satisfy the commutator relation = if"'^'^t^. In the fundamental representation 

they are given by the Gell-Mann matrices = X"'/2. 

The non-Abelian nature of the color group SU(3)c induces gluonic self-interaction 
terms ~ A^ and ~ A^ encoded in Ffj^jy which lead to significant complications compared 
to the Abelian gauge theory QED. They are commonly believed to be the origin of 
phenomena such as confinement and dynamical chiral symmetry breaking in QCD. 

Generating functional. In an Euclidean path-integral approach, the quantum field 
theory which corresponds to the action (2.1) is defined by the generating functional 

Z[J,rj,fj,a,a]= V, V^, c, c] e" ^QcdIA.V.V^] ~ 5gp[A,c,c] + 5c (2.3) 



from which all physical quantities can be derived. An integration over infinitely many 
physically equivalent gauge field configurations A'-' , and the emergence of zero eigen- 
values of the perturbative inverse gluon propagator as obtained from (2.1), are avoided 
by adding the gauge-fixing term 



Sgf[A,c,c] = j d!^x 



^%|^ + (a,c-)«c^) 



(2.4) 



to the classical action^. It introduces unphysical auxiliary fields, the scalar anticom- 
muting ghost fields c and c, where D'^^ = 5°'^d^ + gf°'^'^A'^ is the covariant derivative in 
the adjoint representation. The expression (2.4) arises from evaluating the gauge-fixing 
condition 5(9^^^ — S^B) det M[74] in the path integral which involves the determinant 
of the Faddeev-Popov operator 

M[Ar\x,y) = -d,D';,''5\x-y). (2.5) 



Gauge fixing is not necessary for a direct calculation of gauge-independent quantities via Eq. (2.7) 
as is for instance routinely done in lattice QCD. 
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The source term 



Sc 



(fx 



A^JI, + fjij) + r] + a c + ca 



(2.6) 



contains the external currents which are auxiliary quantities to enable the derivation 
of Green functions in terms of functional derivatives of Z. 

Green functions. All physical properties can be extracted from the Green functions 
of a quantum field theory, and the theory is solved once all of them are determined. 
Abbreviating the fields generically by (p{x) and the associated sources by J(x), the 
Green functions are defined as the time-ordered vacuum expectation values of products 
of fields, denoted by G[ip\: 



_5_ 
Jj 



=0 m 



(2.7) 



For instance, a two-point function is associated with G[ip\ = ip{x)ip{y). In the functional 
formulation the Green functions correspond to averages over all field configurations with 
a probability distribution e"'^!'^^. As indicated in (2.7), they are most conveniently 
obtained as moments of the generating functional Z[J] by taking functional derivatives 
with respect to the sources J. The same procedure yields connected Green functions 
as derivatives of the functional W[J] := — lnZ[J] and one-particle irreducible (IPI) 
Green functions from the effective action which is related to W[J] via a Legendre 
transformation: 



Z[J] 



-W[J] 



Dip e 



-S[<f]+Jf{x) J(x) 



= : e 



(2i 



The averaged field if is the expectation value of if in the presence of the source J: 
6W[J] _ /P(/7e-^M+/¥'WJW(^(a;) 



{(fix)) J = ip{x) 



SJ{x) j-p^e-5M+/^(x)j(x) 
With the following shorthand notation for functional derivatives: 



W'Mxy 



5W[J] 



6J{x) 5J{y) 



r"[o] 



6m 



J=0 



xy 



6ip{x)6ip{y) 



(2.9) 



(2.10) 



the propagator associated with the field ip is given by r"[0]~j^^, the respective three-point 
vertex is r"'[0]2.y^, etc. IPX vertices involving different types of fields are obtained upon 
differentiating with respect to each corresponding variable ifi{x). 

Dyson-Schwinger equations. Dyson-Schwinger equations [22, 23] are the quantum 
equations of motion. They follow from an invariance of the generating functional under 
a variation ip{x) — > ip{x) + e{x) of the fields. Assuming that the integral measure is 
invariant under such a transformation, the condition Z'[J] = Z[J] leads to the relation 



\Sip{x)/ J 



5ip 



6J{x) 



Z[J] = J{x) 



(2.11) 
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Bound states in QCD 



for each field 92 that appears in the action. This is the defining relation for the infinite 
tower of Dyson-Schwinger equations which relate the Green functions of the field theory 
to each other. Inserting Z[J] = e~^^'^ yields the corresponding equation for connected 
Green functions: 



6J{x) 



Sip 

With W'[J]^ = -^{x), r'[^^ = J(x), and 

6 f 5^{y) 5 f „ 6 



1 = J{x) . (2.12) 



J 6J{x) difiy) J ^ difiy) J ^^^^y 6ip{y) ^ ' 



5J{x) 

one arrives at the generating DSE for IPI Green functions: 



y 




Further (n — l)-fold differentiation and finally setting if = yields the system of DSEs^ 
for the IPI n-point functions r(")[0]a;i,...x„- 

Dyson-Schwinger equations in QCD. In the context of QCD, the action whose 
functional derivative appears in Eqs. (2.11-2.14) is the sum of the classical action and 
the gauge-fixing contribution (2.4). Dyson-Schwinger equations for IPI Green functions 
are obtained via functional derivatives of Eq. (2.14) with respect to the gluon, quark 
or ghost fields A, -0, tjj, c and c, where the latter four are anticommuting Grassmann 
variables. A detailed derivation of these relations can be found in Refs. [17, 18,25]. 
The DSEs for quark, gluon and ghost propagators and quark-gluon and ghost-gluon 
vertices are illustrated in Fig. 2.1. 

The Dyson-Schwinger approach provides an appealing tool for several reasons. DSEs 
operate in fully relativistic quantum field theory; they provide access to both perturba- 
tive and non-perturbative regimes of QCD; and they represent a continuum approach 
which is able to cover the full quark-mass range between chiral limit and the heavy- 
quark domain. Of course the main caveat concerns the complexity of this framework: 
for numerical studies one relies upon a truncation of the infinite systems of equations to 
a subset that captures the physical content and is solved explicitly, combined with the 
use of ansatze for those Green functions that enter the equations but are not solved for. 
These ansatze are constrained by symmetry properties, multiplicative renormalizability, 
perturbative limits, etc. During the past years a cross-fertilization between functional 
methods and lattice QCD has provided further insight into the non-perturbative struc- 
ture of Green functions. 

Gauge invariance. While Green functions depend on the gauge, physical observ- 
ables must be gauge independent. Similarly as the classical equations of motion are 



Mathematica package which enables an automated derivation of DSEs is described in Ref. [24]. 
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Figure 2.1: DSEs for quark, gluon and ghost propagators and the ghost-gluon and 
quark- gluon vertices. White blobs denote dressed propagators, filled blobs represent 
IPI Green functions. 



implemented in the functional description through equations for the IPI vertices, the 
requirement of gauge invariance of the generating functional leads to Ward-Takahashi 
(WTIs) and Slavnov- Taylor identities (STIs) which interrelate the Green functions of 
a gauge theory. Although the gauge-fixed action of QCD is no longer invariant with 
respect to local gauge transformations, it still satisfies global gauge symmetry and 
BRST symmetry [26, 27] . The latter is formally equivalent to a gauge transformation 
by a ghost field and can be used to derive the STIs. Moreover it can be shown that 
requiring BRST invariance of a gauge theory generates both the ghosts and the gauge 
fixing while ensuring gauge independence of physical observables [28]. 
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Bound states in QCD 



2.2 Phenomenological aspects of QCD 

Confinement. The absence of colored states in the physical spectrum is one of the 
most characteristic non-perturbative phenomena in QCD, and its underlying mecha- 
nism is still not fully understood. Several different ideas have emerged which are likely 
to explain different aspects of confinement; for reviews, see Refs. [18,21,29-32]. A 
common line of thought builds upon the idea of certain topological gauge field con- 
figurations dominating the path integral near the Gribov horizon which corresponds 
to a vanishing eigenvalue of the Faddeev-Popov operator [33] . This entails a linearly 
rising quark-antiquark potential in Coulomb gauge [34-39] and a strongly infrared- 
divergent ghost propagator and infrared-vanishing gluon propagator in Landau gauge 
which trigger the infrared behavior of other Landau-gauge Green functions [40-42] . 

A somewhat different point of view concerns the identification of physical subspaces 
from the asymptotic state space of QCD. The Kugo-Ojima scenario [43] which relies 
on BRST symmetry describes the cancellation of longitudinal and timelike gluons with 
ghost and antighost fields similar to the Gupta-Bleuler mechanism in QED. On the 
other hand, the Osterwalder-Schrader axiom of reflection positivity [44] implies that a 
certain degree of freedom whose propagator violates positivity and thereby prevents a 
Lehmann representation cannot describe an asymptotic physical state. 

The manifestation of confinement in the bound-state framework is elusive. In the 
particular approach employed in this work, ansatze are employed for the gluon prop- 
agator and quark-gluon vertex such that, out of the coupled system of DSEs, the 
quark propagator is the only Green function which is explicitly solved for. While a 
rainbow-ladder truncation indeed induces complex conjugate poles which ensure posi- 
tivity violation and describe a confined quark, this result is truncation dependent and 
sensitive to the structure of the quark-gluon vertex [45,46]. 

On the other hand, the solution of the quark DSE is insensitive to the precise shape 
of the interaction in the deep infrared region since the self-energy integral is dominated 
by momenta larger than the quark mass. As a consequence, the impact of infrared 
physics upon hadronic ground states is at best modest. It should certainly become 
important in the context of highly excited states with a large spatial extension or 
small-x physics. 

Dynamical chiral symmetry breaking and bound states. A phenomenon which 
has a direct impact upon the spectra of light hadrons is the spontaneous breaking of 
chiral symmetry. An unbroken chiral SUiiNf) x SUpi^Nf) symmetry, realized in the 
massless QCD Lagrangian, would imply mass-degenerate meson parity doublets in the 
chiral limit whose remnants should be visible in nature. The surprisingly small mass 
of the pion compared to its parity partners indicates that chiral symmetry is broken 
spontaneously, and the pion is identified as the massless Goldstone boson of the two- 
flavor case in the chiral limit of massless quarks. It acquires a small mass due to the 
explicit breaking of chiral symmetry at small non-zero current-quark masses, a behavior 
which is described by the Gell-Mann-Oakes-Renner (GMOR) relation [47]. 
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An order parameter of chiral symmetry breaking is the scalar quark condensate. It 
modifies the structure of the QCD vacuum and, through its interaction with the quarks, 
equips them with a large dynamical constituent mass. This behavior is reflected in the 
quark mass function M(p^) that appears in the dressed quark propagator: it connects 
the perturbative ultraviolet momentum region with the nonperturbative infrared and 
thereby communicates the transition from a current quark to a dynamically generated 
'constituent quark'. Since quarks are the building blocks of hadrons, the masses of 
light mesons and baryons are generated dynamically as well. DxSB therefore offers an 
explanation for the quark mass-generation mechanism, and it illustrates why the simple 
constituent-quark reckoning mMeson ~ 2 M(0), rriBaryon ~ 3 M(0) works reasonably well 
except for the light pseudoscalar mesons. 

D^SB becomes manifest not only in the quark propagator but also in the quark- 
gluon vertex whose scalar Dirac structures are dynamically generated together with 
the quark mass function [42]. Such an effect is naturally missing upon employing the 
rainbow-ladder truncation. Here the quark-gluon coupling a{k'^) associated with a bare 
vertex enters the quark DSE as a parametrization including a scale Air as its input 
(Sections 3.2 and 3.3). A non-zero quark mass function in the chiral limit only occurs 
if the coupling exhibits enough strength in the infrared. Above a critical threshold. Air 
is directly proportional to all mass-dimensionful quantities in the chiral limit, hence all 
of them represent a 'scale of D^SB'. 

Pion cloud. It is a longstanding prediction that, induced by D^SB, the low-energy 
and low quark-mass behavior of hadrons is modified by their interaction with pseu- 
doscalar mesons, i.e. the long-range part of qq correlations. Established in the cloudy 
bag model [48-50], where the pion field is coupled to a constituent-quark bag [51], 
meson-cloud effects have been studied in a number of quark models, e.g. [3,52-56]. 
They are systematically implemented in chiral effective field theories [6, 57] which, in 
combination with lattice simulations, provide an efficient tool for describing masses and 
electromagnetic properties of hadrons [8,58,59]. 

In these frameworks hadrons consist of a 'quark core' that is augmented by a pseu- 
doscalar meson cloud which mediates a stronger binding, decreases the hadron's mass 
and increases its size. Chiral effective field theories combined with lattice techniques 
typically predict a reduction of 20 — 30% for dynamically generated hadron masses by 
chiral corrections, an effect which is suppressed with increasing distance from the chiral 
limit. 

In the covariant bound-state approach used herein, the freedom of choosing a current- 
mass dependent rainbow-ladder coupling strength can be exploited to construct a 
hadronic quark core which subsequently needs to be dressed by meson-cloud effects [60] . 
Following up on previous quark-diquark model investigations [61-63], we will frequently 
present results in such a 'core model' and compare to those of chirally extrapolated 
lattice calculations. A natural extension is the explicit implementation of pionic ef- 
fects in the covariant bound-state equations; corresponding results have been recently 
reported [64]. 
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2.3 Bound- state equations 

T-matrix and scattering kernel. As bound states of valence quarks, hadrons cor- 
respond to poles in the quark 4- or 6-point functions G*^^^ and G^^^ or, equivalently, in 
their amputated connected parts, the scattering matrices T^^^ and T^^-* defined via 

^(n) ^ ^ c^^^T^''^G^^\ (2.15) 

(n) 

where Gq is the product of n dressed quark propagators Si. We dropped all Dirac, 
color and flavor indices of each quark leg in the above quantities for brevity. The 
product in Eq. (2.15) is understood as a summation over all occurring indices as well 
as integration over all internal momenta. The full notation is given at the end of this 
section. 

The n-quark T-matrix is related to the n-quark scattering kernel K^"'^ {n > 2) by a 
non-perturbative Dyson sum. In the context of bound-state equations one frequently 
encounters renormalization-group invariant combinations of T^"-* or K^'^'^ with n quark 
propagator legs attached on their right-hand sides. Denoting them by 

f (") := t(")gJ"\ ^(n) — ^HG-W^ (2.16) 
the defining relations for the scattering kernels ET^") read 

J'{n) _ _|_ j^(n) _|_ j^{n) j^{n) j^{n) _|_ _ _ (2.17) 

which may be rephrased as an integral equation, namely Dyson's equation (also referred 
to as inhomogeneous BSE) for the T-matrix: 





1 = 1 


(i + r(")^ 


1 . (2.18) 



This equation provides the central foundation of the approach and is depicted in the 
upper part of Fig. 2.2. It allows for a derivation of bound-state equations for qq, qq 
and qqq systems together with their canonical normalization conditions, cf. Eq. (2.23). 
Moreover, it is of virtue when constructing an off-shell ansatz for the 2-quark T-matrix 
(App. A. 4), and we will resort to Eq. (2.18) for deriving an electromagnetic current 
operator (Section 6.1). Schematically, its inverse form reads: 

(^r(")) = (/^(")) - gJ") ^ (f = - 1 . (2.19) 

The scattering kernels X^") consist of /-quark irreducible components, with / = 
2 . . . n. For instance, the three-body kernel K^^'f which appears in the bound-state 

(3) 

equation of a baryon is the sum of a 3-quark irreducible contribution K>J and three 
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Figure 2.2: Schematic derivation of a two-body bound-state equation. The first row 
illustrates Dyson's equation (2.18). The behavior at the mass pole defines the bound- 
state amplitude and leads to the corresponding bound-state equation (second row). 



permuted 2-body kernels ® [65-67]. With the notation of (2.16), the kernel 
K^'^^ reads 

^(3)=i^(3) + ^Kf), (2.20) 

i=l 

where the subscript i identifies the spectator quark. K^^"^ is illustrated in Fig. 2.3. 



Bound-state equations. At the pole corresponding to the bound-state mass M, 
bound-state amplitudes ^ are introduced as the residues of the scattering matrix via 



rp(n) 



(2.21) 



where P is the total momentum of the n quarks. The possibly dimensionful constant 
Af accounts for the dimensionality of T*^") and depends on the spin of the resulting 
particle. For instance, the propagators of free spin-0 and spin-i/2 particles are given 
by: 



J = 0: 



1 



P2 + M2 ' 



J = 1/2 : 



-if + M 



2M 



A+(P) 

P2 + M2 



(2.22) 



For a scalar or pseudoscalar particle: M = 1. In the spin-1/2 case, the matrix- valued 
amplitude ^' includes the positive-energy projector A+(P) = (1 + ^)/2 (cf. Section 4), 
where P denotes the normalized total momentum; this yields M = 2M. 

Inserting the pole condition (2.21) into Dyson's equation and comparing the residues 
of the most singular terms leads to a bound-state equation at the pole P^ = —M"^, cf. 
Fig. 2.2. An examination of the relation T' = —T{T~^)'T at the bound-state pole, 
where ' denotes the derivative d/dP'^, yields the associated canonical normalization 
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Figure 2.3: Three-quark kernel K^^^ of Eq. (2.20). 



condition. Their combination completely determines the amplitudes ^ on the mass 
shell: 



^ = i^W^, (2.23) 



^' = 1. (2.24) 



Eq. (2.23) is a fully relativistic linear homogeneous integral equation. It is the Bethe- 
Salpeter equation in the two-body case (n = 2) and its quantum-field theoretical ana- 
logue for a three-body system (n = 3). Its solution necessitates knowledge of the quark 
propagator Si and the kernel K^^\ Both ingredients can in principle be obtained 
from the infinite coupled set of Dyson-Schwinger equations, cf. Section 2.2. Feasible 
present-day numerical DSE solutions include 2- and 3-point functions within certain 
truncations, but the complexity of DSEs for 4-point functions has so far prevented a 
sophisticated numerical treatment. 

The construction of appropriate kernels is restricted by the underlying symme- 
tries of the theory. Symmetries in quantum-field theory are implemented by Ward- 
Takahashi (WTIs) and Slavnov- Taylor identities (STIs) which relate different n-point 
functions to each other. A prominent example is the axial-vector Ward-Takahashi iden- 
tity (AVWTI) [68] which relates the two-quark kernel K^"^^ to the kernel of the quark 
DSE. It is imperative to satisfy these identities in the truncation used in a numeri- 
cal study, and it will be the guiding principle which motivates an application of the 
rainbow-ladder truncation in Section 3.1. 

The bound-state approach as an eigenvalue problem. The linear homogeneous 
integral equation (2.23) can be viewed as an eigenvalue problem for the kernel i^^"^: 

^(")(P2)M/. = A,(P2)M/., (2.25) 

where P is the total momentum of the n-quark bound state and enters the equation as 
an external parameter. Upon projection onto given quantum numbers, the eigenvalues 
of K^^^ constitute the trajectories \i{P'^). An intersection Ai(P^) = 1 at some value 
= —Mf reproduces Eq. (2.23) and corresponds to a potential physical state^ with 
mass Mi. 

^ In this context one has to keep in mind the possibihty of anomalous states in the excitation spectra 
of BSE solutions, see e.g. [69]. 
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To enable an iterative solution, the spectrum of i^^") must be bounded. The largest 
eigenvalue Aq of i^^") represents the ground state of the quantum numbers under con- 
sideration and the remaining ones Ai>i its excitations; the associated eigenvectors 
are the bound-state amplitudes. A repeated multiplication by itself is equivalent to a 
projection onto the largest eigenvalue Aq, i.e., the ground state. Solving for the first 
excited state involves a subtraction of the ground-state contribution from the kernel. 
Hence the are obtained upon applying the projectors 



lim 

m— >oo 



Vi := lim 



- Ao|^o)(^c 



(2.26) 



onto a general amplitude ^. As an alternative, one may identify the ground state and 
excited states from their poles in the off-shell vertex whose quantum numbers coincide 
with those of the bound-state amplitude. This is realized by solving an inhomogeneous 
Bethe-Salpeter equation for the respective vertex [68, 70, 71]. 

Solution strategy. In this thesis we will merely be concerned with ground states. 
To solve Eq. (2.23), one must specify the color, flavor and spin structure of the bound- 
state amplitude under investigation. The latter is constructed to be Poincare covariant 
and expressed through a certain number of matrix-valued Dirac amplitudes which 
incorporate all involved momenta. The respective dressing functions depend on the 
Lorentz-invariant scalar products of these momenta. For instance, a two-body ampli- 
tude can be characterized by a relative momentum p and a total on-shell momentum P: 



*(p,p) 



fi{p^,p- P,P^)Ti{p,P) (g) Color (g) Flavor, 



(2.27) 



Projecting the bound-state equation onto its color and flavor quantum numbers and 
the orthogonal Dirac basis results in coupled integral equations for the components /j. 

According to Eq. (2.25), the equation can be solved via iteration within a 'guess 
range' P^ £ {~^mm^ ~^ma.x}i where Mmax is determined from the singularity struc- 
ture of the equation's ingredients, for instance the quark propagator (see App. B.3). 
The eigenvalue A(— M^) = 1 determines the bound-state mass M. These procedures 
are explained in detail, e.g., in [72] in the context of a quark-diquark Bethe-Salpeter 
equation. 

The current-quark mass is an input to the quark DSE, cf. Section 3.1, and can be 
mapped onto the pion mass upon solving the pseudoscalar meson BSE. This allows for 
a determination of all subsequent results as a function of m^, where the physical point 
is characterized by = 138 MeV. Varying the current mass, and thus the pion mass, 
enables a direct comparison to lattice data and their chiral extrapolations. 



Remarks on the notation. In the current section we have dropped almost all po- 
tential occurrences of notational inconvenience, a strategy which we will continue to 
pursue when deriving formal relations. In this formal notation, internal loop integrals. 
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summations or four-dimensional delta functions are not retained; possible color-flavor 
prefactors or symmetrization factors are not considered either. They will be stated in 
subsequent chapters when the expressions are explicitly evaluated. 

Consider for example the quark 6-point function G*-^-*: is defined as the vacuum 
expectation value of time-ordered quark and antiquark field operators, 

3 

G'^^)(a;i,a;2,X3;2/i,y2,y3)aia2a3;/3i/32/33 = (0| 7" J| qa,{xi) qfS.iVi) |0) , (2.28) 

1=1 

where the Euclidean momentum-space representation is given by 



j=i 



(2.29) 



/33 



Greek subscripts collect Dirac, color and flavor indices of the quark legs. A (5-function 
was extracted in the Fourier transform due to translation invariance and thus total 
momentum conservation. The dressed propagator of a single quark i is given by 

{27:)H\h - Pi) S{h)^fs := / d% [ d\ e^(^— G^^\x,, yi)^p (2.30) 



(3) 

such that Gq , the disconnected product of three single propagators, is written as 



{27:)H^ [Y^Ah - Pi)) Gf\ki, . . . ,P3)ai...A = ll{27T)^S\ki - k) S{ki)a,p^ . (2.31) 

i=l 



By virtue of the delta functions, the 6-point function effectively depends on 5 momenta, 
the quark propagator on 1 momentum and the three-propagator product on 3 momenta. 
In the three-body case, the pole assumption (2.21) is written as: 

^(3)/> ^ P2^-Af2 o,^ ^(fcl,fc2,fc3)aia2a3^(Pl^P2,P3)/3i/32fe Q^^ 
^(A;i,...,P3jai.../33 ' 2^1^ p2 _^ ' '^^■^^) 

where the three-quark 'wave function' <I> := Gq'^''^' corresponding to the three-quark 
amplitude ^ is defined to be the transition matrix element between the vacuum and 
the bound state with momentum P = "^^Pi = kf. 



(2vr)^(5^(X;iPi - P) '^{Pl,P2,P3)aia2a3 - 



3 

JJ / d'^Xi e'P'"'' (0| q{xi)a, q{x2)a2 q{x3)a3 \P) 

i=l •' 



(2.33) 



The products in, e.g., Eq. (2.15) are understood as summations over all occurring 
indices as well as integrations over all internal momenta. For instance, the expression 
appearing in Eq. (2.17), 

^(3)^(3) ^ i^(3)G.(3)^(3)^(3) ^ (2.34) 
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involves nine 4-dimensional integrations in tlie first place. Since aside from total mo- 
mentum conservation also the total momenta of K^^^ and G'q^ are conserved on their 
own, the original nine integrals are reduced to six. Four of those are canceled by the 
two delta functions in every occurrence of such that two integrals (thus integra- 
tion over 8 real variables) according to the two internal loop momenta remain in the 
final expression. 



Chapter 3 

Mesons 



The simplest bound states in QCD are those composed of a quark and an antiquark. 
The corresponding Bethe-Salpeter equation has been formulated in [73] and relies upon 
non-perturbative expressions for the involved quark propagator and qq kernel. First 
reliable numerical solutions were obtained in the context of a rainbow-ladder truncation 
[74-78], i.e. an iterated vector- vector gluon exchange between quark and antiquark, 
which has become a standard approach since then. 

A variety of results related to meson spectroscopy [70,78-81], electromagnetic prop- 
erties [82-86], decays [82,87,88] and scattering processes [89,90] have been reported in 
such a setup; see Refs. [19,91] for an overview. The overall tendencies may be summa- 
rized as: (isovector) pseudoscalar and vector-meson ground state properties agree well 
with experiment starting from the chiral limit up to the bottom quark; radial excita- 
tions are sensitive to the model parameters in the interaction; axial-vector mesons are 
too light compared to experimental data. 

The apparent reason for the popularity of rainbow-ladder is tied to its nature of 
being the simplest truncation of the qq kernel that implements the correct scheme of 
chiral symmetry and its spontaneous breaking. Quark-antiquark interactions beyond a 
simple gluon exchange must be consistent with the truncation of the quark propagator 
to maintain the GMOR relation for the pion. Efforts to go beyond rainbow-ladder 
have been made, and are underway, but typically suffer from a drastic amplification in 
complexity. Phenomenologically important corrections beyond RL involve pseudoscalar 
meson-cloud effects: the attractive nature of the pion cloud should induce a sizeable 
decrease in the vector-meson mass toward the chiral limit and similarly affect related 
observables. 

In the current chapter we introduce the basic relations which will frequently be 
referred to in subsequent parts of this thesis. Those are: the quark DSE, meson BSE, 
and the effective coupling a{k'^) which provides the common link in our studies of 
99) QQi QQQ and q{qq) systems. We will investigate the properties of pseudoscalar and 
vector mesons using different inputs for the effective coupling, compare their current- 
mass evolution to lattice data, discuss pionic corrections, and extract simple relations 
which can be used to describe the results. 



3.1 Quark DSE and meson BSE 
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3.1 Quark DSE and meson BSE 



Quark propagator. The basic quantity which appears in any of the bound-state 
equations (2.23) is the renormahzed dressed quark propagator S{p, fi). It consists of two 
dressing functions and as which correspond to a general fermion propagator's vector 
and scalar Lorentz structures. They can be expressed via the quark renormalization 
function Zf and the quark mass function M, where the latter is independent of the 
renormalization point /i: 

Sip,^) = -z^a,{p^^') + as{p',f,') = pf^^^(^2)2 {-^^ + M{p')) . (3.1) 

Another frequently used notation involves the functions A{p'^,ij?) = 1/ Z f{p'^ , n^) and 
B{p\^^) = M{p^)/Z^{p^^J:'). 

These dressing functions represent the solution of the quark Dyson-Schwinger equa- 
tion, also known as the QCD gap equation (see Fig. 3.1) 

S-~\p,^i) = Z2{fi^A^) (i|^ + mo(A2)) +S(p,/i,A) , (3.2) 

where Z2(/i^,A^) = A{p'^ = A^,//^) is the quark renormalization constant and A an 
ultra-violet regularization scale, mo is the cutoff-dependent bare current-quark mass. 
The quark self-energy S is defined via 

^ A 

S(p,^,A) = -^Zi^(/.2,A2) J i^>^S{q,fi)D'^'^{k,fi)r'^{l,k,^i) , (3.3) 

q 

where the prefactor (A'^^ — l)/(2A'^c') = 4/3 stems from the color trace. S involves 
the gluon propagator D^^'^ with gluon momentum k = q — p, and a bare {g Z\p i^f^) 
and dressed (gT^) quark-gluon vertex with renormalization constant Zip, where we 
introduced the average momentum / = (q + p)/2. These Green functions either need to 
be known a priori or determined in the course of a self-consistent solution of the DSEs 
of the quark and gluon propagators together with the quark-gluon vertex. 

Technical details of the solution of Eq. (3.2) are sketched in App. A.l. We note 
that the current-quark mass dependence which appears in Eq. (3.2) in terms of the 
bare mass mo (A) can alternatively be expressed via the mass function M(/i^) at the 
renormalization point through Eq. (A.9). The asymptotic form of M(p^) is given in 
Eq. (A. 3) and defines the renormalization-point independent current mass m. Given 
a sufficiently large renormalization point /i, m can be determined from M{fj?) via 
one-loop evolution, cf. Eq. (A. 5), whereas in the chiral limit: m = 0. 

Gluon propagator and quark-gluon vertex. The dressed gluon propagator, char- 
acterized by a dressing function Z, is in Landau gauge transverse to the gluon momen- 
tum k = q — p: 
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Figure 3.1: The quark DSE (3.2) in pictorial form. 
The dressed quark- gluon vertex consists of 12 tensor structures and can be written as 

4 

i=l 



T^^{i,k,f,) = ^(fl'ur + fPi^ + frk^^) n{i,k), (3.5) 



where the f^^\P,l-k,k'^ , fj,'^) are Lorentz-invariant dressing functions. A possible rep- 
resentation of the Dirac basis elements is given by 

n{l,k) = {l,^,f,[f,^]}. (3.6) 

The four longitudinal basis elements ~ k^ do not survive in the quark-DSE integral 
because of the transversality of the gluon propagator. Likewise, only the transverse 
projections of the remaining ones provide a non- vanishing contribution. In accordance 
with the notation of the quark propagator's dressing functions, the two covariants i'j^ 
and P are referred to as the vector and scalar components, respectively. 

~ ~ 1/2 ~ 

Using the STIs in Landau gauge. Zip = Z2/Z3 and Zg =1, where Z3, Z3 

and Zg are ghost, gluon and charge renormalization constants, the quark self-energy 
integral of Eq. (3.3) becomes 

E(p,^,A) = -^Z|| z^M5(g,;,)^^(af)i7'' + afr) n{l,k), (3.7) 



9 



1=1 



(i) 

where we defined the coefficients af as combinations of the gluon dressing function 
and the vertex dressings: 

a^\lM-k,k') = f-^Z{k',f^')f^'\lM.k,k^t,'). (3.8) 
47r Z2Z3 

1/2 

They are independent of the renormalization point, as can be inferred from Zg Z3 Z^ = 
1 and the renormalization-scale dependence of the quantities g ^ l/-^g) -2^^ l/-^3 and 
fi ~ Z2/Z3. 



Solution of a coupled DSE system. Both gluon propagator and quark-gluon vertex 
satisfy their own DSEs. Progress on a consistent solution of this system of DSEs has 
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Figure 3.2: The meson BSE (3.9) in RL truncation, Eq. (3.11). 



been made both analytically in terms of infrared exponents of the Green functions, e.g. 
[41,92-96], as well as numerically for general momenta in certain truncations [42,97,98]. 
The implementation of such a scheme is however beyond the scope of the present study. 
In the following we will motivate the use of underlying symmetry properties of QCD 
to arrive at a consistent truncation whose numerical implementation is feasible in our 
bound-state approach. 

Meson BSE. The quark-antiquark bound state amplitude T(p,P) with relative mo- 
mentum p, total momentum P, and mass M (at = —M"^) is the solution of a 
homogeneous Bethe-Salpeter equation (see Fig. 3.2) 



where Greek indices represent Dirac, color and flavor indices, and the quark and an- 
tiquark loop momenta are qj^ = q + aP and q- = —q -|- (1 — (t)P. The conventions 
for the momenta were chosen to comply with the diquark case, cf. Eq. (5.6). The 
momentum partitioning parameter a G [0, 1] is arbitrary since in a covariant descrip- 
tion there is no frame-independent definition of a relative momentum. Translation 
invariance implies that for each BSE solution T{p, P; a) a family of solutions of the 
form r(p + {a — <j')P, P; a') exists [99]. For equal quark and antiquark masses a value 
£7 = 1/2 maximizes the calculable meson mass (see App. B.3) and simplifies its bound- 
state amplitude. 

The arguments of the quark propagator's dressing functions in (3.9) are complex 
for timelike P^; methods to evaluate the quark propagator in the complex plane of 
Euclidean four-momentum-squared are discussed in App. A.l. The amplitude's depen- 
dence on p, P can be formulated in terms of the three Lorentz invariants p'^, P^, and 
z = p ■ P. Calculations are simplified if the dependence on the angular variable z is 
expanded in Chebyshev polynomials (see App. B.2). The spin structure of the ampli- 
tudes for different sets of quantum numbers and the method for solving the BSE are 
presented in App. A. 3. 



A 




(3.9) 



g 
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Figure 3.3: The axial-vector Ward-Takahashi identity (3.10) relates quark self-energy 
and quark-antiquark kernel. Crossed circles denote a 7^ insertion. 



The kernel K is the amputated quark-antiquark scattering kernel which is irreducible 
with respect to a pair of qq lines. Along with the quark propagator, it provides the 
physical input to the meson BSE and must be known in advance to obtain a solution 
for the meson's amplitude and mass. 

Rainbow- ladder truncation. The central identity which ensures the correct imple- 
mentation of chiral symmetry and its dynamical breaking in a bound-state approach is 
the AVWTI [68]. It provides a relation between the quark self-energy and the quark- 
antiquark kernel, the latter of which appears in the meson's bound-state equation. The 
identity can be expressed as 

1^5 + ^5|^^ ^ _ I K^^^^^ip, q, P) S{-q_) + S{q+) 7^}^, (3.10) 

and is sketched in Fig. 3.3. A qq kernel which preserves the AVWTI ensures a massless 
pion in the chiral limit as the Goldstone boson related to dynamical chiral symmetry 
breaking. In addition, Eq. (3.10) leads to a generalization of the Gell-Mann-Oakes- 
Renner relation for all pseudoscalar mesons and current-quark masses [68]. In this 
respect it is imperative for any meaningful truncation of the system of DSEs and BSEs 
to satisfy this identity. 

A systematic procedure to formulate a qq kernel which preserves the AVWTI through 
functional derivatives of the quark self-energy has been introduced in [100]. Following 
this prescription, several such constructions have been devised in the literature [42, 
101-110]. 

The simplest setup which corresponds to the lowest order in such a symmetry- 
preserving truncation scheme is the rainbow-ladder (RL) truncation. In this framework 
the qq kernel is expressed by a gluon ladder exchange, including the gluon propagator, 
one bare and one 'dressed' quark-gluon vertex. To satisfy Eq. (3.10), the dressed vertex 
may however only involve Dirac basis tensors with an odd number of gamma matrices, 
and it can only depend on the gluon momentum k. This leaves a vector part 7^^ with 
a purely fc^-dependent vertex dressing the only option in both quark DSE and meson 
BSE. The resulting ladder kernel is written as 
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where the SU(3)c Geh-Mann matrices Aj are exphcitly stated. The flavor structure in 
the BSE gives no contribution in an equal-mass system with isospin symmetry. 

In the same way as ear her, the gluon propagator and quark-gluon vertex dressings 
have been combined into an efl'ective coupling a{k'^). By virtue of the RL truncation, 
it is related to the of Eq. (3.8) via 

^(/^ / • k, k"^) = a{k'^) , and all other Jf^ = . (3.12) 
It poses the single unknown function in the rainbow-ladder approach. 



3.2 Effective quark-gluon interaction 

By virtue of the RL truncation, the entire framework rests upon a choice for the effective 
coupling a{k'^). Rainbow-ladder represents the perturbative remainder of both quark- 
gluon vertex and the qq kernel. To satisfy the one-loop relations of perturbative QCD, 
aik"^) must approach the asymptotic behavior of QCD's running coupling: 

"<^- > • STWAJ^ ' 

where 7m = 12/(llA'^c' ~ 2A^j) is the anomalous dimension of the quark propagator (in 
our calculation we use 7m = 12/25 which corresponds to Nf = 4). On the other hand, 
the interaction should exhibit sufficient strength at small gluon momenta to enable 
dynamical chiral symmetry breaking and the generation of a constituent-mass scale for 
the quark — a feature which would be the result of a combined DSE solution. This 
translates into strong non-perturbative enhancements of the quark dressing functions 
A{p'^) and MijP') at infrared momenta, see, e.g., [17]. 

Several models for a(/c^) combining the UV limit with an ansatz in the infrared 
have been employed in the past and applied to detailed studies of meson physics [75, 
76, 78-80, 111, 112]. In the present study we implement the interaction of Maris and 
Tandy [79] which reads 



cvr 



2 




lnJe2-l+(l + fcVA|^^ 



«(fc^) = ^ ( p ) e-^V(-^Ag) + , (3.14) 



where we use Aqcd = 0.234 GeV and Aq = 1 GeV. The first term of (3.14) charac- 
terizes the infrared properties, expressed through the two parameters c and oj which 
will be discussed in detail below. It provides the characteristic infrared strength which 
is crucial for a dynamical quark mass generation. The second term accounts for the 
ultraviolet behavior of (3.13) and is thus constrained by perturbative QCD. 

We note that the ansatz (3.14) behaves as a{k'^) fe^ for k'^ — > 0. This facilitates the 
numerics as the self-energy integral is of lesser divergence in the infrared, see App. A.l. 
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Nevertheless, since the quark DSE is dominated by the interaction at intermediate 
momenta, the exphcit behavior in the deep infrared does not play an important role for 
hadronic ground states. This has been explicitly verified in [108] with a parametrization 
different from (3.14). 

Physical input. The two parameters c and to which appear in the ansatz (3.14) mod- 
ulate the coupling's strength and width in the infrared. Studies using this parametriza- 
tion have demonstrated that pseudoscalar and vector-meson ground state properties 
are insensitive to a variation of the coupling width u; in a certain range [79], i.e. inde- 
pendent of the detailed shape of the interaction in the infrared. Such an insensitivity 
has been found for the ground state of the nucleon in the quark-diquark model as 
well [113]. To highlight this property we will present subsequent results in terms of 
'w-bands' which denote a variation uj = u) i: Auj. 

In the context of the aforementioned hadron masses, the infrared coupling strength 
c is accordingly the only active physical parameter in the entire setup. In contem- 
porary meson studies it has usually been fixed to reproduce the experimental value of 

= 131 MeV at the u/d current-quark mass. Guided by results from experiments and 
lattice QCD, one may choose a current-mass dependent coupling strength c(m), where 
m is the renormalization-point independent current-quark mass, to implement phe- 
nomenologically reasonable assumptions that comply with the heavy-quark sector. In 
the following we will give a short overview of the strategies that have been undertaken 
in the literature and will be used as an input to our calculations. 

(CI) Current-mass independent coupling strength. A setup which has been ex- 
tensively used in the literature [70,71,78,79,81,83-88,91,114-116] involves a coupling 
strength that is independent of the quark mass, chosen to reproduce the phenomeno- 
logical quark condensate and experimental pion decay constant /tt = 131 MeV at the 
u/d current-quark mass associated with m-,^ = 138 MeV. The corresponding value of 
the coupling strength is c = 0.37. It enables a reasonable description of masses, de- 
cay constants and electromagnetic properties of ground-state pseudoscalar and vector 
mesons up to bottomonium [115]. 

(C2) Fit to the lattice mass function. The availability of lattice data for the quark 
propagator, gluon propagator and quark-gluon vertex, e.g. [117-124], has provided a 
means to test the properties of Dyson-Schwinger solutions, in particular with regard 
to their larger quark-mass behavior. Several parametrizations for the quark-gluon 
interaction have been employed in the literature and compared to lattice results, either 
in a rainbow-ladder-like context [46, 125-127] or more general setups [42,45, 104, 108, 
128]. 

A strategy explored in Ref. [129] was to adopt the ansatz (3.14) in the rainbow- 
ladder approach, where its strength c(m) is kept intact in the chiral region to recover 
/tt at the physical pion mass, but diminished for larger masses such that the resulting 
quark mass function M(p^) agree with quenched lattice results of Ref. [117]. A current- 
mass dependent coupling therefore emulates to some extent a quark-mass dependent 



3.2 Effective quark-gluon interaction 



27 



structure in the quark-gluon vertex beyond rainbow-ladder. Upon solving the respec- 
tive meson and quark-diquark BSEs, corresponding results underestimate nip and Mjsr 
obtained from lattice calculations at larger quark masses, e.g., by ~ 10% at the strange- 
quark mass [129]. This result either suggests inconsistencies between different lattice 
techniques or an inaccuracy of rainbow-ladder well beyond the strange-quark mass. 

(C3) Fit to the /j-meson "quark core". A third practical strategy is to implement 
phenomenological assumptions about the nature and possible impact of corrections 
beyond RL upon hadron properties, and to adjust the current-mass dependence of the 
effective coupling in a way where these corrections are explicitly missing. 

One important effect in the chiral and low-energy regime of QCD is imposed by 
pseudoscalar meson-cloud contributions, cf. Section 2.2. Such corrections provide a 
substantial attractive contribution to the 'quark core' of hadronic observables in the 
chiral regime whereas they vanish with increasing current-quark mass. Their impact 
on the chiral structure of the quark mass function and condensate, /tt, rrip, and nucleon 
and A observables has been demonstrated in the NJL-model [136,137], DSE studies 
[61,64, 108, 138], and chiral extrapolations of lattice results [58]. A sizeable reduction 
of p-meson, nucleon and A masses is paired with an increase of hadronic charge radii 
towards a logarithmic divergence in the chiral limit. The latter effect is clearly missing 
in a RL truncation as can be inferred, e.g., from the BSE result for the pion charge 
radius displayed in Fig. 3.7. 

A further resummation of non-resonant Abelian diagrams in the quark-gluon ver- 
tex and qq kernel provides additional attraction in the vector-meson channel which 
decreases with increasing current-quark mass [103, 105, 107]. On the other hand, an in- 
clusion of the three-gluon vertex exhibits a substantial amount of repulsion [110] which 
suggests a non-perturbative cancellation mechanism beyond RL [102]. Nonetheless one 
may construct a rainbow-ladder 'quark core' which overestimates the experimental p- 
meson mass, most noticeably towards the chiral limit, and resembles the hadronic quark 
core of chiral effective field theories which is subsequently dressed by chiral corrections. 
Such an inflated quark core mass for rrip has been used in [60] via (see Fig. 3.5) 



X 



2 



I + xl/{^.Q + xl) , Xp:=mp/m^ := (3.15) 



with the chiral-limit value = 0.99 GeV. The sum of corrections beyond RL would 
then reduce nip in the chiral limit by ~ 25% whereas the quark core contribution to 
nip approaches lattice results above the s-quark mass. To reproduce Eq. (3.15) upon 
solving the p-meson BSE, the coupling strength c of Eq. (3.14) must be equipped with 
the following current-mass dependence (see Fig. 3.4): 



0.86 6( Aw) 
1 + 0.885 Xq + (0.474 Xg)^ ' 



c(u;, m) = 0.11 + ^ , ^ ^ , ^ ^y , Xg := m/(0.12GeV). (3.16) 



At each value of the current-quark mass ifi, the parametrization 



6(Aw) = 1 - 0.15 Aw (1.50 Aojf + (2.95 Aujf 



(3.17) 
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Figure 3.4: Left panel: Coupling strength c(a;,m) in setup (CI) (dashed line) and 
(C3) {solid band), cf. Eq. (3.16). The current-mass evolution is expressed in terms 
of the squared pion mass obtained from the pseudoscalar-meson BSE. Right panel: 
Coupling strength in (C3) as a function of the width parameter uj for five different 
current-quark masses. The current-mass dependent to plateaus are reflected in (3.18). 
The dotted vertical line corresponds to Aa; = 0, the shaded region to \A.uj\ < 0.1. 




[GeV^] 



Figure 3.5: Solution for mpim^) in setup (CI) [dashed line) and (C3) (solid line), 
where the latter is identical to the parametrization of Eq. (3.15). The curves may be 
regarded as the input which defines the coupling strength c(rn) in both setups. We 
compare to a selection of lattice data [130-132] extracted from Ref. [133], and the 
results of [134] together with a chiral extrapolation [135] (blue band). 
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eliminates the residual to dependence of the BSE result for rup and ensures the validity 
of Eq. (3.15) in the range u) = u){m) it |Aa;|, with the central value given by 

w(m) = 0.38 + 0.17/(1 + |Au;| < 0.1 . (3.18) 

As demonstrated in Refs. [60, 113] and this thesis, the procedure induces consistent 
overestimated values for a range of vr, p, nucleon and A observables as obtained from 
their respective meson BSEs, quark-diquark BSEs and Faddeev equation. The results 
are only weakly w- independent within the range (3.18). 

By implementing either of the discussed models for the current-mass dependent 
coupling strength, the combination of quark DSE (3.2), meson BSE (3.9), rainbow- 
ladder kernel (3.11), and effective coupling (3.14) completely determines the subsequent 
results. In the course of this thesis, these results will be primarily presented in the 
'quark core' model (C3) and occasionally compared to those obtained with a fixed 
coupling strength (CI). Shaded bands denote the sensitivity to a variation of oj which 
we consider for the (C3) setup only; all results corresponding to (CI) are plotted with 
a current-mass independent value u) = 0.4. A connection between the two setups will 
be established in the next section. 

3.3 Analysis of rainbow-ladder meson results 

Quark propagator. A characteristic feature of QCD is the dynamically generated 
enhancement of the quark mass function M{p^) at small momenta, visible in Fig. 3.6. 
The resulting constituent-mass scale, e.g. M{p'^ = 0), is typically several hundred MeV 
larger than the current-quark mass which is the input of the DSE. In addition. Fig. 3.6 
illustrates the impact of the overestimated quark-core model (C3) on the quark prop- 
agator dressing functions M(p^) and Zf{p'^). The resulting mass function at infrared 
momenta is considerably larger than the respective lattice results. An implementation 
of pionic effects in the qq kernel reduces this difference [108, 128]. 

The singularity structure of the resulting quark propagator, i.e., of its denominator 
1 / {p^ + M^(p^)) , is that of complex conjugate poles, a feature which may be an artifact 
of the rainbow truncation but has been found in more general truncations beyond 
rainbow-ladder as well [42, 128]. The pole positions depend on the infrared width uj, 
see Fig. 3.6: a larger value uj forces the poles closer to the timelike axis while a smaller 
value shifts them further into the complex plane. The resulting trajectory is similar to 
the boundary of the parabolic integration domain which is needed in the subsequent 
bound-state equations. It constitutes the singularity limitations encountered in the 
calculation of form factor diagrams, cf. App. B.3. 

The vacuum quark condensate can be calculated from Eqs. (A.3-A.5). The value 
shown in Table 3.1 is obtained from the perturbative tail of the chiral-limit mass func- 
tion and evolved to the scale /i = 1 GeV via the one-loop formula (A. 5). As discussed in 
connection with Eq. (3.21), the quark-core setup (C3) operates with an inflated scale 
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Figure 3.6: Left panel: Quark mass function M{p'^) and wave-function renormaliza- 
tion Zf{p'^) in setup (C3), at a current-mass value m which corresponds to rrin = 0.63 
GeV, compared to lattice results [120]. A renormalization point fj, = 2.9 GeV was 
chosen. The lower (upper) edge of the bands is related to the largest (smallest) value 
of u) for M{p'^), and vice versa for Zf{p'^). Right panel: Complex conjugate poles of 
the quark propagator in the complex p'^ plane for different values of uj. The apex of 
the corresponding parabola (B.17) defines the w-dependent 'pole mass' m-g. 
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Table 3.1: Comparison of the quark condensate, vr and p decay constants, /)-meson 
mass, and pion charge radius in setups (CI) and (C3), characterized by the parameters 
Air, and rj defined in (3.19). A variation of rj = 1.8 it 0.2 in (C3) coincides with 
a; ~ u) lb 0.06. The results correspond to a current mass m = 6.1 MeV which is related 
to the physical pion mass ttItt = 138 MeV. Experimental or phenomenological values 
are quoted in the first row. In the first three rows, ry is dimensionless and r^^ is given 
in fm, all other units are GeV. The last row plots the ratios of sets (CI) and (C3). 
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Figure 3.7: Pion-mass evolution of pion decay constant {upper panel) and squared 
pion charge radius {lower panel) in setup (CI) {dashed lines) and (C3) {shaded bands). 
The compilation of lattice results for [130,131,139-144] was extracted from Ref. [133]. 
Lattice data for are from Refs. [145-147] . 
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of dynamical chiral symmetry breaking. This overestimation compared to the phe- 
nomenological value is visible in the quark condensate and consistent with the solution 
for the quark mass function. Pionic contributions to the condensate are attractive, cf. 
Refs. [108,136]. 

Meson properties. A solution of the pseudoscalar and vector-meson BSEs provides 
the respective bound-state amplitudes which are subsequently used to calculate physical 
meson properties. Table 3.1 collects results for vr and p decay constants, obtained 
from Eqs. (A.44), and the pion charge radius, Eq. (A.48). One observes a uniform 
response to the inflated quark core for vrip: expressed in units of mass, the tabulated 
quantities related to the input (C3) consistently overestimate their experimental values 
by ~ 30 — 35%. Moreover, fj^ and r^r tend to approach lattice results for heavier quarks 
(see Fig. 3.7). This validates the notion of a pseudoscalar meson cloud which increases 
a hadrons charge distribution towards the chiral limit where the charge radius would 
diverge. 

Relation between the models. Simple relations between the setups (CI), (C2) 
and (C3) can be established in the chiral limit. To this end it is beneficial to rewrite 
the infrared part of the coupling (3.14) in a more suggestive way. Upon replacing the 
infrared parameters c, uj and Aq by two new parameters Air and r/, defined via 

the infrared contribution to the effective coupling a^k"^) is expressed through 

aij^(fc2) =7r7?Ve-'''^ x = A;VA?r. (3.20) 

Air only appears in the denominator of x, hence it is the only dimensionful scale of the 
rainbow-ladder truncated DSE-BSE system in the chiral limit as long as the UV part 
of the coupling is not taken into account. It represents the scale of dynamical chiral 
symmetry breaking which, in a coupled solution of quark, gluon and ghost DSEs, would 
be generated self-consistently. As a consequence, mass-dimensionful quantities which 
are sensitive to the infrared properties scale with Air. The chiral- limit values of Air 
are 

(CI): Air = 0.72 GeV, (C3): Air = 0.98 GeV, (3.21) 

which makes clear that the model (CI) for a given set of observables will, upon entering 
its 'core version' (C3), produce results which are overestimated by the same percent- 
age. Dimensionless chiral-limit ratios between the two setups will be equal, namely 
(0.72)7(0.98) = 0.73. Table 3.1 demonstrates that this scaling property still persists at 
the physical u/d mass where the disturbance from a non- vanishing current-quark mass 
is negligible. 

The second parameter replaces the coupling width oj. An insensitivity of ob- 
servables with respect to a variation of the width u at a certain coupling strength c 
translates into an invariance with respect to jy at a fixed scale Air. The combined 
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change of c and oj from (CI) to (C3) according to Eqs. (3.16-3.17) is equivalent to a 
rescaling of Air where is essentially unchanged. 

A rainbow-ladder 'mass formula'. The analysis can be extended to finite current- 
quark masses where the additional scale rh or, equivalently, the pion mass m-,^ is in- 
troduced. An investigation of the numerical results exhibits that quantities with the 
dimension of a squared mass, generically denoted by U , allow for a scale separation in 
the following way (see Fig. 3.8): 

Ml{ml) ^ afjiv) A?R + bfjiml) ml . (3.22) 

For a vanishing pion mass, the first term in (3.22) reproduces the observation from 
above: masses scale with Air, and the dimensionless chiral-limit values au may depend 
on the infrared width parameter rj. For a finite current-quark mass, the second term 
is found to be insensitive with regard to the infrared parameters Air and rj: it only 
depends on itItj and the remaining scales Aqqd and Aq in the coupling. 

The dimensionless coefficients bjj are generic functions of which account for any 
additional pion-mass dependence of the results. In the heavy-quark limit they would 
represent the dimensionless values of Mjj in units of mT^. A BSE analysis of heavy- 
meson observables [116, 148, 149] confirms that meson masses become proportional to 
m^, such that e.g. bmp{oo) 1, whereas decay constants behave as ~ ^/rn^ which 
implies that bf^ and bf vanish with an inverse square root of the pion mass. 

We tested Eq. (3.22) up to m 200 MeV for a range of quantities 

U = M(0), ruq, nip, fn, fp, l/r^r, rrisc, mav, Mn, Ma, (3.23) 

where the diquark and baryon masses are determined in Chapter 5. The relation 
breaks down below Air < 0.5 GeV which is related to the threshold of dynamical 
chiral symmetry breaking in the quark DSE. In the domain of its validity, a rescaling 
of Air — > Ajpj^ in the infrared part of the effective coupling (3.14) produces a current- 
mass independent additive contribution to squared masses. This infrared component 
may however depend on the parameter r]. In the 'core model' (C3), Ajj^ was chosen 
pion-mass dependent; nonetheless Eq. (3.22) provides an estimate for its results: 

Mi\ml) « Ml,{ml) + al{rj) [A'^mD - A^r) . (3.24) 

Moreover, if two quantities U and V independently satisfy au « bu and ay « 6y, 
then their ratio Mu /My ~ au jay will only weakly depend on and Air. This is most 
remarkably realized for the dimensionless product fT^r-,^ (Fig. 3.9): for the pion-mass 
range under consideration, it is constant in m^, identical in (CI) and (C3) (i.e., inde- 
pendent of Air), insensitive to i] (since both /,r and r^r are very weakly ry-dependent, 
see Fig. 3.7), and it agrees with the experimental value and lattice-QCD results [137]. 
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Figure 3.8: Left panel: and the difference m^p—a^p Ajj^ for both model versions (CI) 
and (C3), compared to m^. Right panel: The values of bu extracted from Eq. (3.22) 
for the four calculated quantities /tt, /p, mq and nip. The meaning of the quark pole 
mass mq = Im\/??polc clarified in Fig. 3.6. Shaded bands: (C3); dashed curves: 
(CI) for fixed r] = 1.8. The horizontal lines (or shaded areas, respectively) depict the 
chiral-limit values au{r]) for comparison. Qualitatively similar features are found for 
diquark and baryon masses. 




Figure 3.9: Left panel: The dimensionless product f-j^r^^ as a function of m^, compared 
to the experimental value and lattice results; figure taken from [137] . Right panel: Pion- 
mass evolution of several ratios, normalized by their chiral-limit value. 



Chapter 4 

Baryons: Three-body equation 



The three-body bound-state problem has a longstanding history which dates back to 
the original work by Faddeev [150]. Non-relativistic Faddeev equations have found 
widespread application in the description of three-nucleon systems, see Ref. [151] for 
an overview. The covariant generalization of the Faddeev equation to the three-body 
analogue of a Bethe-Salpeter equation was formulated in Refs. [65, 152]; a comprehen- 
sive introduction can be found in [67]. Within the framework of Section (2.3), the 
equation describes the baryon as a bound state of three spin-i/2 valence quarks where 
the interaction kernel comprises two- and three-quark contributions. 

A solution of the covariant three-body equation requires knowledge of the dressed 
quark propagator and the three-quark kernel; and a specification of the Poincare- 
covariant baryon amplitude. The relativistic spin structure of the latter has been 
explored in [153,154] and described in the light-front formalism in [155-158]. A com- 
plete classification according to the Lorentz group and the permutation group S3 was 
derived in [159] in terms of covariant three-spinors. Their analogues in the form of 
Dirac tensors, kindred to the decomposition of Green functions and meson amplitudes 
encountered in previous chapters, will be stated below (cf. Table 4.1). 

The complexity of the three-body bound state equation has so far prevented a direct 
numerical solution. Upon implementing perturbative quark propagators it has been 
studied, for instance, in the works of Refs. [160, 161], in the context of a three-body 
spectator approximation [162], or a Salpeter equation with instantaneous forces [67]. 
The corresponding equation of a scalar three-particle system with scalar two-body 
exchange based on the Wick-Cutkosky model [163, 164] was recently investigated and 
compared to the light-front approach [165]. An appealing strategy to simplify the 
three-quark problem while maintaining full Poincare covariance is provided by the 
quark-diquark model which will be discussed in Chapter 5. 

The present chapter is devoted to a novel solution of the three-quark bound-state 
equation of the nucleon where the full Dirac structure of the covariant amplitude is 
taken into account. The numerical computation is performed upon truncating the 
kernel to a rainbow-ladder gluon exchange which allows for a direct implementation 
of the effective quark-gluon coupling a{k'^) introduced in Chapter 3. The resulting 
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current-mass evolution of the nucleon mass, when compared to lattice data, exhibits a 
behavior which is qualitatively similar to that of the vector meson in the same approach. 

A future extension to more sophisticated interaction kernels is certainly possible. 
This involves an inclusion of irreducible three-body forces, supposedly dominated by 
a three-gluon coupling to any of the three quark lines, and a generalization towards 
a quark-quark interaction beyond rainbow-ladder. In the context of meson studies, 
the latter has attracted a considerable amount of attention in recent years, and an 
implementation of the findings in the three-body approach will certainly provide further 
insights into the physics of baryons. 

4.1 Faddeev amplitude and equation 

Baryons appear as poles in the three-quark scattering matrix. The derivation pursued 
in Section 2.3 leads to a relativistic three-body bound-state equation which is the direct 
analogue of the Bethe-Salpeter equation in the quark-antiquark channel: 

* = ^(3) ^ ^ ^(3) ^ ^(3) + Kf^ , (4.1) 

1=1 

where ^' is the baryon's bound-state amplitude. The respective 3-body kernel ET^^^, 
stated in Eq. (2.20), comprises a 3-quark irreducible contribution and the sum of per- 
muted two-quark kernels, where the subscript i denotes the respective accompanying 
spectator quark. The quark-antiquark analogue of the two-body kernel X^^-* appears in 
the meson BSE (3.9) where we employed a ladder truncation together with a rainbow- 
truncated quark propagator. 

It has been the guiding assumption for the quark-diquark model that correlations of 
two quarks provide the dominant binding structure in baryons. This was inspired by 
noticing that colored two-quark states can appear in an SU{2>)c anti-triplet or sextet 
configuration where the former, in combination with a color-triplet quark, allows for 
the formation of a color-singlet nucleon. The observation motivates the omission of 
the three-body irreducible contribution from the full 3-body kernel. The resulting 
relativistic Faddeev equation includes a permuted sum of two-body qq kernels: 









(4.2) 


i=l 





which enables an adoption of the formalism established in Chapter 3 in the three- 
body problem. Eq. (4.2) represents a technical simplification as well since the equation 
merely includes one momentum loop and becomes computationally tractable. 



Nucleon amplitude. The bound-state amplitude "^ap-fS of a nucleon carries 3 spinor 
indices 7} for the involved valence quarks and one index 5 for the spin-1/2 nu- 

cleon. It depends on three quark momenta p\, p2, Ps which may be reexpressed in 
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terms of the total momentum P and two relative Jacobi momenta p and q. They are 
related via: 

P={^-V)P3-VPd, Pi = -9+ y = -g- 2 2~ ' 

P2-P1 , Pd P , '^-V r> (4.3) 

g=^^, P2 = q+-=q-^+—P, ^ ^ 

P =Pl+P2+P3, P3=p + 7]P , 

where we abbreviated pd '■= Pi + P2- The relations between the momenta can also 
be inferred from Fig. 4.1. We have chosen equal momentum partitioning 1/2 for the 
relative momentum q and use the value r/ = 1/3 in connection with the momentum p 
which maximizes the upper boundary for the nucleon mass with respect to singularity 
restrictions (see App. B.3). 

The nucleon amplitude can be decomposed into a certain number of Dirac structures: 

64 

*a/375(P> Q^P)=Y1 fi^P^^ ^' ^)"/37<5- (4-4) 

i=l 

The amplitude dressing functions ft depend on the five Lorentz-invariant combinations 

, , ZQ=p^-q^ , zi=p- P, Z2 = q- P, (4.5) 

where a hat denotes a normalized 4-vector and pj, = Tp^p'^ a transverse projection 
^rpfiu — ^fiu _ pMpi'^^ \Yg abbreviated the angular variables by the shorthand notation 
{z} = {zq, zi, Z2}. The total momentum-squared = — is fixed since the nucleon 
is on its mass shell. The Dirac structures Ti{p,q,P) will be explained in Section 4.2. 

Faddeev equation. Using the above kinematics, the full Dirac and momentum de- 
pendence of the Faddeev equation (4.2) reads (see also Fig. 4.1): 



'^al3^5{P,q,P) = I Kaa'l3l3'{k) Sa'a"{kl) Si3'f3"{k2)'ifa"l3"'y5{P^^\q'^^\P) 

+ KfSf3>^-y'{k) Spipii{k2) Syy^ks) ^af3"-y"S (p«,<7«,P) (4.6) 
+ Ky^'aa'{k) S^i^"{k3) Sa'a"{ki) '^a"l3y"S 
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where we already anticipated a two-quark kernel K that only depends on the gluon 
momentum k which is also used as the loop-integration variable. The quark propagators 
S depend on the internal quark momenta ki = pi — k and kj = pj + k, and the internal 
relative momenta are given by 

pW=p + k, q^^^=q-k/2, 

pi^)=p-k, q(^^=q-k/2, (4.7) 

p(3) = p ^ qi3) _ _|_ ^ _ 

Specifically, we use the rainbow-ladder kernel of Eq. (3.11): 

K^aw'ik) = zl 1^ rr iL' ih' ' (4-8) 

which involves two bare quark-gluon vertices and a gluon propagator, subsumed into 
the eflFective coupling a{k'^) of Eq. (3.14). Tracing the color structure induced by the 
quark-gluon vertices leads to a color factor 2/3 in front of the integral (4.6). 



4.2 Dirac basis covariants 

A general Green function with 4 fermion legs which depends on 3 independent momenta 
involves 256 independent components. The subspaces corresponding to positive and 
negative parity consist of 128 Dirac structures each. A possible linearly independent 
basis for the positive-parity and positive-energy nucleon is given by the 64 components 

S±(p, g, P) := (r, A±) (75C) (r, A+) , 
Pg(p, g, P) := 75 (r. A±) (75C) ® 75 (r, A+) , 

where the Ti{p, q, P) with i = 1 . . . 4, given in Eq. (4.15), carry the relative-momentum 
dependence, A+(P) = (1 -|- /*)/2 is the nucleon's positive-energy projector and C = 
7^ 7^ the charge-conjugation matrix. The notation, e.g. {Sfi)ap,'y5 = (A+75C)q,/3(A+)^5, 
refers to two outgoing quark legs with indices q,/3 on the left-hand side of the tensor 
product (hence the (7^C) insertion, cf. Eq. (A. 34)) and an outgoing quark and incom- 
ing nucleon leg (indices 7, 5) on the right-hand side, where the positive-energy projector 
is attached to the latter. The ordering of indices in the tensor product is insignificant 
since, for a complete basis, any permuted version can be expressed by a sum of basis 
elements in the above "canonical" ordering using Fierz identities, see e.g. [157]. 

Each subspace of definite parity and sign of energy, corresponding to the nucleon's 
(1/2) + , (1/2)+, (l/2)~, (l/2)~ states, includes 64 covariants. The negative-parity 
basis elements are obtained by attaching a factor 7^ to either left or right-hand side 
of Eq. (4.9); the negative-energy (antibaryon) structures by replacing A_|_ on the r.h.s 
with the negative-energy projector A_ = (1 — ^) /2. The symbols S and P in Eqs. (4.9) 
were chosen to reflect the combination of two scalar or pseudoscalar covariants whose 
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product again exhibits positive parity. All possible further basis elements, e.g. 



7^75 (riA±)(75C)® 7^75 (r,A+) 
7^ (riA±)(75C)® 7^(rjA+) 

4" (r,A±)(75C7)® 4^ (r,-A+) 



(4.10) 



linearly depend on the ones of Eq. (4.9); respective relations are given in Eq. (A. 76). In 
the actual calculation we use an orthogonalized set of covariants constructed from the 
basic structures (j = 1 . . . 4) 




where each carries a superscript (it). Here we have exchanged the elements {S2j, P2j} 
by {Aij,Vij} to facilitate a comparison with the multispinor notation used in the 
literature, e.g. Ref. [159]. 

Momentum-dependent covariants. The momentum-dependent Dirac structures 
which appear in Eqs. (4.9-4.10) involve the following 8 basis elements: 



Ti{p,q,P)K±{P) G {1, [;^,^], A ^} X {1, f } 



(4.12) 



It is convenient to choose a set of momenta {py,gi,P} which are orthonormal with 
respect to the Euclidean metric, i.e. 



2 ^2 

PT = qt 



(4.13) 



This is realized via 



= q>^ = T^^ q^ = T^^ q^ , (4.14) 



where Tj^'^ = S^'^ — k^k^ denotes a transverse projector with respect to any four- 
momentum k. The four covariants Ti{p^ q, P) are those of a fermion-scalar vertex and 
effectively depend on two momenta: 

ri(p, (z, P) = {i, I itl it] ■ (4.15) 

Without loss of generality one may choose the momentum alignment (cf. Eq. (B.5)) 
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(4.16) 
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which yields in the baryon's rest frame: 
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(4.17) 



and hence h.± = ^(1 it 7'^) and Ti{p,q,P) = |l, 7^7^, 7'^, 7'^}- The loop momenta 
and covariants inside the integral of Eq. (4.6) will natm'ally have a more complicated 
appearance. 

Angular momentum decomposition. The Dirac basis elements (4.11) can be clas- 
sified with respect to their quark-spin and orbital angular momentum content in the 
nucleon's rest frame. In general only the total angular momentum j = 1/2 of the 
nucleon is Poincare-covariant while the interpretation in terms of total quark spin and 
orbital angular momentum will differ in every frame. The spin is described by the 
Pauli-Lubanski operator: 



1 



(4.18) 



where we chose the total momentum P to be normalized. J^'^ and are the generators 
of the Poincare algebra satisfying the usual commutator relations. The square of the 
Pauli-Lubanski operator, 

I.,,., .„„ ^„^„ 



j/iu jfiu _|_ pfj.pujij.ajc 



is one of the two Casimir operators of the Poincare group. Its eigenvalues are given by 
W"^ — > j{j + 1), where j describes the spin of the particle. 

For a system of three particles with total momentum P and relative momenta p 
and q, the total angular momentum operator consists of the total quark spin S and 



the relative orbital angular momentum L 
Lorentz-covariant operators 



+ L 



(</)• 



Upon subsuming them into 



(p) 



1 fiua/S pu 
4t ^ 



a 



af3 



101-^10 o-"^®l-hl0l0 d"^) 



yuuappu faQl3 _ ^(3 go 



1 (g) 1 1, 
1 (g) 1 (gi 1, 



(4.20) 
(4.21) 
(4.22) 



one may verify that the basis covariants (4.11) are indeed eigenfunctions of the square 
of the Pauli-Lubanski operator = + L|^^ -|- L^^-^ with j = 1/2. Here the following 
identities prove to be useful: 



52 = 1 1 1 1 + i (cj^'^ C7^^ 1 + perm.) , 
S-L^p) = ^pt^ {dxYp (cj^'' 1 1 + perm.) , 

= [dTYp+P^TPT (dTYidT); - p't {dT)AdT. 



P ' 



pt^ (9t)^ {dTY, -PT-qr {dT)p-{dT)g • 
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1/2 
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1/2 
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Sl3 


Pis 


Al3 
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Si4 


Pl4 


Al4 


Vl4 


8 



s 
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V6X± 


V6X% 


# 


3/2 


2 


3S33-V11 


3 P33 - An 


4 


3/2 


1 


3 S34 — 3 S43 — 2 V12 


3 P34 — 3 P43 — 2 A12 


4 


3/2 


1 


3S31-V13 


3 P31 - Ai3 


4 


3/2 


1 


3 S41 — Vi4 


3 P41 — Ai4 


4 



s 


/ 


V2X^^. 


V2X± 


# 


3/2 


2 


2 S44 + S33 — Vn 


2 P44 + P33 - An 


4 


3/2 


2 


S34 + S43 


P34 + P43 


4 


3/2 


2 


—2 S42 + S31 — Vi3 


-2 P42 + P31 - Ai3 


4 


3/2 


2 


2 S32 + S41 — Vi4 


2 P32 + P41 - Ai4 


4 



Table 4.1: Orthonormal Dirac basis constructed from Eq. (4.11) via a partial-wave 
decomposition. The superscripts it are not displayed. 



The basis states can furthermore be classified with respect to the eigenvalues of 
5^ — > s(s + 1) and — > l{l + 1), which, in the nucleon's rest frame, assume the 
interpretation of total quark spin and intrinsic quark orbital angular momentum. Such 
a partial-wave decomposition allows for an arrangement of the 64 basis covariants (32 
for total quark spin s = 1/2 and s = 3/2, respectively) into sets of 8 s waves (/ = 0), 
36 p waves {I = 1), and 20 d waves (/ = 2). We denote the resulting orthonormal basis 
states by the symbol X^., with i = 1 . . . 8 and j = 1 ... 4, and collect them in Table 4.1. 

Orthogonality. The combinations X^- satisfy the following orthogonality relation: 

\ T:{X^ X^;,,} = \ (XI,)^^ (Xi;,)„,,,. = ^n' Srr' , (4.23) 
where the charge conjugation X is defined as 

X(p, q, P) = (C C) X{-p, -q, -Pf{C^(d C^) . (4.24) 
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As an example, consider: 

Sj.(p, q,P)=C {r,(-p, -q, -P) A^ ^5Cfc^ C {T,i-p, -q, -P) A.fc^ = 

= (c^75)(A±r,) ®(A+r,), 

where we used Ti{p, q, P) = CTi{—p, —q, —P)^ C'^ and the relations A±(— P) = A=p(P) 
and A±(P) = C A±(-P) = A±(P). 

Multispinor notation. An alternative way to construct a basis for the nucleon am- 
plitude is to use the Dirac spinors U'^{P), V{P) = 75 [/"^(-P) which satisfy the free 
Dirac equation for a spin-i/2 particle, i.e. which are eigenspinors of the positive and 
negative energy projectors A-t. Normalized to U^{P) U^[P) = (5po-, they are expressed 
as 

The 64 linearly independent basis states for the nucleon wave function are obtained 
from the multispinors UUU := U®U®U , VVU, VUV, UVV and their parity-reversed 
counterparts VVV, UUV, UVU, and VUU by equipping them with the 8 possible spin- 
up/down arrangements and so on. A corresponding basis has been explicitly 
constructed in Ref. [159]. Using identities such as 

A+ = UW^ + U^U^ , A+ (75C) = U^U^ - U^U^ (4.26) 

leads to the following relations for the s-wave basis states of Table 4.1: 

-S+ W = {Wu^ - U^W^) , A+ W = {Wu^ + U^W) W -2 WWu^ , 
-s-^ W = (y - W) ui ^ _A^^ w = (y + v^V^) W -2 v^V^u^ , 

-p+ [/T = (yTf/i _ yi[/T) ^ v+ = {V^U^ + V^U^) - 2 V^U^V^ , 

-p-^ [/T = (^jTyi _ [/iyT) ^ _v-^ W = {WV^ + U^V^) -2 WV^V^ . 

Due to the Poincare-covariant construction of Eqs. (4.9-4.10), the remaining 56 covari- 
ants depend on the relative momenta. For instance, the Dirac structure S^g satisfies 
the relation 

iS^^U^ = {V^V^ -VW^) (([i5?]i + i[p?]2)W + [p?]3FT) . (4.27) 
Using the special momentum alignment (4.17) leads to 

iS^^U^ {V^V^ -V^V^)V^ (4.28) 

and similar relations for the remaining basis elements in Table 4.1. At the same time 
the expression (4.28) is the 'parity- flipped' counterpart of S^^, 

- (75 » 75) Sti ill ® U^) = {V^V^ - V^V^) , (4.29) 
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which appears in the amphtude of a negative-parity nucleon (l/2~) and that of a 
positive-parity antinucleon (1/2+) . Nevertheless the seemingly odd-parity structure 
(4.28) in the multispinor basis still contributes to the (1/2"*") state via Eq. (4.27), where 
the parity flip induced by the spinor replacement U ^ V is saturated by an odd parity 
introduced by the relative momentum pr- As a consequence, indeed all 64 three-spinor 
combinations contribute to the nucleon's amplitude. 



4.3 Quark exchange and Pauli principle 

The full Dirac-flavor-color amplitude of the nucleon reads (cf. Eq. (A. 74)): 

^{p, q, P) = {^{p, q, P)m^ ® + ^(P, q, P)ms ®^Ms}® ^ • (4-30) 

The Pauli principle requires the Faddeev amplitude to be antisymmetric under exchange 
of any two quarks, as can be inferred from Eq. (2.33). The color singlet wave function 
^ABC is totally antisymmetric, hence the Dirac-flavor amplitude must be symmetric. 
Since the two isospin-1/2 flavor tensors '^Mai ^x^, given in Eq. (A.75), are either 
mixed-symmetric or mixed-antisymmetric, the same feature must hold for the Dirac 
remainders Ma ^^"^ ^ Ms the nucleon amplitude. 
The Faddeev kernel ^(3) 

is per construction invariant under the permutation group 
: it commutes with any permutation of two quark legs. The Dirac parts of the 
solutions to the Faddeev equation can be arranged into irreducible multiplets 

("1) 

of which the first two (totally symmetric or antisymmetric solutions) are unphysical 
while the mixed-symmetry doublet constitutes the nucleon amplitude according to 
Eq. (4.30). The Faddeev equation will in general mix the two linearly independent 
solutions ^Ma ^^"^ Ms- However, since the rainbow-ladder kernel presently employed 
is flavor-independent and the two flavor tensors '^Ma ^"^^ '^Ms orthogonal to each 
other, the equations for the Dirac amplitudes '^Ma-> '^Ms decouple: 

^ = K(3)^ ^ ^Ma=K^^^Ma. (432) 

These two states do indeed emerge as independent solutions upon solving the equa- 
tion. The dominant amplitudes in either case are the mixed-antisymmetric and mixed- 
symmetric covariants defined in Eqs. (4.9) and (4.10): 

"^Ma ~ S+ = A+(75C) A+ , (4.33) 
"^Ms - A+ = 7^75 A+(75C) ® 7^75 A+. (4.34) 
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s-wave components. Upon restricting oneself to the eight momentum-independent 
covariants of Table 4.1 with s = 1/2, I = and applying the symmetrizers/anti- 
symmetrizers, these 8 basis elements can be rearranged into three mixed-symmetry 
doublets and two symmetric or antisymmetric singlets: 

^^^^ -S+ 

„/,(!) _ A + 

y^Ms ~ '^11 

V'S^5 = v^i-vri-Ari, 
V'£^^ = 2Sn + v+ +vri-(p+ +pri), 

^Ms = 2 An + V+ - Vn + 3 (P+ - Pn) , 
i;s = 3 (Pi - Pn) - (V+ - V^) - 2 . 



(4.35) 



Via exchange of the involved momenta, the dressing functions /j(p^, Q^, {z}) defined in 
Eq. (4.4) transform as irreducible representations of the permutation group as well. If 
those coefficients fi were totally symmetric, e.g., by being constant or by depending only 
on certain symmetric combinations of p^, and {z} as derived in [159], the nucleon 
amplitude would be a linear combination of the six and Ms basis elements in 
Eq. (4.35): 



(4.36) 



Since the coefficients fi can appear in all symmetry representations the inclusion of 
the remaining Dirac covariants ^-^d ips (and, in general, all 64 basis elements) is 
however necessary. 



4.4 Results 

The explicit numerical implementation of the Faddeev equation is described in App. A. 5. 
We solve for all 64 dressing functions fiip^ , q'^ ,0, zi, Z2) but omit the dependence on 
the angular variable zo = pt'Qt for the sake of numerical efficiency. In the context of a 
quark-diquark model, the dependence on zq is excluded a priori due to the separability 
assumption of the amplitude-'^, cf. Eq. (5.9). 

^ Note however that this assignment is not strictly vahd as the quark-diquark amphtude needs to 
be symmetrized to obtain the corresponding Faddeev amphtude, thereby changing the interpretation 
of the involved momenta. 
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Figure 4.2: Chebyshev moments of the dressing functions corresponding to the 
dominant covariants in the Faddeev amphtude ^Ma^ plotted as a function of p and q. 
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As a consequence of Eq. (4.32), the two states "^Ma "^Ms emerge as indepen- 
dent solutions of the Faddeev equation, where by virtue of the iterative method the 
symmetry of the start function determines the symmetry of the resulting amplitude. 
Both separate equations produce approximately the same nucleon mass, where the de- 
viation of ~ 2% is presumably a truncation artifact associated with the omission of 
the angle zq. For each of the two solutions, typically only a small number of covariants 
are relevant. Comparing the relative strengths of the amplitudes (cf. Fig. 4.2 for the 
mixed- antisymmetric solution) allows to identify the dominant contributions: 

^1*3' ^13' '^n ~ '^ii' (4.37) 

^53 + ^53 ^63 + ^63 

which indicates a sizeable admixture of p waves to the dominant s-wave components 
(cf. Table 4.1). 

Fig. 4.3 displays the angular dependence in the variable zi, in terms of the first 
few Chebyshev moments, of the amplitudes S^^ which contribute to ^Ma- '^^^ -^i 
dependence is much more pronounced than that in the variable Z2 where already the 
zeroth Chebyshev moment provides a satisfactory approximation. This is again kindred 
to the quark-diquark model, where q is related to the relative momentum between the 
two quarks in a diquark amplitude and the dependence on the associated angle Z2 is 
small, cf. Eq. (A.51). 
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Figure 4.4: Current-mass evolution of and nip compared to lattice results: [166- 
168] for Mtv and [134, 135] for m^. The dashed lines are related to setup (CI). The 
solid line for rup is the input of Eq. (3.15) in setup (C3); the band for M^r is the result of 
the Faddeev equation (mixed-antisymmetric solution). Both correspond to a variation 
of u. Dots demarcate the experimental values. 



The resulting nucleon masses at the physical u — d quark value in both setups (CI) 
and (C3) of Section 3.2 are: 

(CD- ^"^^ • fC3)- "^-^^ ■ ^-^^^^^ 

^X^: 0.97GeV, ^a^^ : 1.31(2) GeV, 

where the u dependence is explicitly taken into account in the 'core' setup (C3). The 
current-mass evolution of M^v is plotted in Fig. (4.4) and compared to lattice results. 
The findings are qualitatively similar to those for mp-. setup (CI), where the coupling 
strength is adjusted to the experimental value of 7,^, agrees with the lattice data while 
the input of (C3) provides a description of a quark core which consistently overestimates 
the experimental values while approaching the lattice results at larger quark masses. 
The sensitivity of M^r to the width parameter u) is more articulate than that of the 
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quark-diquark model result for the nucleon mass, cf. Section 5.3, which might also 
be a side effect of neglecting the angular variable zq. For a detailed discussion of the 
current-mass dependence of M^v we refer to Section 5.3 where the quark-diquark model 
results for M^v and Ma are compared to estimates from chiral effective field theory. 
In this context it will turn out that the quark-diquark solution for the nucleon mass 
provides a quite reasonable approximation to the result of the three-body calculation. 



Chapter 5 



Baryons: 

The quark-diquark picture 

The underlying assumption of the Faddeev truncation in Chapter 4 was the identifica- 
tion of two-quark correlations as the dominant structure which binds a baryon. Upon 
implementation of a rainbow-ladder truncation — i.e., an iterated gluon exchange, the 
same mechanism which has been used to describe qq bound states — a direct numerical 
solution of the relativistic Faddeev equation (4.2) was obtained. 

The same premise can be implemented with less numerical effort in a quark-diquark 
model. It treats such two-quark correlations as a separable pole sum in the qq scattering 
matrix and leads to a description of baryons as bound states of effective quarks and 
diquarks. In the simplest version of the model the nucleon is made of a quark and a 
scalar (0"*") diquark whereas the A comprises a quark and an axial-vector (1^) diquark. 
A strong attraction in the color-antitriplet diquark channel has also been proposed to 
explain missing exotic states in the hadron spectrum and the masses of light scalar 
mesons [169,170]. Further support for the diquark concept has recently been drawn 
from lattice calculations, cf. Table 5.2. 

The simplest realization of a quark-diquark description of baryons is the Nambu- 
Jona-Lasinio (NJL) model [171, 172] which allows for a formation of bound states 
of quarks and point-like diquarks via quark exchange [173-178]. In this context it 
was soon realized that axial-vector diquarks provide substantial attraction in the nu- 
cleon and should be taken into account as well. An extension of the model to in- 
clude diquarks of finite width led to a series of studies investigating nucleon and A 
masses [61,179,180], nucleon electromagnetic form factors [62,63,99,181-185], and the 
nucleon's pseudoscalar, scalar and axial-vector form factors [184,186]. 

The present chapter extends the quark-diquark model insofar as the dynamics of 0^ 
and 1^ diquarks are determined from their underlying quark and gluon constituents. 
Parametrizations for the diquark amplitudes are removed and replaced by solutions 
of the corresponding diquark Bethe-Salpeter equations. The identification of colored 
diquarks as poles in the qq scattering matrix becomes possible within a rainbow-ladder 
truncation: it induces timelike 0"*", 1+, ... diquark poles whose mass scales play an 
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important role in the description of Ught baryons. These poles, however, correspond 
to unphysical asymptotic states and disappear from the spectrum when going beyond 
rainbow-ladder. The quark-diquark BSE is derived from the relativistic Faddeev equa- 
tion and will be applied to calculate the masses of N and A. The current-mass depen- 
dence of these masses will be compared to lattice-QCD results, and we will discuss the 
impact of chiral corrections. 

5.1 The diquark ansatz for the qq scattering matrix 

Faddeev equation revisited. As a first step in deriving the quark-diquark BSE from 
the relativistic Faddeev equation (4.2), one introduces the Faddeev components ^'j of 
the baryon amplitude via 

3 3 
^ = ^^f)^=:^vl/.. (5.1) 

i=l i=l 

The are identical to the three graphs on the right-hand side of Fig. 4.1. Again, the 
index i refers to the non-interacting spectator quark which contributes to each diagram. 
The equation is subsequently rewritten in a way where the 2-quark scattering matrices 

~(2) ~ (2) 

T- appear instead of the 2-quark kernels . They are related to each other through 
Dyson's equation (2.18) which entails 

^ + f .(2)^ ^(2) ^ f{2)^ ^f^^^ ^(2)^ ^. ^^ 2) 

The relativistic Faddeev equation is thereby transformed into a set of coupled integral 
equations for the components 

= ) (^ - ^i) = ) (^, + ^fc) , (5.3) 

where {i, j, k} is an even permutation of {1, 2, 3}. This equation is depicted in Fig. 5.1. 

No new information has been gathered by this transformation. Instead of directly 
implementing an ansatz for the kernel in the original equation (5.1), the modified 
version necessitates an expression for the qq scattering matrix T^'^^ which must be 
determined from the kernel in the intermediate step (5.2). In particular, the Faddeev 
components are still three-body amplitudes which depend on three independent 
momenta and carry the same Dirac structure as the amplitude ^. A solution of the 
equation relies upon the solution techniques described in Section 4. 

An investigation of the structure of the T-matrix in this setup would certainly be 
an interesting issue in itself. The basic intention of Eq. (5.3) is however to replace that 
solution by an ansatz which is adequate to reduce the complexity of the equation, i.e., 
by simplifying the Faddeev approach to a two-body problem. 
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Figure 5.1: Relativistic Faddeev equation (5.3) which involves the 2-quark scattering 
matrix. 



Diquark correlations. Up to this point we have assumed that correlations between 
two quarks provide the dominant attraction not only in the meson but also in the 
baryon channel; and that these are dominated by an iterated gluon exchange between 
any pair of two quarks which constitute the baryon. In the following discussion we 
will further exploit this concept and aim at an ansatz for the T-matrix which contains 
diquark poles at timelike values of the total quark-quark momentum but still retains 
some of the characteristic features of Dyson's equation (2.18). This is realized through 
a separable sum over diquark correlations. We restrict ourselves to the O"*" and 1^ 
channels, i. e. to scalar and axial-vector diquarks for reasons explained below. The 
corresponding expression for T^^-* reads: 




where the notation (/xi^) implies that the Lorentz indices can be dropped in the scalar- 
diquark contribution. F'^'^^ are the diquark analogues of the respective pseudoscalar 
and vector meson amplitudes and the D^^^'^^ denote scalar and axial- vector diquark 
propagators. The assumed poles of the T-matrix in the ansatz (5.4) are embedded in 
the diquark propagators and define the associated diquark masses: 

^^^^ P^ + Ml' ^ ^ P^ + Ml' ^^-^^ 

where Tp' again denotes a transverse projector, now with respect to the diquark mo- 
mentum P. In the same manner as described in Section 2.3, this leads to a homogeneous 
diquark BSE, F = K^'^^T, for a diquark bound-state on the mass sheh p2 = which 
resembles the meson BSE given in Eq. (3.9). It was used for detailed studies of diquarks, 
e.g., in [187, 188]. The equation is shown diagrammatically in Fig. 5.2 and reads 

A 

^a^iP^P) = J K^yMP^<l^P){siQ+)'^^^KQ,P)S^iq^)}^^ , (5.6) 

where the momenta have been defined in the discussion of Eq. (3.9). The replacements 
S{—q-) — > S^{q-) and Ka-y^sfS ~^ ^0^,135 amount to a substitution of an antiquark- 
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p- q- 




FlGURE 5.2: The Diquark BSE (5.6) in RL truncation. 



with a quark leg. For the sake of consistency the kernel K is identified with the 
rainbow-ladder kernel (3.11). 

By working out the color traces of Eqs. (3.11) and (A. 34) one finds that the resulting 
equation for T with quantum numbers is identical to that of a color-singlet 
meson except for the diquark's coupling strength which is reduced by a factor of 2. 
This confirms that the interaction in the color anti-triplet diquark channel is strong 
and attractive. The same analysis entails that the interaction is strong and repulsive 
in the color sextet channel [74, 189]. Comparison with meson phenomenology hence 
suggests that the lightest diquarks are the scalar diquarks (the parity partners of the 
pseudoscalar mesons), followed by axial- vector, pseudoscalar and vector diquarks. This 
was also observed in Bethe-Salpeter [187] and lattice studies studies (cf. Table 5.1) and 
justifies the restriction to the scalar and axial-vector diquark channels for describing 
light baryons composed of quark and diquark. 

Diquark poles and diquark confinement. A necessary prerequisite to justify the 
particular ansatz (5.4) is that the rainbow-ladder kernel, via (5.2), indeed generates 
timelike scalar and axial-vector diquark poles in the quark-quark scattering matrix. 
While this has not been explicitly studied in the context of Dyson's equation, the 
existence of a solution to Eq. (5.6) provides a satisfactory indication. 

Asymptotic diquark states correspond to timelike poles in the diquark propagators, 
hence one might suspect a violation of diquark confinement. However, the absence 
of a Lehmann representation of a certain propagator is a sufficient but not neces- 
sary criterion for confinement of the corresponding state due the associated violation 
of refiection positivity [18,44,190,191]. Two-point correlations of colored fields may 
contain real timelike poles in momentum space without contradicting confinement, a 
statement which is also true for the quark propagator [45, 192, 193]. In addition, free- 
particle quark and diquark propagators can yield quantitatively meaningful results for 
hadronic observables (see, e.g., [184]). 

On the other hand, the introduction of beyond-RL interaction terms in the skeleton 
expansion of the qq kernel which appears in the diquark BSE removes diquark states 
from the physical mass spectrum due to large repulsive corrections [101, 102, 105, 194]. 
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In this respect, kernels that do not produce diquark bound states but induce a more 
complicated singularity structure in the qq T-matrix can still support a physical inter- 
pretation either in terms of mass scales or inverse correlation lengths of diquark inter- 
actions inside a baryon. Further motivation for the significance of the diquark concept 
has come from investigations of diquark confinement in Coulomb-gauge QCD [195]. 

Diquark masses. While diquark masses are gauge-dependent, gauge-independent 
mass differences can be determined from lattice calculations. Several such investigations 
have been performed in different approaches and with various fermion actions [196-201]. 
While they exhibit quite different quantitative results, the common qualitative feature 
is that the mass splitting between scalar and axial-vector diquark in the chiral limit 
is of the size of hundred to several hundred MeV (see Table 5.1) and decreases with 
increasing current quark mass. 

Fig. 5.3 shows scalar and axial- vector masses together with their mass and squared- 
mass difference as a result of the diquark BSE (5.6). The diquark masses exhibit large 
sensitivities to the width parameter a;, a feature which has previously been observed in 
Ref. [187]. The scalar-axial-vector mass splitting is comparatively weakly dependent 
on model details, cf. Table 5.2. The deviation in the squared-mass difference between 
setups (CI) and (C3) has its origin in the current-mass dependence of the scale Air in 
(C3) and can be removed by examining the ratio 



al) + (C - hi) . (5.7) 



Ml^ - (3^2) , 2 _ 2\ (l,^ 

^^IR ^^IR 

Since Air is constant in setup (CI), a weak current-mass dependence (dashed line in the 
upper right plot of Fig. 5.3) implies that the squared-mass difference is predominantly 
generated by the infrared contribution to the effective coupling (3.14) which owes to 
dynamical chiral symmetry breaking. 

OfFshell behavior of the T-matrix. Upon solving the scalar and axial-vector di- 
quark BSEs (5.6), the quark-quark scattering matrix is determined at the diquark mass 
poles = —M^^ and = — M^^. Since the diquarks in the nucleon are off-shell, the 
description of baryons as composites of quark and diquark requires knowledge of the 
T-matrix for general diquark momenta as well. The ansatz (5.4) dictates its off-shell 
behavior to be inherited from the separable structure in the vicinity of the poles in 
terms of off-shell diquark amplitudes and propagators. 

Given a certain ansatz for the off-shell amplitudes, one may exploit Dyson's equation 
(2.18) to obtain an expression for the diquark propagators. This procedure is detailed 
in Appendix A. 4 and yields propagators of the form 

D{P^) = R{P-) + (5.8) 



where the finite parts R{P^) emulate off-shell contributions which are suppressed at 
the mass poles but determine the ultraviolet behavior of the T-matrix. 
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Figure 5.3: Left panel: Scalar and axial- vector diquark masses from the diquark BSE 
vs. squared pion mass [113]. The bands denote a variation of a; in setup (C3); the 
central lines correspond to u). For comparison we show the (C3) input for rup {solid 
line), together with its experimental value {star), lattice data and a chiral extrapolation 
{dashed line). Respective references are given in Fig. 3.5. Right panels: Scalar-axial- 
vector mass splitting {lower panel) and squared-mass splitting {upper panel). Dashed 
lines are the results of (CI); bands those of (C3). 





Air 


7? 






(CI) 


0.72 


1.8 


0.81 


1.00 


(C3) 


0.98 


1.8(2) 


1.08(10) 


1.35(8) 



Table 5.1: Comparison of diquark masses Msc, -^av as obtained from the diquark 
BSE in setups (CI) and (C3), characterized by the parameters Air and r] defined in 
(3.19). A variation of r/ = 1.8 it 0.2 in (C3) is equivalent to tj ~ u; it 0.06. The results 
correspond to a current mass m = 6.1 MeV which is related to the physical pion mass 
= 138 MeV. The units of Air, M^c and Mav are GeV; r] is dimensionless. 



(CI) (C3) 


Ref. [196] 


Ref. [197] 


Ref. [198] 


Ref. [199] 


0.20 0.27(3) 


0.10(5) 


0.14(1) 


0.29(4) 


0.36(7) 



Table 5.2: Scalar-axialvector diquark mass splitting in the chiral limit, cf. Fig. 5.3. 
The BSE values are compared to several lattice-QCD results. The units are GeV. 
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5.2 Quark-diquark BSE 

Separable ansatz for the amplitudes. The relativistic Faddeev equations in the 
form (5.3) represent coupled integral equations for either of the Faddeev components 
All elements of the equation, i.e., the dressed-quark propagator and the diquark 
ingredients of the quark-quark scattering matrix, have been specified in the previous 
sections. The separability of the diquark ansatz (5.4) for the T-matrix is one of two 
requirements to reduce the original three-body equation for the baryon to a two-body 
Bethe-Salpeter equation for a quark-diquark bound state. The second prerequisite is 
the separability of the components in terms of the relative momenta between quark 
and diquark (p) and within the diquark, where the latter is now denoted by pr, cf. 
Fig. 5.4. This is realized via the following ansatz for the ^f. 

^af3^5{p,Pr,P) = Y,T%^{Pr,Pd) D''\pd) ^""^^iP^ P), (5.9) 
a,b 

which involves a combination of scalar and axial-vector diquark amplitudes and propa- 
gators, pd and P are the total diquark and nucleon momenta, respectively. The super- 
scripts a, 6, c collect the diquarks' Lorentz indices: a = 5 denotes scalar and a = 1 ... 4 
axial- vector quantities; the diquark propagator is either scalar (a = 5 = 5) or axial- 
vector (a,b = 1 ... 4). The Dirac, color and flavor decomposition of the quark-diquark 
amplitudes $|^^(p, P) thereby defined is described in App. A.6. The final baryon am- 
plitudes, as stated for the nucleon in Eq. (A. 74), are constructed by symmetrizing the 
Faddeev components in Eq. (5.9). 

The spin- and isospin-i/2 nucleon is a sum of scalar and axial-vector diquark corre- 
lations, whereas the SU{2) flavor algebra excludes scalar diquarks from participating 
in the spin- and isospin-3/2 A baryon: only an isospin-1 axial-vector diquark can be 
combined with an isospin-i/2 quark to obtain / = 3/2. The quark-diquark amplitude 

thus appears in the following manifestations: 

AT: ^-{p,p) a>5(p,P), $'^(p,P), 

A : $"(p,P) — > <^^''{p,P). ^ ■ ' 

The final N and A quark-diquark spinors are obtained upon contraction with the Dirac 
and Rarita-Schwinger spinors us{P), u'^{P) which describe free spin-i/2 or -^2 particles 
with momentum P: 



N: <^is{p,P)m{P),KsiP^P)MP), . . 

A: $^^(p,P)<(P). ^ • ' 



Two-body equation. Inserting the ansatz (5.4) for the T-matrix together with (5.9) 
into the relativistic Faddeev equation (5.3) yields a quark-diquark Bethe-Salpeter equa- 
tion on the baryon's mass shell [180, 184]: 
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Figure 5.4: The quark-diquark BSE (5.12) in pictorial form. 




where the kernel of the equation is given by 



K%{v,k,P) 



(5.13) 



It couples scalar and axial- vector diquarks (i.e., a = 1 ... 5) in the case of the nucleon 
whereas only the axial-vector index a = 1 ... 4 is required for the A. The kernel (5.13) 
describes an iterated exchange of roles between the spectator quark and the quarks 
which constitute the diquark; this quark exchange generates the attractive interaction 
that binds quarks and diquarks to a baryon. 

The momenta in Eqs. (5.12-5.13) read (cf. Fig. 5.4): 



kn 



P + 7] P , 

k + r]P, 



Pd 
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-p+{l-rj)P, 
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(5.14) 



(l=Pd 

kr = {1 — a) Pq — a q . 
Pr = {1 - a) Hq - a q . 

p and Pr are the relative momenta between quark and diquark and within the diquark; 
Pd, P are total diquark and nucleon momenta. The respective momenta k appear 
inside the loop integral. Again, the momentum partitioning parameters cr, rj ^ [0, 1] for 
diquark and quark-diquark amplitudes are arbitrary. In contrast to the analogous case 
in the three-body equation, r] = 1/3 is no longer the optimal choice since the pole limits 
now result from the combined singularity structures of quark and diquark propagators 
(and, theoretically, also of the diquark amplitudes). We set cr = 1/2 but keep as a 
variable since it can be used to ease these constraints (cf. App. B.3). 

Upon working out the color and flavor factors of the quark-diquark amplitudes, 
given in Eqs. (A. 87) and (A. 92), and of the diquark amplitudes (A. 34), the BSE kernel 
picks up a color-flavor factor 

1/-1 V3\ (0 0- 

2 I 73 1 i ' I -1 



N : 



(5.15) 
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where the first row (column) represents the scalar part of the kernel and the second 
row (column) the axial-vector part. 

5.3 Nucleon and A masses 

The results of the quark-diquark BSE (5.12) for nucleon and A masses in the 'quark 
core' setup (C3) are shown in Fig. 5.5. The left panel depicts the calculated values 
for nip, Mat and Ma, each together with a selection of lattice results and their chiral 
extrapolations (if available). The corresponding abscissa values are obtained from 
the pseudoscalar meson BSE. The solid curve for nip is the input defined in Eqs. (3.15- 
3.17) and completely fixes the parameters in the interaction. The bands represent 
the sensitivity of the results for Mjy and Ma on the width parameter uj. At larger 
quark masses the deviation from the lattice data diminishes, in accordance with the 
assumption of Eq. (3.15), namely that beyond- RL corrections to hadronic observables 
become small in the heavy-quark limit. 

The resulting values at the physical pion mass in both setups (CI) and (C3) are 
displayed in Table 5.3. Again, Eq. (3.22) can be used to relate both models to each 
other. The intention of the former was to reproduce vr and p properties; in addition, 
the model also yields and A masses that are close to the experimental values. This 
is quite remarkable since, upon fixing the coupling strength in the effective quark-gluon 
coupling to meson phenomenology, no further parameters have been used as an input 
of the calculation. 

A comparison with the results of Section 4.4, obtained through a solution of the 
relativistic Faddeev equation, shows that the corresponding nucleon mass is larger 
than the result of the quark-diquark approach. Fig. 4.4 indicates an approximately 
uniform shift by 50 — 100 MeV throughout the examined current-quark mass range. 
This illustrates that the quark-diquark approach is somewhat too attractive in the 
nucleon channel but still a quite reasonable approximation. 

Pionic corrections. It is instructive to compare our results to baryonic core masses 
estimated from chiral effective field theory. In this framework the mass of a baryon B 
is obtained from the expression (e.g., [58]) 



where the "baryon core" includes the bare (and a priori unknown) parameters Oj ap- 
pearing in the effective Lagrangian: 



The inherent assumption of a momentum cutoff regularization, expressed through the 
scale A, allows for a physical interpretation in terms of a scale separation between quark 
core and pion cloud. The cutoff A regularizes the short-distance divergences associated 
with pointlike baryons and pions in the effective field theory and retains the long- 
distance or low-energy part of the self-energy integrals. A non-pointlike nucleon-pion 



MBiml) = Mr%ml,A) + ^^(m^, A) 



(5.16) 



Mr%ml,A) = ag)(A) + ag)(A) ml + ag)(A) + . . . 



(5.17) 
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Figure 5.5: Evolution of N and A masses {left panel) and the mass sphttings Ma — 
Mjv and Mav — Msc {right panel) vs. pion mass squared [113,202], obtained in the quark- 
diquark model using setup (C3). The bands denote the sensitivity to a variation of 
w. We compare to a selection of dynamical lattice data and their chiral extrapolations 
{dashed lines) for M^r [166-168] and A/a [203]. The depicted data for nip are identical 
to those of Figs. 3.5 and 5.3. Stars denote the experimental values. 





Air 
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Mn 


Ma 


Exp. 






0.94 


1.23 


(CI) 
(C3) 


0.72 
0.98 


1.8 
1.8(2) 


0.94 
1.26(2) 


1.28 
1.73(5) 



Table 5.3: Comparison of Mjv and Ma for both model inputs (CI) and (C3) at 
the physical point. The parameter rj is defined in (3.19). The first row quotes the 
experimental values. The units of Air, Msc and Mav are GeV; r] is dimensionless. 
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Figure 5.6: Expansion of the nucleon propagator in chiral perturbation theory. 
Dashed, solid, and thick solid lines correspond to pseudoscalar-meson, nucleon and 
A degrees of freedom. 
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interaction requires a cutoff of the order of the baryon size (A > 0.2 GeV ~ r < 1 fm). 
Of course one may alternatively employ a dimensional regularization technique [204] . 

in Eq. (5.16) denotes the sum of meson-loop self-energies, e.g. to 1-loop order 
the sum of intermediate A^vr and Avr states (cf. Fig. 5.6). In heavy baryon chiral 
perturbation theory [205,206], the respective contributions are given by [58,207]: 



EA(m2,A) = -AA3 



dx 



(x) + 



32 



25 LO + 6 
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Una 



(x) 



(5.18) 



dx 



25 uj- 5 



u%^{x) 



where x, uj{x) and 6 are pion momentum, intermediate pion energy and the physical 
A^-A mass splitting, normalized by the cutoff A, and the prefactor A includes the 
experimental values for the pion decay constant = 131 MeV and the axial coupling 
9A = 1.26: 
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uj{x) 
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0.29 GeV 
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(5.19) 



The cutoff-dependent NNtt, A Avr and A^Att vertex dressings are denoted by ubb'{x). 
A dipole form factor mimics the physical shape of the meson-baryon vertex. Choosing 
all dressings identical via u{x) = (1 -|- x^)~^ with a regulator A = 0.8 GeV yields for 
= 140 MeV: 

E^r = -0.30GeV, Sa = -0.27 GeV, (5.20) 

and hence a similar reduction of both nucleon and A masses. In combination with the 
experimental numbers for and Ma, (5.20) provides the simple estimates 

M^™^ ~ 1.24 GeV, M^'"" ~ 1.5 GeV. (5.21) 

These values are roughly consistent with those obtained in, e.g., the cloudy-bag model 
[208], NJL model [209], and nucleon-pion Dyson-Schwinger studies [61,210] and make 
clear that nucleon-pion loops are attractive and the binding energy reduces the baryon's 
mass. 

The separation into a "core" and a "meson cloud" contribution is cutoff-dependent. 
Due to the approximate A^-dependence of the self-energy integrals, different values of 
the regulator can vary the above results considerably. Naturally, the chiral extrapola- 
tions of the physical N and A masses depicted in Fig. 5.5 must be independent of the 
cutoff A. Via expansion in m^, the self-energy contributions are split into cutoff- 



independent non-analytic terms 



~ m^, ~ m^lnm,r and cutoff-dependent terms that 

(i) 
B 



are even in [168,211,212]. The latter are combined with the bare parameters a 
to renormalized coefficients such that the chiral expansion of Mjv and Ma reads [213] 

MB{ml) = {Mb{^) + cf ml + c 
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(5.22) 
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It is expressed in terms of renormalized low-energy constants Mb{0), and which 
may be determined by a fit to lattice data [58, 168]. Through this renormalization 
procedure the regulator dependence is removed [204,214]. 

Discussion. The previous considerations illustrate that both the quark-diquark model 
and the three-body result for M^r in the 'core' setup (C3) are roughly consistent with a 
pseudoscalar-meson dressing providing the dominant correction to the nucleon's quark 
core. In particular, it complies with the assumptions which motivated the introduction 
of (3.15) in the context of pion and p-meson observables. 

Nevertheless one has to keep in mind that the identification of the baryonic quark 
core (5.21) with the quark-diquark or three-quark 'core' is more complicated than in 
the meson case. Eq. (3.15) assumes that corrections to rrip are partly induced by a 
pseudoscalar-meson cloud and to a similar extent owe to non-resonant corrections to 
RL truncation. In the baryon one has additional lines of improvement: an inclusion of 
irreducible 3-body interactions, which were neglected in the derivation of the Faddeev 
equation (4.2), can still describe a quark core in the sense of Eq. (3.15). Morever it 
is conceivable that contributions beyond rainbow-ladder induce different ramifications 
for rup and M^r. 

An interpretation of the quark-diquark result for Ma is more difficult. In view 
of Eq. (5.21), the solution depicted in Fig. 5.5 indicates larger corrections beyond the 
quark-diquark approach than encountered for the nucleon. The Ma solution addition- 
ally exhibits a sizeable uj dependence, a feature which is less pronounced in vr, p and 
nucleon observables. The scalar and axial- vector diquark masses exhibit particularly 
large sensitivities to uj (see Fig. 5.3) which apparently cancel upon constituting the 
nucleon mass. In the Ma case, the same consideration could suggest that taking into 
account only an axial-vector diquark may not be sufficient for describing the A, and 
that a possible further isospin-1 (tensor) diquark component with a mass large enough 
to be irrelevant for the nucleon could diminish the A 'core' mass. In this respect it 
is highly desirable to extend the three-body study of the nucleon to the A. Such a 
framework automatically implements the effect of inserting further diquark channels 
which could affect N and A properties differently. 

The right panel of Fig. 5.5 displays the A^-A mass splitting. According to (5.20), 
the pseudoscalar meson contribution to the experimental value Ma — M^y = 0.29 GeV 
is small and positive: Ea — Sjv = 0.03 GeV, cf. Ref. [58]. This is not the case in 
our calculation, where at the u/d mass (Ma — M^Y"^'^ = 0.48(4) GeV and therefore 
predicts a negative correction to the full splitting. 

We also compare Ma — Mjy with the diquark mass splitting Mav — Mgc- Both 
decrease with increasing current-quark mass; nevertheless there is no direct relationship 
between the two quantities, since the axial- vector diquark contribution to the mass of 
the nucleon does not vanish [202]. 



Chapter 6 

Nucleon: 

Electromagnetic form factors 

The nucleoli's charge and magnetization structure is encoded in its electromagnetic 
form factors which, for space-like values of the photon momentum, are experimentally 
determined via elastic nucleon-electron scattering. Nucleon electromagnetic form fac- 
tors have been studied in a variety of approaches; an overview on the experimental and 
theoretical progress can be found in the recent review articles of Refs. [215,216]. 

The results presented in this chapter are rooted in a long tradition of nucleon 
form factor studies within the quark-diquark model [99,181-185]. These calculations 
share some common caveats. First, pionic contributions play an important role in the 
low-energy and small-quark mass behavior of the nucleon's electromagnetic structure. 
Such effects are not included in a quark-diquark 'core' and must be added on top of 
it [61-63,210]. Second, access to the large-Q^ region and thereby to the truly pertur- 
bative domain is so far only feasible upon implementing pole-free model propagators 
which, in turn, exhibit essential singularities at timelike infinity. The problem is not 
of fundamental concern; it merely awaits a thorough numerical treatment. Third, the 
quark-mass dependence of magnetic moments and charge radii, while emerging natu- 
rally in lattice calculations, is practically inaccessible in a quark-diquark model due to 
the unknown mass dependence of the modeled ingredients. 

The rainbow-ladder based quark-diquark approach, introduced in Chapter 5, re- 
moves the latter obstacle. Upon resolving the diquarks' substructure, the form factors 
are immediately related to the parameters in the effective quark-gluon coupling a(A;^), 
in particular: its quark-mass dependent coupling strength. In the following we will 
restrict ourselves to the quark 'core model' which represents a quark-diquark core that 
needs to be dressed by meson-cloud effects. A comparison of the core's static properties 
with lattice results is appropriate at larger quark masses; form factors depending on the 
photon momentum may be compared to experiment at > 2 GeV^ where pion-cloud 
effects are diminished. 
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6.1 Electromagnetic current 

Nucleon-photon vertex. The vertex which describes the coupling of the nucleon as 
a spin-1/2 fermion to a photon is constructed from 12 Dirac covariants: 

{7^p^Q^}x{l,f,^, [f,^]}, (6.1) 

where the Breit momentum P and photon momentum Q are combinations of incoming 
and outgoing nucleon momenta: 

P=^{P, + Pf), Q = Pf-Pi. (6.2) 

A nucleon on the mass shell satisfies Pf = P'j = — and therefore = — — 
and P • Q = 0; hence the only independent variable is the photon momentum-squared 
Q'^. The nucleon's electromagnetic current J^{Q'^) is obtained by sandwiching the on- 
shell vertex by the nucleon spinors C/^(P/), V^iPi) of Eq. (4.25) which are solutions of 
the Dirac equation: 

A+(P,) C/"(P,) = C/°(P,), Uf'iPf) A+{Pf) = U^iPf) . (6.3) 

Equivalently, one may construct a matrix- valued current by taking the spin sums, i.e., 
by contracting the vertex with the positive-energy projectors 



^ ^ a=l 



\f) ■ (6.4) 



They reduce the 12 basis elements to three: {7'^, P'^, Q'^}. Current conservation 
Qf'Jf'iQ'^) = eliminates the structure Qf" via A+(P/) ^ A+(P) = and Q • P = 0. 
The most general electromagnetic current of the nucleon is therefore given by 



Ji^iQ') = A+(P;) ((Pi + P2) ij'^ - F2 ^) A+(PO , 



(6.5) 



thereby defining the Dirac and Pauli form factors Pi((5^) and F2{Q'^). Using the 
Gordon identity 

A.(P/)(n'^ + ^-^)A,(^0=0 (6.6) 
with a^'^ = — J [7^,7*^] transforms Eq. (6.5) into: 

J^iQ') = A+{Pf) (^Pi i^-^i cr^-'Q"^ A+(PO . (6.7) 

Pi and P2 are dimensionless; for = they reduce to the proton and neutron charges 
= {1)0} and anomalous magnetic moments expressed in nuclear magnetons 
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eh/{2M). Charge conservation -Fi(O) = A is automatically satisfied if the nucleon 
amplitudes are canonically normalized via Eq. (2.24) [183]. 

Electromagnetic current of a composite system. To relate the electromagnetic 
current of Eq. (6.7) to the underlying description of the nucleon as a composite object, 
the baryon must be resolved into its constituents to each of which the current can 
couple. The construction of an electromagnetic current operator in the framework of 
the Bethe-Salpeter equation was first treated by Mandelstam [217]. In the same way 
as an n-particle T-matrix reduces to the form (2.21) at the bound-state pole with mass 
M, thereby defining the bound-state amplitude, the electromagnetic current matrix 
is the residue of the (2n -|- 1)— point function T'^ which describes the photon's coupling 
to the T-matrix at the bound-state mass: 

^ ' (P2 + + Q2/4)2 • (^-^^ 

The Breit momentum P is defined in (6.2) and yields at the hadron pole: 

[P] + M'^){Pf + M^) = (p2 + + q2/4) = 0. (6.9) 

A systematic procedure for the construction of a hadron-photon vertex based on 
electromagnetic gauge invariance is the "gauging of equations" prescription [218-220] 
which represents a generalization of the normalization condition in Eq. (2.24) to finite 
photon momenta. In this context, "gauging", formally denoted by T ^ T'^, is a 
derivative: it is linear and satisfies Leibniz' rule. The gauged n-quark scattering matrix 
satisfies 

^ _j, (rp-iyj, ^ _j, ^^-1 _Q^Yj,^ri.^Q,i^ K-^K^R-^) T, (6.10) 

where, according to Section 2.3, K is the interaction kernel. Go is the product of 
n propagators, and Dyson's equation has been implemented. The gauged T-matrix 
depends on incoming and outgoing total momenta Pi, Pf (or, equivalently, on P and 
Q) and two further relative momenta which are not relevant for the following discussion. 

Eq. (6.10) entails in combination with the pole behavior (2.21) and the bound-state 
equation ^ = KGq^: 

^!=HM^ A^pISp (GS + 0„ G„, J%. (0.11) 

and comparison with (6.8) yields 

= (2^-i)M ^. ^ + j^i. ^. _ (e_i2) 

The impulse approximation Gq involves, e.g., the quark-photon vertices (—5"^^)^ where- 
as K^^ represents the photon's coupling to the kernel K. For instance, in the case of 
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Figure 6.1: The two types of diagrams which contribute to the nucleon's three-body 
current in rainbow-ladder truncation, Eq. (6.16) 



a two-body system with a rainbow-ladder kernel, K does not depend on the total 
momentum, hence = 0: the electromagnetic current of a meson consistent with 
a rainbow- ladder truncation is an impulse- approximation current. This is no longer 
true for a baryon, either described in terms of three valence quarks or a quark-diquark 
system. 

For vanishing photon momentum, 'gauging' is the derivative with respect to P^: 



_ dT 



7 

2P^^ =: 2P^T' 



(6.13) 



Comparison with the normalization condition N (T^^^'^' = 1 of Eq. (2.24) yields 
NJ^{Gi) = 2P^ (modulo the bound-state wave functions, e.g., the nucleon spinors). 
Examples for this relation are given in (A. 46) for the pion's electromagnetic current 
(M = 1) and (6.7) for that of the nucleon [M = 2M). At the level of the bound-state 
constituents which contribute to the total current this relation is ensured by differential 
Ward identities. 

Three-body approach. The kernel which appears in the bound-state equation of a 
baryon described by three valence quarks (cf. Section 4) is given by 



(2) 



i=l 



The corresponding gauged kernel reads 

3 



i^(3) 



(3) 



i=l 



(2) 



i=l 



(6.14) 



(6.15) 



Neglecting the three-body irreducible contribution and using a rainbow-ladder two- 
body kernel therefore leads to the electromagnetic current 



G(; + Go^i^f)0(5ri)^Go 



i=l 



(6.16) 



6.2 Electromagnetic form factors 



65 



which is depicted in Fig. 6.1. We will not further investigate this equation in the present 
work and focus instead on the respective current operator in the quark-diquark model. 

Diagrams in the quark-diquark model. The general discussion of the last section 
applies to the quark-diquark model as well. In this context, the incoming and outgoing 
baryon states are described by quark-diquark amplitudes ^>j, $ / introduced in Eq. (5.9). 
Upon interaction with the external current, the baryon is resolved into its constituents: 
quark and diquark and the interaction between them to each of which the current can 
couple [183]. 

In the context of Eq. (6.12), T is the quark-diquark scattering matrix, Gq = S D the 
product of a dressed quark and diquark propagator, and K = T ST the quark-diquark 
kernel describing the quark exchange. The quark-photon and di quark-photon vertices 
are defined as the gauged inverse propagators: Tq := — (^S~^)^ and T'^^ := — (^D~^)^. 
The gauged diquark amplitudes T'^ =: M^, also referred to as seagull vertices, describe 
the photon coupling to the diquark amplitudes. The ingredients of the current matrix 
are then written as 



These diagrams are worked out in detail in App. A.9 and illustrated in Fig. 6.2. 

By virtue of Eq. (6.17), the electromagnetic current of a baryon in a quark-diquark 
framework is completely specified by identifying the quark-photon vertex, the scalar 
and axial-vector diquark-photon vertices and an ansatz for the seagull terms. These 
quantities are constrained by Ward-Takahashi identities and thereby related to quark 
and diquark propagators and diquark amplitudes which have already been determined 
previously. The corresponding vertices are collected in Apps. A.l, A. 7 and A. 8. 

6.2 Electromagnetic form factors 

Sachs form factors. Defining r = Q^/(4M^), the correspondence between the electric 
and magnetic Sachs form factors to the Dirac and Pauli form factors which appear in 
(6.7) is given by 



= {S Df = {3^^8)0 + S{D r;^q D) , 

= [T STY = Sf + T [ST'^S)f + T SM^ . 



(6.17) 



Ge = Fi- tF2 
Gm = Fi+F2 



Ge + T Gm 



1 + r 



F2 



Gm — Ge 



1 + r 



(6.18) 



They can be extracted from the current via 



p2 



Ge, -T^'TriJf'i^''} =tGm. 



(6.19) 



M 
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Figure 6.2: The five diagrams tliat constitute the nucleoli's electromagnetic current 
in the quark-diquark model, corresponding to Eq. (6.17). 



Unless further specified, the symbol F{Q'^) denotes a generic form factor in the following 
discussion, either Dirac/Pauli or Sachs. The electromagnetic radius rp corresponding 
to F{Q'^) is defined as the slope at zero momentum transfer via the Taylor expansion 



m 



A 



1 
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rp 



rp 



-6 



dF 
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dF 



F{0) dQ^ 



(6.20) 



(6.21) 



where F(Q'^) in the first row denotes an 'electric' form factor (Fi or Gp) with the 
electric charge F(0) = A, in the second row a 'magnetic' form factor (F2, Gm) with 
the (anomalous) magnetic moment -F(O) = k or /i such that Gm(0) = fi = X + n. With 
Eq. (6.18) the correspondence between Sachs and Dirac/Pauli radii is 



r-2 



3k 
2M2 



r2 



r J + K r2 
1 + K 



(6.22) 



where the second term in r|; involving the anomalous magnetic moment is the so-called 
Foldy term. Experimental values for electromagnetic radii and magnetic moments are 
collected in Table 6.1. 
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[fm] 


/ \ 2 


rP = 0.886(15) 
(r|)2 = -0.115(4) 
= 0.855(35) 
= 0.873(11) 


= 1.793 

Kn = -1.913 

Ks = -0.120 
Kt, = 3.706 


= 0.82(2) 
r'1 = 0.11(2) 
= 0.88(6) 
= 0.88(1) 


(rf)2 = 0.68(3) 
{r\f = 0.65(3) 
(r|)2 = 0.78 ±2 
{r^^f = 0.77(6) 



Table 6.1: Experimental numbers for radii and magnetic moments. The radii in 
the left column are quoted from Table 1 of Ref. [225]; for references, see therein. 
Magnetic moments are the PDG values [226]. The remaining quantities are inferred 
from Eqs. (6.18) and (6.39), ignoring possible correlations between the statistical errors 
of the data. 



Form factors are Lorentz invariant, but their interpretation depends on the reference 
frame. Nonrelativistically, the Sachs form factors Gj(Q^) are the three-dimensional 
Fourier transforms of the nucleon's spatial charge and magnetization distributions pch, 
Pm in the Breit frame where the energy of the transferred photon is zero [221]. In this 
case the interpretation of the inferred Sachs radii as charge and magnetization RMS 
radii is valid: 

r| = 47r j drpch{r)r'^, rfj = An j drp^{r)r^, X = 4tt j dr p^hir) r"^ , (6.23) 

where A is the nucleon's electric charge. Relativistic boost corrections obscure this 
interpretation, and the extraction of these densities from experimental form factor 
data becomes model-dependent [215,222]. Such analyses implicate a positive central 
charge density of both proton and neutron [223] , in agreement with the notion of the 
neutron being sometimes a proton surrounded by a negatively charged pion cloud. 
On the other hand, the Dirac form factor Fii^Q"^) non-ambiguously corresponds to 
the nucleon's transverse charge density in the infinite momentum frame [224] , and the 
experimentally observed negativity of -^"(Q^) induces a negative central charge density 
of a fast-moving neutron. It has been argued in Ref. [113] that a negative Dirac form 
factor for the neutron can be explained in the context of a quark-diquark model which 
includes axial- vector diquark degrees of freedom (see also Section 6.3). 

Form factor measurement. Electromagnetic form factors are experimentally studied 
via elastic nucleon-electron scattering. Due to the smallness of the electromagnetic fine- 
structure constant a = 1/137, the Born approximation which describes the scattering 
process in terms of a single exchanged photon has been commonly employed. The 
scattering off a spin-1/2 target with extended structure and an anomalous magnetic 
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moment is then given by the Rosenbluth cross section [227]: 

(6.24) 



da _ f da\ eG\ + T 



dfl Vc^^/Mott + 

where the prefactor represents the Mott cross section of a spin-0 point particle, e G [0,1] 
is the virtual photon polarization, and 9 the scattering angle: 

da\ _ a^E' I 

dnj^^^^ 1 + 2 (1 + r) tan2 0/2 " ^" ^ 

E and E' are the initial and final electron energies. The Sachs form factors are extracted 
from the e dependence of the quantity e G\ + tG\j at fixed which entails a reduced 
sensitivity to at large and G m at small photon momenta. Measurements of the 
neutron form factors, carried out via electron-deuteron scattering due to the lack of a 
free neutron target in nature, suffer from systematic uncertainties. 

Recent polarization-transfer experiments allow for an increased accuracy through a 
direct extraction of the form factor ratio Ge/Gm via [228] 

Ge PtE + E' ^, , , 

where Pi and Pt are the longitudinal and transverse polarization components of the 
recoil proton, transferred from the longitudinally polarized electron. 

Surprisingly, the corresponding measurements [229-231] indicated a linear decrease 
in that ratio with increasing . This is incompatible with the result obtained by the 
Rosenbluth separation technique where the ratio is roughly ~ 1 (both are compared in 
Fig. 6.10). The discrepancy is currently believed to originate from two-photon exchange 
corrections which have a minimal impact upon the polarization results but significantly 
affect the Rosenbluth cross section [232,233]. 

Phenomenological aspects. Dimensional counting rules of perturbative QCD pre- 
dict the following behavior of the Dirac and Pauli form factors at large photon momen- 
tum transfer [234]: 

Fi~l/Q\ F2~l/Q^ Q^F2/Fi^ const., (6.27) 

where logarithmic corrections [235] have been neglected. Correspondingly, the Sachs 
form factors scale as Ge,m ~ which implies that the ratio Ge/Gm becomes 

constant. In this respect, dipole-like parametrizations for the Sachs form factors were 
found to provide a reasonable description of the experimental data: 



Gi{Q'^) = ^'1, 2 ' A = 0.84GeV, (6.28) 



(1 + Q2/A2)^ 

with gi{Q^) = const, except for the neutron electric form factor where the Galster 
parametrization [236] 5^(Q^) = —finT/il + 5.6 r) has been frequently employed to fit 
the data. 
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The deviation from the dipole form, expressed through the dependence of the re- 
mainders gi{Q'^) on the photon momentum, is nevertheless sizeable for > 2 GeV^. 
It is especially pronounced in the electric form factor of the proton: the polarization- 
transfer data for the form factor ratio G^/G^j (so far only available for be- 
low 6 GeV^) show a linear fall-off in and even point towards a zero crossing at 
~ 8 GeV^, implying the presence of a further small scale ~ 0.07 fm. The dis- 
crepancy between the perturbative prediction and the experimental data has been 
attributed to the presence of non-zero quark orbital angular-momentum content in the 
proton [185,237,238]. Possible evidence for an onset of the perturbative scaling can be 
observed in the data above ^ W GeV^ [215]. 

Pion-cloud effects are expected to play an important role in the form factor struc- 
ture below ~ 2 GeV^. Such contributions are suppressed at large where the 
photon probes the nucleon's quark core. Attempts have been made to attribute low-Q^ 
systematics in the form factor data to pionic effects, e.g. through a phenomenological 
double-dipole fit [239], in a dispersion-relation approach [240], or via implementation 
in a chiral quark model [3]. In analogy to the discussion of the pion's charge radius in 
Chapter 3.3, pion-loop contributions are believed to provide sizeable additions to the 
nucleon's quark core radii and magnetic moments, where the overall strength of these 
effects can be estimated from chiral perturbation theory. 

Flavor contributions to form factors. Proton and neutron form factors can be 
combined to study the flavor dependence of the nucleon's charge and magnetization 
structure. Each diagram of Fig. 6.2 can be split into terms where the photon either 
couples to a u or a d quark inside the nucleon. The form factors (again generically 
denoted by F) are therefore linear combinations of these u- and d-quark contributions: 



FP = 2q^ FS + q, Fl = \ [AF^ - Fj) = \ (aF^ - F^) , 
F- = q^F: + 2qaF^ = \ (F„" - F^) = \{f' - 



(6.29) 



where because of charge symmetry the u{d) contribution in the proton equals the d{u) 
contribution in the neutron: 

F^:=FP = FS, F'':=FI = F:. (6.30) 

The definition (6.29) implies F"(0) = F'^(O) = 1 for the 'charge' form factors Fi and 
Ge which entails 

{rlf = {<f + \ {rlf , {rff = {r{f + 2 {r^,f (6.31) 

for the charge radii. The flavor contribution to the magnetic form factors i<2, Gm is 
usually separated via /i" = 2 FS{0), /J-d = QdF^{0), and charge symmetry implies 

/Xp = //" + /, ^„ = -2/-l/x". (6.32) 
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The contribution from strange quarks, implicit in the seaquark content, is missing 
in Eq. (6.29) since we are working in isospin-S'C/(2). Nevertheless it has been argued 
that the present world data are consistent with the proton's strange form factors being 
zero [241]. Under this assumption one may extract the u and d contributions from the 
experimental quantities. Using the values of Table 6.1 yields: 

r]^ = 0.82(2) fm, rf = 0.83(2) fm , = 2.44 , / = 0.34. (6.33) 

6.3 Results and discussion 

This section provides results for the nucleon's electromagnetic form factors in the quark- 
diquark calculation where the current is constructed from the diagrams in Fig. 6.2. 
Apart from the ingredients which have been determined in previous chapters (nucleon 
amplitude. Section 5.2; quark propagator. Section 3.1; diquark amplitudes and prop- 
agator. Section 5.1 and the respective appendices), the calculation involves a quark- 
photon vertex (App. A. 2) and the effective diquark vertices (Apps. A. 7 and A. 8). Each 
vertex satisfies a Ward-Takahashi identity which, in total, ensures conservation of the 
nucleon's charge: Ff (0) = 1, F^{0) = 0. 

Form factor contributions. In Fig. 6.3 the proton's and neutron's electromagnetic 
Sachs form factors at the physical point are compared to experimental data. The main 
difi'erence between the results presented herein and those of Ref. [129] originates from 
the inclusion of the p-meson pole in the quark-photon vertex, cf. Eq. (A. 28). In accor- 
dance with its efi'ect on the pion charge form factor [60,83], it reduces the form factors 
and contributes 50% to their charge radii. Specifically, it cancels the previously positive 
result for to zero (within the considered domain of the uj parameter). Note that the 
form factors in terms of the underlying quark-photon and diquark-photon vertices are 
well constrained up to ~ 2 GeV^; the arbitrariness introduced in connection with 
the transverse seagull term of Eq. (A. 138) becomes important only at larger photon 
momenta. 

Fig. 6.5 depicts the contributions to the form factors G^" and Fg'" from the quark- 
photon coupling, the diquark-photon coupling, and the exchange and seagull diagrams. 
The strongest contribution at = is the direct coupling of the photon to the quark 
line. This feature may change with the evolution: for instance, the quark-photon 
contribution to exhibits a zero crossing at Q'^/M'^ PS 2 which was observed in the 
model of Ref. [63] as well. The diquark contribution to the magnetic form factors (in 
particular, the scalar-axial-vector transition) provides only a small fraction of the total 
result. The unit charge of the proton results from the canonical normalization of the 
quark-diquark amplitude. Current conservation ensures a vanishing neutron charge in 
terms of a cancelation of the components in at = 0. 

The relative strengths of the scalar-scalar, axialvector-axialvector and scalar-axial- 
vector components, according to the type of the incoming and outgoing nucleon ampli- 
tudes, can be read off from Fig. 6.6. The dominant contributions are provided by the 
scalar-diquark contributions. 
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Figure 6.3: Results for the Sachs form factors at the physical point. The selection of 
experimental data is based on Ref. [239] (data compiled by P. Grabmayr). The bands 
correspond to a variation of uj in setup (C3). 





Figure 6.4: Contributions to from the Ball-Chiu (BC) and vector-meson part (p) 
in the quark-photon vertex. The overall effect on the magnetic form factors is similar 
to that in G^. 
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Figure 6.5: Individual contributions to Ge and F2 according to Fig. 6.2. Depicted 
are the form factor sum (SUM), photon-quark (Q) and photon-diquark couphng (DQ), 
and the combined exchange and seaguh terms (EX + SG). The strength of these con- 
tributions in the remaining form factors can be extracted from the plot via Eqs. (6.18), 
(6.39) and (6.29): Fi = Ge + tF2 and Gm = F1 + F2; for the isoscalar/isovector combi- 
nations: s = p + n, V = p — n; and the up- and down-quark contributions: u = p + n/2, 
d = p + 2n. The iv variation is identical to Fig. 6.3. 



Comparing to experimental and lattice data. Since the nucleon mass as a 
result of the quark-diquark calculation usually deviates from its experimental value 
-^exp = 0.94 GeV, some precaution must be taken when comparing form factor results 
to experimental data. There are basically two scales which enter the electromagnetic 
current and form factors: the photon momentum and the nucleon mass M, the 
latter of which absorbs both the inherent scale Air and the current-quark mass de- 
pendence of the system. As form factors are dimensionless they can only depend on 
the combination Q'^/M'^. This is accounted for in Figs. (6.3-6.6) where the results are 
plotted as a function of this dimensionless variable. 

The same effect is achieved by an appropriate rescaling of the photon momentum 
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Figure 6.6: Contributions to G^; and F2 according to the type of incoming and 
outgoing quark-diquark amplitudes ^. The form factor sums (SUM) are decomposed 
into scalar-scalar (S), axial- vector-axial- vector (A) and the sum of scalar-axial- vector 
and axial- vector-scalar contributions (M). 



Presenting results as a function of removes the trivial scale dependence induced 
by the overestimated quark-core result for M. When investigating the current-mass 
dependence and comparing to lattice-QCD results, the experimental mass in Eq. (6.34) 
should be replaced by that obtained from the lattice: Mexp — > M\a,t. This affects the 
comparison of charge radii as well: instead of comparing the dimensionless combinations 
rpM'^ (e.g., as in Ref. [113]) we rescale the radii accordingly, 

(^'^')calc = ^^l'at> (6-35) 

and plot the quantities instead of r^. These issues were not relevant in the context 
of pion observables /tt, in Chapter 3.3 since the mass of the pion is not expected to 
drastically change upon implementing beyond rainbow-ladder contributions. 

A second remark concerns the comparison of magnetic moments to lattice-QCD 
results where F2 and are usually expressed in terms of experimental (i.e., fixed in- 
stead of current-mass dependent) nuclear magnetons [59] . The corresponding quantity 
which appears in the electromagnetic current (6.7) is i<2/(2Mexp). Again, our calcu- 
lated result for F2 (the value in units of 'running' magnetons) is dimensionless, hence 
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we leave it unchanged but remove the current-mass dependence when comparing to 
lattice data via _ 

where M^^^^ is the calculated nucleon mass at the physical u/d point. 

Magnetic moments. The results for the nucleon's magnetic moments at the phys- 
ical u/d mass are collected in Table 6.2. The current-quark mass dependence of the 
anomalous magnetic moments Kp, k„ is shown in Fig. 6.7. At large quark masses, 
where possible pion-cloud effects should be diminished, the result still overestimates 
the lattice data. This suggests further missing corrections beyond pionic contributions. 

In analogy to the discussion of the nucleon mass in Section 5.3, the pion-cloud 
effect to magnetic moments can be quantified by examining the loop contributions in 
heavy-baryon chiral effective field theory: 

fip,n{ml) = f,;%%ml,A) + fi^^imlA) . (6.37) 

Retaining only those loop diagrams which describe a direct coupling of the photon to 
the intermediate pion, where the accompanying baryon is either a or a A, yields 
[242,243]: 



oo 

S , 2 *^ , AMjv f , 
^P,n("^^>A) = ±^-2^2 J dx 





^4 



■ 2 2 / ^ UJ{UJ + 26) 2 r ^ 



(6.38) 



where C = —2- (0.76), and the remaining quantities are explained in the context of 
Eq. (5.18). Using the same input as there, i.e. a dipole regulator unn[x) = una[x) = 
l/(l-|-x^)^ with A = 0.8 GeV, yields at the physical point: /Up,j = ±0.61. The inclusion 
of further meson- loop diagrams can diminish this value [243] . 

Since the simplest pion-loop contributions to proton and neutron carry an opposite 
sign, their total cancels in the nucleon's isoscalar and doubles in its isovector combina- 
tion, defined by 

This implies hence isoscalar quark-core quantities should be comparable 

to their experimental or lattice counterparts. The experimental isoscalar magnetic 
moment is small and negative: Kg = —0.12. Our result at the u/d current mass is 
Kg = —0.03(3); it is consistent with zero throughout the calculated pion-mass range. 

Electric and magnetic radii. Table 6.3 shows the results for the charge and magnetic 
radii together with the experimental values. In agreement with the interpretation in 
terms of an hadronic quark core, and similar to the result for the pion charge radius in 
Fig. 3.7, the radii are sizeably underestimated. 

By virtue of Eq. (6.39), the isoscalar and isovector Dirac and Pauli radii read 
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Figure 6.7: Anomalous magnetic moments of proton and neutron vs. pion mass. The 
results of setup (C3) are compared to quenched lattice data and their chiral extrapo- 
lations {dashed curves) [59]. Dots denote the experimental values. 









Kg 


(C3) 


2.56(5) 


-1.58(3) 


-0.03(3) 


Exp. 


2.79 


-1.91 


-0.12 



Table 6.2: Results for proton, neutron and isoscalar magnetic moments in setup (C3), 
compared to experiment. 





p 




' M 




(C3) 


0.79(2) 


0.00(1) 


0.73(2) 


0.72(2) 


Exp. 


0.89 


-0.12 


0.86 


0.87 



Table 6.3: Electric and magnetic radii of proton and neutron in setup (C3) compared 
to experimental values. All units are fm except (r^)^ which is fm^. 
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Because of the smallness of Hs both experimental and calculated values for r| suffer 
from a large statistical uncertainty. Lattice studies usually compute isovector form 
factors where numerically expensive contributions from topologically disconnected di- 
agrams cancel. Pion-cloud effects are increased in the isovector channel; corresponding 
lattice results are compared to the quark-diquark calculation in Fig. 6.8 where a mutual 
agreement at large pion masses is clearly visible. 

In contrast to the magnetic moments, the pionic components of electric and mag- 
netic charge radii r|., r^^ which correspond to intermediate A^vr and Avr states carry 
opposite signs in the flavor-separated channels, i.e., for the u and d-quark contribu- 
tions [246]. The pion contributions to the charge radii logarithmically diverge in the 
chiral limit since a massless pion has an infinite Compton wavelength. The quark- 
diquark result yields (r^)^ « throug 
and, using Eq. (6.29), therefore implies: 



diquark result yields (r^)^ « throughout the inspected current-quark mass range 



{r^Ef - {4)' - {4)"- (6.41) 

A chiral expansion for the electromagnetic radii of proton and neutron similar to (5.22) 
can be found in [57]. 

As discussed in Section 6.1, the Dirac form factor Fi(Q^) is the Fourier transform 
of the nucleon's transverse charge distribution in the infinite-momentum frame [224] , 
hence the corresponding Dirac radii feature a direct interpretation in terms of transverse 
charge radii. Experimentally: FJ^{Q'^) < 0, and the neutron's Dirac radius is positive: 
(rj*)^ > 0. With the convention of Eq. (6.30), where denotes the radius of the up- 
quark contribution in the proton and, via charge symmetry, that of the down-quark in 
the neutron, Eq. (6.29) implies that the charge radius of the d-quark in the neutron is 
smaller than that of the u quark: 



2 



3 



{r^f >0. (6.42) 



The extracted transverse charge densities confirm that the central charge distribution 
of a fast-moving neutron (proton) is negative (positive) [222]. This result is at odds 
with the traditional view of a zero-charge neutron whose pir^ pion-cloud component 
generates a non-zero charge distribution which has a negative long-range tail but is 
small and positive at its core. 
The negativity of Ff (Q^) 

is a natural feature of a quark-diquark model [63, 113]: 
it can be explained by axial-vector dd diquark correlations (corresponding to s'^ in 
Eq. (A. 35)) which, other than the scalar-diquark ud and axial-vector ud contributions, 
induce a localization of the d-quark in the neutron. The quark-diquark model result 
for the scalar-diquark contribution to (r")^ at the light quark mass is 0.00(1) fm-^ 
(Fig. 6.9). Adding axial-axial and scalar-axial correlations yields (r")^ = 0.11(1) fm^ 
which basically cancels the neutron's Foldy term to obtain ~ 0, cf. Eq. (6.22). This 
result is large compared to the experimental value (r")^ = 0.01 fm^ and indicative of 
further destructive interference with pion-cloud corrections in the axial-vector diquark 
channel. 
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Figure 6.8: Squared isovector radii corresponding to the Dirac and Pauli form factors 
Pi 2 = Pi2 - ^12 compared to lattice results [244,245]. 
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Figure 6.9: Left panel: Electric charge radii for proton and neutron vs. pion mass- 
squared. Right panel: Pion-mass evolution of the neutron's Dirac charge radius. The 
full result is compared to the scalar-diquark contribution, cf. Fig. 6.6. Dots with 
statistical errors denote the experimental values. 
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Figure 6.10: Q^-evolution of the proton's form factor ratio ^Ge/Gm in setup (C3), 
where the band denotes a variation of u). The result is compared to experimental 
data obtained from Rosenbluth-separation [247] and polarization-transfer measure- 
ments [229,230,248]. 



Large-Q^ behavior. At larger momentum transfers ^ 2 GeV^, pion cloud effects 
must vanish since the structure of the nucleon at small distances (< 0.15 fm) is probed. 
The large-Q^ behavior of the form factors therefore reflects genuine properties of the 
nucleon's quark core. Since quark and diquark propagators obtained from the quark 
DSE or via the T-matrix ansatz necessarily exhibit singularities, the applicable 
range of the form factor calculation is bounded from above (see App. B.3) and typically 
limited to few GeV^. Access to the large-Q^ domain, where the rainbow-ladder result is 
expected to become increasingly accurate, necessitates appropriate methods to evaluate 
these Green functions beyond their dominant singularities, i.e., to include the respective 
residue contributions in the form factor integrals. 

The solution for the proton's form factor ratio /XpG^/G^^ is shown in Fig. 6.10 and 
comparable to the quark-diquark model results of Refs. [62,63]. The small-Q^ structure 
is a further indication of missing pion-cloud effects. A Taylor expansion at small 
entails 




(6.43) 



Experimentally: ~ r^, whereas the quark-diquark calculation yields > r^. 
Pionic contributions affect electric and magnetic radii differently [63]. 
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The evolution beyond ^ 2 GeV^ is mainly dominated by the proton's quark 
core. We note that the phenomenological transverse term introduced in the seagull 
vertices, Eq. (A. 138), has a sizeable impact on the form factor ratio in that region. The 
parameters (A. 140) optimize mutual agreement with the polarization-transfer data for 
Ge/Gm and those for G^j and G^^. Nevertheless a variation of these parameters 
only minimally affects the small-Q^ region which therefore only depends on the input 
introduced in connection with the effective coupling a{k'^) of Section 3.2. We note 
that the effective diquark vertices would naturally no longer appear in a form factor 
calculation within the three-body approach according to Eq. (6.16). 



Chapter 7 

Conclusions and outlook 



The covariant bound-state formalism, based on the Dyson-Schwinger equations of QCD, 
provides a versatile tool for the calculation of hadron properties. It is formulated within 
QCD, fully relativistic, and represents a continuum approach; it provides access to 
infrared and ultraviolet momenta; and, with appropriate numerical algorithms which 
have become available in the past years, it covers the full quark mass range from the 
chiral limit up to arbitrarily large current masses. At the present stage, especially in 
the baryon sector, we are only beginning to explore its potential and possibilities. 

In this thesis a study of pion, p-meson, nucleon and A properties was presented in 
such an approach, where the inherent link between meson and baryon observables is 
expressed through a rainbow- ladder truncation in the quark-antiquark and quark-quark 
channel. The effect of dynamical mass generation is manifest, while the truncation at 
the same time preserves the nature of the pion as a Goldstone boson of spontaneously 
broken chiral symmetry. 

The rainbow-ladder kernel embodies an effective quark-gluon coupling a{k'^) which 
represents the only phenomenological input in the calculation of meson and baryon ob- 
servables. Two models for the current-quark mass dependence of the coupling strength 
were compared: the first setup is defined by a current-mass independent strength which 
has long been known to provide a good description of pseudoscalar and vector-meson 
properties; in the second setup, chiral corrections — partly owing to pseudoscalar 
meson-cloud effects — were anticipated to construct an overestimated quark core for 
the p-meson. The distinction between the two models arises in the infrared part of the 
coupling, and its impact upon hadronic observables turned out to be an approximate 
additive contribution to quantities with dimension of a squared mass. 

A self-consistent solution of the three-body equation for the nucleon was presented 
for the first time, where the kernel was given by a rainbow-ladder gluon exchange 
between any two quarks. The result was compared to a quark-diquark calculation 
which is enabled by the occurrence of artificial diquark mass poles in the qq scattering 
matrix induced by the gluon ladder kernel. The result for Mj\f{m'^) is very similar in 
both approaches; the quark-diquark model is slightly too attractive, i.e. by 50 — 100 
MeV. The overall properties of an inflated quark core are clearly reproduced in Mjv 
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and Ma; in particular, agreement with lattice data is obtained for all investigated 
hadron masses when the quark mass becomes large. The mass of the A baryon in the 
quark-diquark calculation appears to be unexpectedly large; it remains to be seen if the 
three-quark approach provides further attraction or, as in the nucleon case, moderate 
repulsion. 

Results for the nucleon's static electromagnetic properties as a function of the pion 
mass were presented and compared to lattice results and chiral extrapolations. Simi- 
larly to the pion charge radius, good agreement was found for the nucleon's isovcctor 
radii above the strange-quark mass whereas missing pion-cloud effects arc clearly visible 
in the chiral region. Except for G^, the overall consistency of the nucleon's electro- 
magnetic form factors as a function of the photon momentum with experiment is quite 
remarkable. Pionic effects are missing at small Q^, and their absence is amplified in the 
proton's form factor ratio Ge/Gm- The implementation of self-consistently obtained 
quark and diquark propagators with a complex singularity structure leads to an upper 
limit at several GeV^ ; together with a model dependence of the diquark ingredients at 
larger photon momentum it impedes an unambiguous signal of a possible zero crossing. 
The latter obstacle can be overcome by a form factor calculation within the three-quark 
approach. 

It is imperative to study the effects of different interactions in the qq and qq channels: 
in particular, further admixture of a scalar-scalar interaction, induced by a scalar part 
in the quark-gluon vertex, is expected to yield notable ramifications for the infrared 
structure of QCD and derived hadron properties; the same is true for pionic effects. 
Nevertheless, with a solution of the three-quark equation available, the journey beyond 
rainbow-ladder is now simultaneously and consistently possible in both meson and 
baryon channels; and the covariant bound-state framework as an ab-initio approach 
evolves one step forward towards an understanding of hadron dynamics. 

The now achieved increase of predictive power finally allows to reach out for more 
observables which are in the focus of current experimental interest. For instance, a cal- 
culation of the A mass, electromagnetic form factors and A'"-A transition form factors 
is a natural application of the three-body approach. Further possible future directions 
are: an extension to the heavy-quark regime, including hadrons with open strange, 
charm and bottom quantum numbers; radial excitations of baryons and a clarifica- 
tion of the role of the Roper resonance; strong and electroweak scattering processes; 
form factors at large Q^; a study of exotic mesons in the bound-state formalism; and 
a description of the hadron's internal momentum and angular momentum structure 
expressed through parton distribution functions and GPDs. 
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A collection of propagators, 
vertices and amplitudes 



In this appendix we collect the basic Green functions and bound-state amplitudes which 
appear in the main text. We state their general properties and specify the expressions 
that are used in the numerical computation. The quark propagator (A.l) and quark- 
photon vertex (A.2) represent the 'elementary' quantities of the approach. Color-singlet 
mesons and colored diquarks are composite objects; their on-shell and effective off-shell 
properties are summarized in (A. 3) and (A. 4), respectively. We discuss the bound-state 
amplitude for the nucleon in the three-body framework (A. 5) and that for N and A 
in the quark-diquark approach (A. 6). The diagrams which contribute to the nucleon's 
electromagnetic current in the quark-diquark model are stated in (A. 9); they involve 
diquark-photon vertices (A. 7) and seagull amplitudes (A. 8). 



A.l Quark propagator 

As a fermionic two-point function which involves one momentum p, the quark propa- 
gator can only depend on the two Dirac structures 1}. The corresponding dressing 
functions and as may be expressed through the quark renormalization function Zj 
and the renormaliization-point independent quark mass function M: 




Another frequently used notation involves the quantities A{p'^,fi^) = and 
B{p'^,fi'^) = M{p'^) /Zf{p'^ , fj?). The inverse propagator reads: 

= Aip^^i^) {z^ + Mip')) . (A.2) 



A.l Quark propagator 
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Asymptotic behavior and quark condensate. Asymptotically, the DSE solution 
for the quark mass function reproduces the behavior predicted from perturbation theory 
(see, e.g., [98,193]): 

, 2^ P^'-oo m ^ 27r^7^ -{qq) 



where T^p^) = ^In yp^ /Aqq^J . The coefficients m and —{qq) defined thereby are 
the renormalization-point independent current mass and chiral condensate. For finite 
current masses, the second term is suppressed by a factor of while in the chiral limit, 
defined by m = 0, it determines the behavior of the asymptotic mass function. 

The cutoff-dependent bare mass mo that appears in the bare quark propagator 
and enters the quark DSE (3.2) is related to the renormalized mass via the mass 
renormalization constant Zm- m-o(A^) = Zm(/i^, A^) m^. For a large renormalization 
point, M{fj?) can be identified with which entails mo(A^) = M(A^) [98]. In the 
following we will always assume to be large, i.e. fi ^ Aqcd; in our calculation we 
use the value fi = 19 GeV. 

The renormalization-point-dependent chiral quark condensate is obtained from the 
trace of the chiral quark propagator: 

A 

- {qq)i, = Z^ifi^A^) Zmifi^A^) Nc J Tin{S,i,i,,i{q, /i)}, (A.4) 

with Zm(/i^, A^) = M [h?) / M [fj?) evaluated at large current masses. For a large renor- 
malization point, the current-quark masses and condensates are related via 

^(^') = :^(^' -{m), = -{m)H^^'V-- (A.5) 

Solving the quark DSE. The quark DSE (3.2), 

S{p, fi)-' = Z2{fi\A^) {ii> + M(A2)) + S(p, fi, A) , (A.6) 

can be rewritten in terms of two coupled integral equations for the quark propagator's 
dressing functions A{p'^,fi^) and M{p'^): 



(A.7) 



A{p\fi^)= Z2{fi\K^) + ^A{p\fi\^^), 

M{p^)A{p\fi^) = M(A2) Z2{fi^A^) + SM(p^/x^ A2), 

where the scalar functions and T,m constitute the quark self-energy via 

S(p,^,A) = i;^SA(p^/i^A2) + SM(/,/U^A2). (A.8) 

Eqs. (A.7) can be solved iteratively for chosen values of Z2 and M(A^). Upon employing 
a renormalization condition, e.g. A{fi'^, fi^) = 1, and specifying the current mass M^ji^) 
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at the renormalization point as an input parameter, both Z2 and M(A^) are determined 
together with A(p^) and M{p'^) in the course of the iteration via 

Z,(;.^A^) = 1-S^(;.^^^A2), M(A^) = ^(^')-^||;^g^^'^ (A.9) 

Omitting the renormahzation-point and cutoff dependence for brevity, the self- 
energy integral (3.3, 3.7) in rainbow truncation reads 

A 

np) = jn''9{k')rS{q)^\ g{k'):=Zl^-f^, (A.IO) 
<? 

where k"^ = + (f' — 2p ■ q =: + (f' — 2 \/p^ \fq^ z is the squared gluon momentum 
and a shorthand notation for the effective coupling. The self-energy coefficients 

become 



E^(p2) = y a,{q^)g{k^)F{p\q\z), Sa/(p^) = 3 / oM^)9{^^) (A.H) 
9 



and involve the quark dressings Ogijf') which depend on Ai^q^) and M(cj^). The 

dimensionless quantity F is given by 



2 2 
= P-Q + -^{P- k){q -k) = 2,p- q - ^{p^ q^ - {p- qf) = (a.12) 



= ~k + = P +3p-fc + 2— 

Since the only Lorentz-invariant combinations which appear in the quark self-energy 
(A. 11) are p^, and z, the loop integral becomes two-dimensional (cf. Eq. (B.6)), e.g.: 



a,iq^)g{k^)F{p^q'',z)r^ j dq\^ a,{q^) j dz^J I - z^ g^k") F[p\q\ z). 

q 0-1 

For a spacelike external momentum p^ G R-i-, the squared gluon momentum k'^ is real 
and positive as well and the coupled system (A. 7) can be solved without complications. 

Quark propagator in the complex plane. If the external quark momentum p^ is 
complex, the argument /c^ of the coupling gik"^) constitutes the interior of a parabola 

fc2 =p2 _^g2 _^2 v^V^= (^^±i|Im\/p2|)^ tGE+, (A.13) 

which passes through the outer point p^ and has its apex at (cf. Fig. A.l) 

t = Q =^ kl = - (lm^/^y = ^ (Rep2 _ \p^\) . (A.14) 



A.l Quark propagator 
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Ug{k^)F{p^,q^,z) exhibits non-analyticities in that k"^ domain, one must either adjust 
the size of the parabola (A. 13) by imposing a hmit on j?^ or resort to refined numerical 
methods. The kinematic singularity in F at A:^ ^ induced by the transverse gluon can 
be compensated by a vanishing g{k'^ 0). The infrared behavior of the parametrization 
(3.14) almost satisfies that criterion: 

-p > -j^, Ao = 1 GeV. (A. 15) 

The remainder stems from the ultraviolet term and is relatively small compared to the 
overall strength of g{k'^). It results in small numerical artifacts which are visible in the 
complex functions (Jv{p^) and (Jsip^) upon a straightforward integration. 

The complex conjugate branch points in the logarithmic tail of a(A;^) lead to a 
theoretical limitation 

ImVp2)^ > -Kl,cD (^Im'\/-1± Ve2-2^ = -(0.31 GeV)^ (A.16) 

which is practically unimportant since they are concealed by the large oscillations of 
the exponential parts in a{k'^). The complex conjugate poles of the resulting rainbow- 
ladder quark propagator are insensitive to these singularities as well. 

A coupling which does not satisfy g{k'^ ^ 0) — > or involves further singularities 
in the integration domain inevitably requires an advanced numerical treatment. The 
same is true for a truncation beyond rainbow-ladder which leads to a more complicated 
structure of the self-energy integrals. We shortly discuss two such strategies: 

• Complex rays. Upon performing all the integrations except the integral, 
singular points in /c^ lead to branch cuts in the complex plane, as illustrated 
in Fig. A.l. For instance, a pole at A;^ = generates a circular branch cut in 
the complex q^ plane with an opening at q^ = (dashed line). The original 
integration contour G R crosses these branch cuts and the resulting numerical 
artifacts become dominant in a straightforward integration. 

The logarithmic one-loop behavior (3.13) entails that any coupling exhibits sin- 
gularities at 7^ as well: those will generally lead to more complicated branch 
cut structures (dash-dotted line) which however still leave the arc q^ = re*^''^^ , 
r G 11+ , unharmed. A possible way to avoid all occurring branch cuts is to de- 
form the integration contour q^ G (0,A2) to a complex arc that passes through 
the point and eventually returns to A^ G IR, in the far spacelike region. Since 
(Jy{q^) and CTs{q^) must already be known on these complex paths, the complex 
DSE solution is therefore obtained via iteration of (A. 7) on a family of deformed 
complex paths in p^. 

This method has been sketched in [129] and is ideally suited if the singularities 
of the resulting quark propagator which are generated during the iteration only 
appear on the timelike p^ axis. Such a singularity structure appears if a Ball- 
Chiu-like ansatz for the quark-gluon vertex is applied [45]. 
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Figure A. 1 : Left panel: Integration domain of quark (q'^) and gluon propagator (/c^) in 
the self-energy integral corresponding to a fixed external point p^. Right panel: Branch 
cuts in the complex plane of the quark propagator and a possible integration path. 



• Complex parabolas. A different strategy to avoid singularities in the inte- 
gration domain has been described in Refs. [110, 128,249]. The basic idea is to 
use the gluon momentum k instead of the quark momentum q as the integra- 
tion variable^. As a consequence, the argument k'^ of the coupling becomes real 
and complex; hence the only relevant singularities are those of the final quark 
propagator which are generated during the iteration. Now the internal quark 
momentum q"^ = p"^ + k^ + 2 p ■ k, with /c^ E R+, is bounded by a parabola 

q%=p'^ + k'^ ±2y^Vk^ = (^t' ±i\lmy^\y , G 1R+, (A.17) 

whose position is characterized by the external quark momentum p'^. The itera- 
tion in (A. 7) takes place on such a parabola, and the Cauchy formula is used to 
obtain the quark propagator in the interior at each iteration step. 

The method is particularly useful if the resulting quark propagator exhibits com- 
plex conjugate poles, a feature which is common in rainbow-ladder studies and 
was also recovered in more general truncations [42, 128]. It furthermore allows for 
a more accurate determination of the complex-plane propagator at larger quark 
masses [249]. 

Analogous procedures can be implemented for evaluating meson and diquark ampli- 
tudes in the complex plane of the relative momentum between the contributing quarks. 

^ The change of integration variables requires a translationally invariant regularization of the self- 
energy integral. This is not realized by a hard cutoff A which leads to (however small) deviations 
between results obtained in both methods. Implementing a Pauli-Villars regulator 1/(1 + / h?) 
attached to the gluon propagator suppresses the UV modes and restores translational invariance but 
ignores the perturbative 1-loop behavior of the effective coupling. 
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A. 2 Quark-photon vertex 

The quark-photon vertex is the central ingredient of any hadronic electromagnetic form 
factor diagram since, of all the fundamental degrees of freedom in the QCD Lagrangian, 
only the quark is electrically charged. With the notation k = (/c+ -|- k-)/2 (relative 
or Breit momentum) and Q = k+ — k- (incoming photon momentum), where k+ and 
k- are outgoing and incoming quark momenta, the general form of the quark-photon 
vertex is 



12 
i=l 



where a possible representation of the Dirac basis elements is given by 

rf (A:, Q) e {7^ Q^} x ft, Q]} . (A.19) 

Induced by electromagnetic current conservation, a vector Ward-Takahashi identity 
constrains the longitudinal contributions ~ of the vertex by relating them to the 
quark propagator: 



Q'^T>'^^^ik,Q) = S-\k+) - S-\k^) 



(A.20) 



Here we used the abbreviations 



F{kl)^F{}l_) F{kl)-F{k^_) 

■= ^ , .= ^2 _fc2 ' (^-21) 

where A{p'^), BljP') are the quark propagator's dressing functions, and for 0: 
F(A;2) and A^;' F'{k'^). The differential Ward identity for ^ reads 



Tfa^C^'O) = —77^ = ^rA[k^) + 2A:^ (iM'(fc') + B'{k^)) . (A.22) 



Implementing both relations yields the most general expression for the quark-photon 
vertex: 



rfq)(A;,Q) =i7^SA + 2F(i^AA + AB) + r^'^r^(A;,Q). (A.23) 



The first part is the Ball-Chiu vertex [250]. The transverse contribution Tj.{k,Q) is 
constructed from the eight basis elements ~ j^,k'^ and must satisfy r^(A;,0) = 
to obey the Ward identity, either by a dependence of the basis elements or by a 
vanishing amplitude at = 0. 
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Figure A. 2: Momentum routing in the quark-photon vertex. 



We note that, due to the transversahty of the photon, the purely longitudinal com- 
ponents of the vertex do not contribute to any hadronic matrix elements: only the 
transverse projection of (A. 23) does. In this sense the WTI alone provides no constraint 
on physics. If, according to the basis (A. 19), one starts from the general expression 

rfq)(A;, Q) = T^" {iY T^'Hk, Q) + F T^^^k, Q)) + T^'^k, Q), (A.24) 

then r*^-*^) and T^^^ determine the physical content, while only the unphysical component 
r*^^) is constrained by the WTI: 

Q^r(3)(fc,Q) =L^'^ (^iY^A + 2k''{i^AA + AB)y (A.25) 

The transverse and longitudinal parts (i.e., the brackets in Eqs. (A.24) and (A.25)) must 
however share the same limit at — > in order to avoid a kinematic singularity. This 
is ensured by the differential Ward identity which subsequently determines Eq. (A. 23). 

Modeling the transverse part. The transverse part must vanish for = due 
to current conservation. In the perturbative limit ^ cxd it is C'(a)-suppressed 
compared to the Ball-Chiu construction which ensures a bare vertex via Z2 and 

Aa^b 0. Several parametrizations for the transverse part have been devised in the 
literature, amongst which is the Curtis-Pennington ansatz that ensures multiplicative 
renormalizability in the Dyson-Schwinger equations of quenched QED [251]: 

rf^{k, Q) = zAa ^ ({k ■Q)Y- k^ q) . (A.26) 

(A: • Q)^ + (SmO ^ ^ 

A modified ansatz which accomplishes the same constraint in unquenched QED was 
recently proposed [252]; see Ref. [253] for a comprehensive overview on the topic. 

A self-consistent solution of the full quark-photon vertex is enabled by its inhomo- 
geneous Bethe-Salpeter equation [68], 

rfq)(A:, QU = Z2 il^ + j Kik, k', Q)^^^sp {s{k'^) T^^^p\ Q) S{k'^)} , (A.27) 

k' 
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which, for consistency with the rainbow-ladder approach, needs to be solved with a 
gluon-ladder kernel as its input. The resulting vertex [83] self-consistently generates a 
vector-meson pole at = — whose contribution significantly increases the charge 
radii of pseudoscalar and vector mesons [83,86]. On the domain — < < 0.2 GeV'^ 
it could be described by a Ball-Chiu vertex together with a phenomenological ansatz 
owing to the yO-meson pole for the transverse part. We do not intend to repeat that 
calculation here and hence adopt a similar ansatz: 

r^(fc, Q) = -1 ^ e-^(-) r(:,(fc, Q), (A.28) 

Qp X+l 

where x = Q"^ /m?p and gp = ^/2mp/ fp. Eq. (A.28) necessitates knowledge of an (un- 
physical) off-shell />meson amplitude T';^^(k,Q). On the mass- sheh (Q2 = -mf,) it 
is determined self-consistently from its homogeneous BSE; for general we use the 
off-shell prescription stated in App. A. 4. 

To compensate for any arbitrariness from the off-shell part, we introduced the func- 
tion e~^^^^ where the choice g{x) = {pi -\- p2x'^){l + x) optimizes agreement with 
the results for the quark-photon vertex and the pion charge form factor obtained in 
Ref. [83] within setup (CI). The corresponding value for is reproduced if a current- 
mass dependent parameter g{0) = pi = 0.001 -|- m|/(3.72 GeV^) is chosen^ [113]. The 
remaining parameter p2 is relevant for the medium-Q^ evolution of the p-meson part 
in the vertex and impacts upon the nucleon form factors of Section 6.3. The value 
P2 = 0.001, together with the transverse part in the seagull amplitudes (App. A. 8) 
maximizes agreement with the polarization-transfer data for the proton's form factor 
ratio Ge{Q^)/Gm{Q^). 

We repeat that the ansatz (A.28) is inspired by low-Q^ phenomenology which, at 
this point, is the only accessible domain in our form factor calculation due to singularity 
restrictions (cf. App. B.3). An exponential suppression of the transverse part clearly 
disagrees with the large-Q^ analysis, and future form factor investigations in that region 
inevitably require an omission of the above ansatz in favor of a self-consistent solution 
from the inhomogeneous BSE (A. 27). 



Note that the different values for pi in Refs. [60] and [113] were obtained through a shghtly 
different off-sheU ansatz for the p-meson amplitude. 
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A. 3 Meson and diquark amplitudes 

Meson amplitudes. The meson bound-state amplitude has been introduced in con- 
nection with its homogeneous BSE (3.9). The Dirac structure of the amplitudes is 
determined by the Clifford algebra of the gamma matrices. A general fermion- scalar 
vertex depending on the momenta q and P allows for 4 basis elements 

ri(g,P)e{l,f, ^, [^,f]}, (A.29) 

and a fermion- vector vertex includes 12 structures: 

T>t{q,P) G {l^q^P^} X {l,f, I [If]}. (A.30) 

The negative parity requirement for pseudoscalar and vector meson amplitudes, 

^(g,P) = -7^^(Ag,AP)7^ 
T^{q,P)=^^K^''T%Kq,kP)-i^ ^ ' ' 

with the parity transformation A = diag(— 1, —1, —1, 1), requires the inclusion of a 7^ 
matrix. The resulting amplitudes, written with full Dirac, color and flavor dependence, 
are given by: 





= Y^fl\q\z,P^) 










k=l 






(A.32) 




= Y,fr{q\z,P') 
k=l 


{<{,,F)}^^« 





The dressing functions fk{q'^,z,P'^) depend on the Lorentz scalars g^, P^ and the 
angular variable z = q-P. Upon solving the meson Bethe-Salpeter equation, they are 
obtained on the domains G 1R+ (by implementing the same methods as discussed in 
App. A.l also for q"^ G C), P^ = — (i.e., on the mass shell), and z G (—1, 1). Greek 
indices refer to the Dirac structure, and the color structure of the meson amplitudes is 
diagonal {A, B = 1,2,3). 

The isospin-triplet flavor matrices r°, corresponding to the states tt^, tt^ and p^, 

are given by 



r+ = ud"!" = i (cti + ia2 
r~ = du"!" = 5 (cTi — ia2 



w • N 71 ("^ '^'^ ) = 71 , (A.33) 



where the dj are the Pauli matrices and u = (1, 0), d = (0, 1). They are normalized to 
unity: Tr{r^tre'} = See'- 
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Figure A. 3: Notational conventions for diquark amplitudes. The quark momenta are 
q± = ±q + P/2 and the indices are Dirac/Lorentz, color and flavor indices. 



Diquark amplitudes. Apart from opposite parity, scalar and axial-vector diquark 
amplitudes exhibit the same Dirac structure as their pseudoscalar and vector meson 
counterparts. The incoming antiquark momentum — q_ is replaced by an outgoing 
quark momentum which is reflected by the charge conjugation matrix C = 7^7^. 
We denote diquark amplitudes by the same generic symbol T: 



4 

T{q,P)=^fnq',z,P')[ij'rk{q,P)c}^^ ^ £^ ^ s° 



fc=l 

12 



(A.34) 



r'^{q,P)=j:m<l',z,P^)hri:{q,P)c} ® ^ C3 sl'^^' 



The isospin singlet and triplet matrices s^j^ and s^'^'^ for the diquarks read 

= Uut = 1 (l + fjs) , 

= ^ (udt - dut) = ^ , = ^ (udt + dut) = -Ij ai , (A.35) 

S3 = ddt = i (1 - CJ3) 

with the Pauli matrices 0", and u = (1,0), d = (0,1). They are again normalized via 
Trjs^^s^ } = See'- Due to the Pauli principle, diquark amplitudes must be antisymmetric 
under quark exchange g+ <-> 

r(5,p) = -r^(-gU„(i_^),p), (A.36) 

where the transposition involves all Dirac, color and flavor indices. Because of the anti- 
symmetry of the color anti-triplet diquark the combination of flavor and spin structure 
must be symmetric. As a consequence, spin and isospin states coincide for the two-flavor 
case: scalar diquarks correspond to an antisymmetric isospin singlet and axial-vector 
diquarks to a symmetric isospin triplet. 
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While the flavor matrices do not give a contribution to the BSE integral (5.6), 
the color factor ~ eabe representing the diquark anti-triplet configuration leads to a 
prefactor 1/2 compared to the meson BSE: diquarks are less bound than mesons. The 
Dirac amplitudes ~ 1, 7^ in (A. 29) and (A. 30) are the dominant ones in a solution 
of the rainbow-ladder BSE for the lowest-mass mesons and diquarks and reproduce 
masses of the full solution within an error of < 20% [78, 79, 187]. 

Charge conjugation, C parity and Pauli principle. The charge-conjugated Dirac 
amplitudes are defined by 

T>^{q,P):=-CT^^'^{-q,-P)C^, 

where the superscript T denotes matrix transposition'^. Pseudoscalar and vector mesons 
with equal quark masses are C-parity eigenstates (e.g.: J^^ = '", 1 for vr and p) 
which entails 

T{q,P)=T{q,-P), f^iq, P) = F^'iq, -P) (A.38) 

for the Dirac part of the ground-state amplitudes, assuming equal momentum parti- 
tioning. On the other hand, diquarks are subject to the Pauli principle which leads 
to 

r(g, P) = -r^(-g, P) , T^^iq, P) = T'^^'i-q, P) . (A.39) 

For an appropriate Dirac basis whose elements satisfy Eqs. (A.38-A.39), both con- 
straints require the dressing functions to be even in the angular variable z = q ■ P. If 
a Chebyshev expansion is employed (for details, see App. B.2), only even Chebyshev 
moments contribute which is ensured by the fact that the BSE kernel decouples even 
and odd Chebyshev moments. Conversely, pseudoscalar and vector-meson states with 
opposite C— parity (associated with the 'exotic' quantum numbers 0~~ and 1 '") carry 
an odd dependence on z. 



Mesons and diquarks on the mass shell. For the actual solution of the meson 
and diquark BSEs, Eqs. (3.9) and (5.6), it is advantageous to construct orthogonalized 
versions of the general Dirac basis elements (A. 29, A. 30) at the respective mass poles 
P^ = —M"^. Suitable choices for the (pseudo-)scalar and (axial-)vector cases are given 
in Table A.l. Their orthogonality relations read 

^ Ti{Ti T,} = 5i, a,{z) , ^ T^"" TV{rf r/} = <5,, bi{z) , (A.40) 

^ Note that this definition differs in two respects from the somewhat standard convention in the 
hterature. First, the left-hand sides of Eq. (A. 37) are usually denoted by T[q, —P), whereas the above 
definition is more convenient when using the abbreviation V. Second, the conjugated vector-meson 
amplitude is usually defined without the minus sign. We chose Eq. (A. 37) to enable a common definition 
of charge conjugation for meson and diquark amplitudes, given that the dominant amplitudes in both 
cases carry a factor i, i.e. 
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T n 

J — u 


J — i 


Ti = 1 

r2 = r 

T3 = Z§T 

r4 = i[ir] 


Tr^ = iq^ Ty = q'^^T — I qrl^ 



Table A.l: Orthogonal basis elements for pseudoscalar mesons/scalar di quarks 
(J = 0) and vector mesons/axialvector diquarks (J = 1), designed such that all corre- 
sponding dressing functions fkiq'^, z, — M^) are real and even in z for equal momentum 
partitioning, qt^ = Tp'^q'^, where Tp'^ is the transverse projector with respect to P. Due 
to the transversality of on-shell vector mesons (and, by analogy, axialvector diquarks) 
an additional projector Tp"^ must be attached to the J = 1 basis elements; hence the 
four longitudinal basis elements ~ do not contribute on the mass shell. 



where 

ai = a2 
as 

Contraction of the BSEs with these basis elements and exploiting the orthogonality rela- 
tions leads to coupled homogeneous integral equations for the coefficients fk{q^i z, — M^) 
which can be solved in any arbitrary frame. Exploiting the 0(4) symmetry of the 
problem, a decomposition into Chebyshev polynomials of the second kind is usually 
employed for numerical convenience (see App. B.2). Typically only a few Chebyshev 
moments fj^{q^) have to be taken into account to match the full solution [78,79]. This 
observation will be used in the context of the diquark amplitudes' off-shell continuation 
where the angular dependence is neglected (App. A. 4). Upon introducing an artificial 
parameter A(p2) in the homogeneous BSE (5.6), the equation becomes an eigenvalue 
problem where a bound-state solution is obtained for A(— M^) = 1. 

Normalization. A bound-state amplitude is normalized by a canonical normaliza- 
tion condition, Eq. (2.24), which is obtained from evaluating the derivative of Dyson's 
equation at the mass pole. Since a two-body ladder kernel is independent of the total 
momentum P only the derivatives of the propagators contribute. Defining the quantity 

g(M^)(p) j {r(^)(g,i^)5(g+)rM(g,K)5(-g_)}^^_^^ , (A.42) 

q 



hi 
3 

T 



62 
' 3 

-63 



64 

■y2 



h = 

"2^2 



-63=3(1 

2 



(A.41) 



1 
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where the notation Q^^^'^'> includes both cases J = and J = 1, and M is the pseu- 
doscalar or vector-meson mass, the normalization condition reads: 



d 



1 . 



(A.43) 



Qt = Tp^Q^'^/S is the transverse component of (A. 42) for J = 1. Tr/j denotes a Dirac 
trace; the color-flavor trace is 1 since we use normalized color and flavor matrices. 
An analogous normalization applies to the di quark amplitudes, except for replacing 
S{—q-) — > S^{q-) in (A. 42) and attaching a symmetrization factor 1/2 in front of 
the integral [177]. In combination with the normalization (A.43), the Bethe-Salpeter 
equation completely determines the meson and diquark amplitudes on the mass shell. 

Electroweak decay constants and pion charge radius. The pseudoscalar and 
vector-meson leptonic decay constants are defined by [79] : 



ifnm^ = y/NcZ2TrD J 7' ^ S{q+)Tiq, P) S{-q.) 

1 



p2=-ml 



A 



7^5((?+)r-(g,P) S{-q^) 



P^=—nn?p 



(A.44) 



where the prefactor \/Nc emerges from the color normalization in (A. 32). 

As a consequence of electromagnetic gauge invariance, the pion's electromagnetic 
current operator which is consistent with the rainbow-ladder kernel is the impulse- 
approximation current (cf. Section 6.1): 



J^(Q^) =Tr^ / r{qf,Pf) S{qX)T''^^^{q+,Q) S{q^)T{q„Pi) S{-q.) , (A.45) 



where r^*^^ is the quark-photon vertex of App. A. 2, and the involved momenta depend 
on the photon momentum Q, average total momentum P = {Pi + Pf)/'^ and loop 
relative momentum q via 



Pij = PTQ/2, q± = ±q + P/2, qf = q+±Q/2. 



The incoming and outgoing pion momenta are onshell: P^ = Pj = —rn^ and hence 
P • Q = 0. The most general Poincare-covariant expression of the current is given by 
(cf. Eq. (A.112)) 

J^(Q2) = 2P'^F^(Q2) , (A.46) 

where Ft^{Q'^) denotes the charged pion's form factor. Charge conservation -F(O) = 1 is 
ensured by the Ward identity (A. 22). The quark-photon vertex (A. 23), together with 
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the phenomenological transverse part (A. 28) which incorporates a vector-meson pole, 
yields for the pion form factor: 

F^{Q^) = F^,Bc(Q') - i^p..(Q') e-3(^) , (A.47) 

QpX + l 

i.e., the sum of a Ball-Chiu piece and the pTrvr triangle diagram whose experimental 
value is Fp-j^^^i^Q'^ = —rr?^ = gfyj^^^ = 6.14. The latter component vanishes for zero 
photon momentum transfer owing to the transversality of the respective vertex. For 
(5^ > it reduces the Ball-Chiu contribution to the form factor and typically supplies 
~ 50% to the squared pion charge radius [60,83], extracted via 

rl = -6F;(0) = r2 Bc + F,..(0) e'^W. (A.48) 



A. 4 OfF-shell structure of the T-matrix 

Diquarks correlations in the nucleon are off-shell: it is the very requirement that the 
integration domain of scalar and axial- vector diquark propagators in the quark-diquark 
BSE (5.12) cannot exceed their pole locations without a systematic inclusion of residue 
contributions. The tool to gather information on the off-shell structure of the quark- 
quark scattering matrix T^^^ is the Dyson equation (2.18). The diquark ansatz (5.4) 
was introduced as a workaround to avoid its explicit calculation; a subsequent solution 
of the diquark BSE however only determines T^^^ on the diquarks' mass shells. As will 
be detailed below, a naive implementation of the on-shell T-matrix for off-shell (and 
in general complex) total diquark momenta ^ — M^, i.e., by using on-shell diquark 
amplitudes and free spin-0 and spin-1 diquark propagators 

^ I -""sc ' av av 

poses several conceptual problems which require closer attention. 

A manifest strategy to bridge the gap between the separable ansatz for T^^) and its 
solution from Dyson's equation is to consult Eq. (2.18) once again, namely to determine 
the diquark propagators D{P'^) upon constructing a sensible analytic continuation of 
the on-shell amplitudes T(q,P). Thereby the Dirac, color and flavor structure of the 
separable ansatz is maintained at off-shell momenta while certain features of the full 
self-consistent solution of the T-matrix are implemented as well. 

Asymptotic behavior I. The asymptotic limit of the kernel K^"^^ is a gluon ladder 
exchange. This can be inferred from the kernel's skeleton expansion: the ladder diagram 
is independent of P^ while all higher-order contributions vanish for P^ — > oo due to the 
appearance of additional quark propagators which behave as 5 — > 1 j^fP^. The Dyson 
sum (2.17) implies that the gluon ladder kernel (3.11) also dominates the asymptotic 
behavior of the T-matrix: 

T.,MV. Q, P) ~ ^ ( y)^^ (y i^^n^s . (A.50) 
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This is apparently not reproducible through the ansatz (5.4). The implementation of 
Dyson's equation in the previously described manner at best guarantees the correct 
power-law behavior T DT ^ const, of the T-matrix in the variable P^: if the diquark 
amplitudes are chosen to become constant for large (e.g., via their dominant Dirac 
structures 7^C, 7^C), the resulting diquark propagators will become constant as well. 
Phrased differently, they pick up finite parts which would appear in a T-matrix beyond 
the diquark ansatz. These contributions are suppressed on the mass shells but dominate 
the ultraviolet region. 

OfFshell ansatz for diquark amplitudes. We start from Eq. (A. 34) and for conve- 
nience discuss the offshell dependence of the diquark amplitudes in terms of the basis 
elements Tk{q, P) and t^(q', P) alone while the dressing functions are left unchanged at 
their mass-shell values: 

4 

rsc(g,P) = ^/r(g2,z = 0,-M2)i7V,(g,P)C, 

(A.51) 

r^..iq,P) = ^/r('?',^ = o,-Mfji<(g,p)c. 

k=l 

A solution of the quark-diquark BSE in the nucleon's rest frame requires boosted 
diquark amplitudes. A boost shifts the angular variable z in the diquark amplitudes' 
dressing functions into the complex plane and even outside the convergence radius \z\ < 
1 of the Chebyshev expansion which therefore no longer poses a sensible procedure. 
Setting z = (only) in the dressing functions is a reasonable approximation on the 
diquark's mass shell and similar to keeping only the zeroth Chebyshev moments. . 

The transversality condition for the on-shell axial-vector amplitude need not be 
included explicitly since it is already ensured by the transverse pole in the axial-vector 
diquark propagator: each diquark amplitude in the subsequent calculations appears in 
conjunction with the respective propagator. Likewise, the purely longitudinal off-shell 
components related to Tg i2iQjP) generated by the longitudinal projection of the 
diquark propagator which is suppressed by a factor of P^ + M^^ on the mass shell. 

In the following we will make use of the dimensionless diquark momentum variables 

Xsc := P^Ml, xav := P^M^, (A.52) 

and abbreviate both by the symbol x, where the context selects either of the two 
possibilities. Specifically, it entails x = — 1 on both scalar and axial- vector diquark 
mass shells. 

The on-shell basis of Table A.l, expressed in terms of a normalized diquark momen- 
tum P, does not provide a unique analytic continuation to off-shell momenta: 

• Using that basis, all amplitudes would be equally important at ^ oo and 
determined by their strengths on the mass shell. This is not a fundamental prob- 
lem if one keeps in mind that, in our context, an off-shell diquark amplitude is 
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no meaningful object in itself but merely an auxiliary device for constructing the 
off-shell T-matrix, especially since the requirement (A. 50) cannot be met any- 
way. Nevertheless we choose the dominant diquark amplitudes Tsc{q,P) ~ ^7^, 
Tav{q, P) ~ i'j^ to prevail in the UV such that the diquark becomes 'pointlike' 
at large diquark momenta. This feature would be explicit in an inhomogeneous 
BSE solution for a scalar and axial- vector vertex. It is accomplished by suppress- 
ing all subleading amplitudes with a factor g{x) which leaves the on-shell value 
unchanged: g{—l) = 1. We use g{x) = {x + 2)"^^/^ and, for the time being, 
keep K > as a variable (later we will set k = 2). For large n the subleading 
amplitudes are suppressed at spacelike momenta and provide support only in the 
neighborhood of the mass shell which resembles the case where only the dominant 
diquark amplitudes are taken into account. 

• The normalized basis of Table A.l includes kinematic singularities at P-^ = 
which, if not dealt with, lead to imaginary parts in the spacelike diquark prop- 
agator calculated from Eqs. (A. 59) below. Contributions ~ P and ~ z should 
behave as in the vicinity of the origin. For this reason we attach a factor 
h{x) ~ ^/x to any occurrence of P^ or z = q ■ P. The choice h(x) = —i\fx 
amounts to the replacement of the normalized by an unnormalized basis, i.e. 
P^ — > P^ /{iM). It would be sufficient to guarantee the correct on-shell behav- 
ior but renders the UV behavior completely arbitrary: some amplitudes would 
become constant, others would rise with powers of y/P^. We therefore choose 
h[x) = —i^Jxjix -|- 2) which entails 



and hence 1) = 1, /i(0) = 0, h{po) = —i. 

The final off-shell ansatz for the scalar and axial-vector diquark bases is given in Table 
A. 2. This is the construction which was used in Refs. [113,202]; it slightly differs from 
that employed in Ref. [129]. In combination with the diquark propagator of Eq. (A. 59) 
it provides a prescription for the quark-quark T-matrix which is unique on the mass 
poles and in the ultraviolet whereas its intermediate momentum behavior depends on 
the parameter k inherent in the definition of ^(x). We note that a further P^-dependent 
function attached to the full scalar or axial-vector diquark amplitude does not change 
the product Y DY since it would also appear in the diquark propagator and leave the 
T-matrix itself (and as a consequence, baryonic observables) invariant. 

Diquark propagators. With an off-shell ansatz for the diquark amplitudes at hand, 
the Dyson equation (2.18) can be exploited to obtain a consistent expression for the 
diquark propagator. Insertion of the diquark pole ansatz for the T-matrix, 




(A.53) 



T = V DV ^V^ D^'^'V'', 



(A.54) 



98 



A collection of propagators, vertices and amplitudes 



J = 


J = l 


ri = 1 
T2 = ghf 
T3 = ghz 
U = ghi[^,f] 


rf = tI^ = ghz It - q^} 

T^ = ghjt^r = ghi[q^^r-^^J'^[lr]} 

T.^=giq^ =g{q^iT-lqT'y^} 

= gh^ z = gh{\eTi^r-\ m n } 



Table A. 2: Scalar and axial- vector diquark off-shell bases, to be used in conjunction 
with Eq. (A. 51). The replacements |t — > (I — h'^z f) and — > (1 — h^z^) are implicit. 



(A.56) 



into Dyson's equation yields: 

TDf + r'^D^'T = K + KGqT Df + KGq r^'D^'^T" . (A.55) 

Successive application of quark propagator pairs Gq and scalar or axial-vector diquark 
amplitudes from the left and the right and closing the loops by integration and tracing 
yields (transitions between scalar and axial-vector amplitudes vanish because of their 
flavor traces): 

(rGor)D(rGor) = {TGoKGoT)+ 

+ {TGoKGoT)D{TGoT), 

(r°Gor^)z)^''(PGor^) = {t'^GoKGoT^)+ 

+ {T''GoKGoT^')D^'^{T''GoT^). 
We simplify the notation by introducing the shorthand notation 

n('^'^) = ^r(^)Gor('^) = -^Tr^y T^^^ SST'-''\ (A.57) 

k(•^'^) = A.fMGoKGoT^"^ = -^^^'dJ j T^^'^ SS K SST^^l (A.58) 

where we chose the letters n for normalization and k for kernel: n^'^^-' is, up to a 
minus sign, the dimensionless version of the normalization integral Q{P^) defined in 
Eq. (A. 42), but now without the dependence of the diquark amplitudes held fixed. 
Upon decomposing each axial- vector quantity D^'^, n^'^, k^^ into a transverse and 
longitudinal part (e.g. D^^'^ = DxTp'^ + Dl L'^), the above equations become 
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Figure A. 4: Defining equation (A. 59) for the diquark propagators. The first graph 
corresponds to n^/A;, the second term to n. 



These are the defining relations for the scalar and axial-vector diquark propagators in 
the quark-diquark model, n and k are scalar functions that depend on x = P"^ /M"^, 
where M = M.. or M.v. 



An ansatz for the two-loop integral. n{x) is a one-loop and k{x) a two-loop 
integral, to be evaluated for x G C since the quark-diquark BSE samples the diquark 
propagators within a parabolic domain in the complex plane. While n(x) can be 
easily evaluated, we try to circumvent the Monte-Carlo calculation of the quantity •n? /k 
for general momenta P^ and model it instead by an ansatz, a procedure which will also 
prove useful in App. A. 7 when constructing a diquark-photon vertex. What conditions 
can be imposed upon these functions? 

• By construction, the diquark propagators should exhibit timelike poles at P^ = 
— which enable the derivation of the diquark BSEs and thereby define the 
diquark masses. Transversality of the axial-vector diquark amplitude on the mass 
shell is ensured by a transverse pole in the axial-vector propagator. The pole 



conditions = -M^J 

amount to the relations 



P2 + M2^ and D-^P^ 



-Mi 



n 



{T)(-l) = fc(T)(-l) 



n 



(T) 



(-1) = 1 + ^(T)(-1 



P' + Mi 



(A.60) 



To avoid kinematic singularities, the transverse and longitudinal parts of the 
axial- vector propagator should be equal at P^ = 0: -Dr(O) = Dl{0), and for each 
single term: nr(0) = nL{0), /ct(0) = A;l(0). 



Defining A 

,2 



ni 



-1) and P := 1 + n'( — 1), one can write 



n 
~k 



\sc + PscF{x), 



n; 



Tl 

7^ = Aav + /3avFT(x), 

kr ki 



Aav + /3avFL(x), (A.61) 



where the 2-loop content has been shifted into the unknown functions F{x). The 
above conditions require 1) = 0, F'^rp^{—1) - 

loop evaluation of the integral k^^^'^^ at specific P^ 
together with Ft(oo) = Fl{oo), are indeed satisfied. The values of A and /3 depend 
on the number of amplitudes that are taken into account. For instance, if only the 



: 1, and Fr(0) = Fl(0). A two- 
E K, verifies that these relations. 
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dominant amplitude is retained: /3 = 0, since by virtue of our off-shell ansatz the 
diquark amplitudes are then independent of and n^^^) becomes the normalization 
integral, with n'^j,^ = — 1. 

An ansatz for F{x) can be constructed in the following way. Formally, if we had 
started from the inverse Dyson relation = — Go, we would have arrived at 
the expression /k ~ T T. Since the rainbow-ladder kernel is independent of the 
total momentum P and the dominant diquark amplitudes are as well, the only P^ 
dependence of that term is introduced by the subleading diquark amplitudes through 
the functions g{x) and h{x). On the mass shell, T K~^T ~ n^/k = A is real, and so 
is each single contribution ~ fiK~^Tj. Hence F(x) can only depend on even powers 
of h{x), otherwise it would be complex on the positive real axis. Consider the scalar 
off-shell basis given in Table A. 2: the general form for F{x) is 

F{x) =Fi + {g{x)h{x)f {F2 + F^hixf + F^h{xf) . (A.62) 

Exploiting the conditions F{ — \) = and F'{—\) = 1, using g'{ — \) = —Hi/2 and 
/i'(-l) = -1, leads to 

+ ig{x) Hx)f {F,ih' - 1) + F(4)(/.* - 1)} , 

which for large x becomes 

I + 2F3+4F4 

F{x ^ 00) = . (A.64) 

The numerical 2-loop result yields F{x — > cxd) ^ 1/(k + 2). F3 and F4 only emerge 
through the transverse projection of ^ in T3, so it is conceivable that they are small. 
We therefore choose the approximation: 

Fix) « '-(^(-^'^i-^^' = —Jl + . (A.65) 



K + 2 K + 2\ {x + 2Y+^ J ' 

A similar analysis for the axial- vector case (here the Monte- Carlo calculation yields 
F(x 00) ~ 1/k) leads to 

, , l-g(x)^ \ ( \ \ , 

The longitudinal parts Fl(x) are poorly constrained. For simplicity we choose Ft{x) = 
Fl{x), i.e. n^/kT = n\/kL- Note that this does not implicate ut = ni, i.e. the axial- 
vector pole still appears only in the transverse part of the propagator. 
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We summarize: the diquark propagators we use are given by 
D-\P^) = Ml { - n{x) + A,c + /?sc Fsc(x) } 

(D-l)'^'^ (P) = Mf, { - n'^'^(x) + (Aav + /?av Fav(x)) 5^^} 

where x = P'^/M'^ (M = Msc or Mav). n{x) and are the numerically computed 

quark-loop integrals of Eq. (A. 57) with the onshell values n('p)(— 1) = A and n'^y^(— 1) = 
/? — 1. The functions F(x) are defined in Eqs. (A.65-A.66), and we use the value k = 2 
in our calculation: 

With Eqs. (A. 67), (A. 51) and Table A. 2, the off-shell behavior of the quark-quark 
scattering matrix is completely determined. 

Asymptotic behavior II. By virtue of Eq. (A. 67), the asymptotic limits of the inverse 
diquark propagators are 

D-\P^ ^ oo) = Ml {Asc + /?sc/4} , 

(D-l)^" (P2 ^ OO) = Mf, {Aav + /3av/2} <5^^ 

i.e. they become constant in the ultraviolet since the quark-loop integrals n vanish in 
the UV due to the P^ suppression of the quark propagators. This is the anticipated 
result which follows from the off-shell ansatz for the amplitudes: 

rsc(g, P ^ oo) ^ /r(g2, 0, -Ml) ij'C, 

Typical calculated values for A, /? and the amplitude normalizations are: 

A,c ~ 1, Psc ~ 1, /f (0, 0, -Ml) ~ 20, 

Aav ~ 0.5, /3av ~ -0.4, /r(0,0, -Ml) ~ 10. 

Consider the propagators constructed from (A. 67) in the limit of a single dominant 
amplitude, i.e. where /5 = 0: 

11 1 / T^'^ L^^ \ 

^^^'^ = Mf Asc-n(x)' ^'^'^^^^ = Ml VAav-nT(x) + -nL{x)) ' ^^'^^^ 

A simplification which reproduces the perturbative and on-shell behavior (the error in 
the intermediate P^ region is < 30%) reads 



which explicitly demonstrates the appearance of non-resonant contributions in these 
'diquark propagators' which serve as a surrogate for a more involved structure of the 
T-matrix. Of course one could have directly extracted a pole contribution from the 
result (A. 67) to arrive at a similar form. 
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A. 5 Nucleon amplitude 

This section collects supplements to the Faddeev equation of Chapter 4. The nucleon's 
three-quark amplitude (depicted in Fig. A. 5) including its full Dirac, flavor and color 
dependence is given by 



^{p, q, P) = {^(p, q, P)ma ® + ^(P, q, P)ms ® } ® ^ • (^-74) 

The Dirac amplitudes '^Ma^ ^Ms carry four fermion indices and depend on two rel- 
ative momenta p, q and the total onshell nucleon momentum P, with = — M^. 
Their decomposition into Dirac basis tensors with Lorentz-invariant coefficient func- 
tions is stated in Eq. (4.4), and an orthogonal 64-dimensional basis X^- is presented in 
Table 4.1. ^^'^ ^Ms mixed- antisymmetric or mixed-symmetric with respect 

to the permutation group S^. The Clebsch-Gordan construction of the corresponding 
flavor tensors reads: 

'^Ma = 4i? «cr2 1 = s° (g) (uut + ddt) , 

Ta^s = --^o-ia2 (gi cr = du"!" - -i^ (g) (uut _ ddt) - J j 5^ (g) ud^ , 

where ai are the Pauli matrices and s*^, S"^'^''^ denote the (anti-)symmetric isospin-0 
and isopin-1 quark-quark representations defined in Eq. (A. 35). A projection onto the 
proton's or neutron's flavor state involves a contraction of the rearmost flavor index 
with either of the two isospin basis states u or d. 

Basis transformations. In the Dirac basis of Table 4.1 we exchanged the basis 
elements S2j, P2j (j = 1 • • - 4) with the vector- vector and axialvector-axialvector com- 
ponents Vij, Aij. The remaining basis elements of Eqs (4.9) and (4.10) are linear 
combinations of the former: 

521 = ±(P34 - P43 - A12) V31 = Si3 ± P3I T Al3 

522 = =F(P33 + P44 — All) V32 = S23 =F P41 + Si4 

523 = ±(P32 + P4I - A14) V33 = Sii T P33 ± All 

524 = ±(P42 - P3I + A13) V34 = S12 T P34 ± A12 



(A.76) 



V21 = S43 — S34 =F P12 V41 = Si4 lb P41 =F Ai4 

V22 = S33 -I- S44 lb Pii V42 = S24 ± P3I — Si3 

V23 = S32 + S41 =F Pi4 V43 = S21 =F P43 ± A21 

V24 = S42 — S31 lb Pi3 V44 = Sii =F P44 ± All 
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Figure A. 5: Momentum routing in the nucleon's three-body amplitude. 



The positive/negative signs refer to the superscripts it which we did not state explicitly. 
The corresponding relations for P, A are obtained by interchanging S ^ P, V <-> A, and 
similar dependencies hold for the Tjj, e.g.: T^^ = — 2 A^^. Eqs. (A. 76) can be verified 
by expressing internal products such as through the respective commutator and 

using 

4^A± = Tie'^'^°^P"7^75A±, with 4^=-l[^A^,^-] , (A.77) 
together with the e-tensor identities 



pa/37 £MP- =S^p{6f,^6^^-6prS. 



■ya J 



+ Oaa [OprOyp " OppOyr) y^p _ 

t t — U Ooifj . 

+ ^ar {SppSjcT — Spcr^^p) , 

Explicit implementation. The Faddeev equation is solved via iteration using similar 
techniques as in the two-body case. A 'wave function' <I> = SS ^ is extracted from the 
equation ^ = KSS^ and evaluated outside the loop integral. The subsequent inte- 
gration ^ = is carried out by calling the wave function ^ with the loop momenta 
as its arguments: 

^^''\p, q, P) = S{pb) S{pc) ^{p, q, P) , 

^{p,q,P) = X{P^)Y, [ K^^\k)(t'^^\p^''\q^''\P). ^^"^^^ 

The index a = 1 ... 3 denotes the three permutations of the Faddeev kernel, and {a, 6, c} 
is an even permutation of {1,2,3}. These steps are repeated for fixed P^ until the 
eigenvalue A(P^) converges. The value A(— M^) = 1 which defines the nucleon mass M 
is obtained upon processing the above procedure for different P^ . 

If the Dirac basis is complete, the wave function can be projected onto the same 
basis elements as the amplitude (cf. Eq. (4.4)): 



^%[p.q:P) = Y.~fl''\v\q\{z})n[p,q,PU,s, (A.79) 

i=l 
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where the Ti are given by the orthonormal basis elements . Using the orthonormality 
relations (4.23) yields coupled equations for the amplitude and wave function dressing 
functions: 

^ °=^fc (A.80) 
where the kernel /C|^'* and quark propagator matrix are the matrix elements 



4f{p, q, P, k) = Up, q, P)p^,s, K^:lpp,,,m T, {p^^\q^^\P)^,^,^^,^ 

Gij\p, Q, P) = Ti{p, q, P)f3a,5'y Gal' Pf3'-yy' 1, P) Tj{p, Q, P)a'f}',j'5 

of the quantities 

(3) 



(A.81) 



(A.82) 



and the remaining expressions corresponding to a = 1,2 are obtained by a cyclic 
permutation of the index pairs and momentum indices. 



Numerical aspects. The Faddeev equation (A.80) is equivalent to an iterated (multi- 
dimensional) matrix-vector multiplication, where in a straightforward implementation 
the kernel /C and propagator matrix Q would be computed in advance. The kernel 
depends on 9 Lorentz-invariant momentum variables: five (p'^ , q'^ , {z}) correspond to 
the outer momenta p, q, P and four (/c^, k ■ P, k ■ p, k ■ q) to the loop momentum k. It 
furthermore involves the amplitude indices i,i = 1 ... 64 and the permutation counter 
a = 1,2,3. Choosing 20 grid points for each momentum variable leads to a memory 
requirement of 20^ (momenta) x 64^ (amplitudes) x 3 (permutations) x 16 (double- 
precision complex number) Byte ~ 90 Petabyte, which is clearly beyond the capacities 
of today's computing facilities. 

The uppermost priority in a numerical optimization thus concerns the reduction of 
memory usage of the kernel /C -"^ . This can be accomplished by splitting the kernel in 
a momentum-independent component which involves all Dirac traces (to be computed 
in advance) and a momentum-dependent part which is evaluated during the iteration 
process. With the momentum alignment (4.16), on whose account pr, qt and P become 
the three orthogonal unit vectors of Eq. (4.17), the conjugate basis elements fi in 
Eqs. (A.81) effectively do not depend on any momentum at all, and the momentum 
dependence is carried by the inner elements Tj (p(") , , P) . One may rewrite the 
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orthonormal basis X^- in the following way: 



1 01 



Ti,{p,q,P) (a±75C0A+), 



/ i = 1,5,7\ 
i = 2,6,8 
i = 3 

Vz = 4 y 



(A.83) 



Abbreviating pr ^ p and qt ^ q, the momentum-dependent parts Tjj for i = 1 ... 4 
(i.e. the s = 1/2 covariants) and z = 5 . . . 8 (s = 3/2) read: 



{i < 4) 


\/6 (i = 5,6) 


V2 T,j {i = 7, 8) 


101 

10^1^'^] 
1 0|^ 

10^ 


3 1^ — 7^ 7^ 
3(/0^-^0;^) - 7^0 7^[|i,^] 
3;^ 01 -7^07^;^ 
3^01- 7^ 7^^ 


|^0;^ + 2j^0^-7^0 7^ 
1^ ^ + ^ 0|^ 



This simplifies the extraction of the momentum dependence, e.g. in the nucleon's rest 
frame via 



1 



<-i<4,2 



and similarly for the remaining ones. Hence the kernel can be written as 
lcf;\p,q,P,k) = [/Cjf (A;)] \gUp^''\q^''\P)] , 



(A.84) 



(A.85) 



where the momentum dependence carried by the Tij has been shifted into the functions 
glj. In the same way one may isolate the kernel's dependence on the gluon momentum 
k: 

47ra(A;^) 



lcl;\p,q,P,k) 



(a) /, 



k^ 



(A.86) 



This strategy greatly reduces the memory demand to < 1 GB. The impact on the 
run time due to the above multiplication (which is now processed on the fly) is still 
slightly outweighed by the time consumed to interpolate the dressing functions inside 
the integral. For this reason we drop the dependence on the angular variable zq = pt'Qt 
which, from analogy of the quark-diquark model, is expected to be small. In addition 
we perform an expansion into Chebyshev polynomials of the first kind (see App. B.2) 
in the remaining angles z["'\ that appear in the wave function coefficients inside 
the integral. The resulting run times are accessible by a parallel cluster. 
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A. 6 Quark-diquark amplitudes for and A 

Nucleon. The matrix- valued quark-diquark amplitudes ^^^p{p,P) which were intro- 
duced in Sec. 5.2 feature the following Dirac, color and flavor decomposition: 



^AB ^ ^1,2,3 
ab 



(A.87) 



p is here the relative momentum between quark and diquark and P the total nucleon 
momentum on the mass shell: = —M^. The number of basis elements can be 
inferred from the Clifford algebra. The positive parity condition for the full Faddeev 
amplitude translates into positive parity of the quark-diquark amplitudes. This entails 
a maximum number of four and twelve Dirac basis elements for the scalar and axial- 
vector quark-diquark amplitudes, with a possible basis for both cases given by 

{i,^,r,^r}, {r,p^,p^}x{i,^,r,^r}- (a.ss) 

The constraint of positive energy for the nucleon, expressed by the positive-energy 
projector A_|_(P) = (1 -|- /')/2, halves the number of basis elements via -f* A_|_ = A_|_ to 
2 Dirac matrices {1,/^} for the scalar quark-diquark amplitude and 6 matrices for the 
axial-vector one. A partial-wave decomposition in terms of quark-diquark total spin 
and orbital angular momentum eigenstates in the nucleon's rest frame [180, 254] leads 
to the orthogonal, dimensionless Dirac basis stated in Table A. 3. Its orthogonality 
relations are 

^Tr{fkTiA+} = {-1)'^Um, ^TYKr,'^A+} = (-1)^5^,. (A.89) 

The dominant amplitudes are those corresponding to the s waves ri and r^, with a 
further addition of s (rg*), p {t2, ^2*45) and d waves (r!^)- 

The flavor wave functions of the full nucleon amplitude, constructed from the 
Clebsch-Gordan prescription, are given in Eq. (A. 75). The flavor matrices in the 
quark-diquark amplitudes are their "quark remainders" upon removing the diquark 
contributions: 

t° = uut + ddt = l, 



ti = ^|dut = -lg(ai-m2 



t2 = -^(uut-ddt) 
t' = -\/i^dt = --)g(ai+zc72) 
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^>^(p, P) 


Tl = 1 
T2 = ^ 





Table A. 3: Orthogonal basis for the nucleon quark-diquark amplitudes. We abbrevi- 
ated the normalized, transverse relative momentum hy := ipr^, where pi^ = Tp^p^ 
and Tp^ is the transverse projector with respect to the total nucleon momentum. 



A projection onto the proton and neutron isospinors u = (1,0) and d = (0,1) yields 
the following flavor wave functions: 



id,-^u,o), 



Delta. Compared to the nucleon, only an isospin-3/2 diquark can contribute to the 
spin-3/2, isospin-3/2 A amplitude which excludes the involvement of a scalar diquark. 
Denoting the total A momentum by P with P^ = — M^, the on-shell quark-diquark 
amplitude of Eq. (5.10) is decomposed into 8 covariant structures [180]: 




where the basis elements include the Rarita-Schwinger projector onto positive-energy 
and spin-3/2 spinors: 

pp'^(p) = A+(P) T^" tI,"" (^5"^ - ^ 7^7'^^ . (A.93) 

A general Green function with two fermion legs, two vector legs and two independent 
momenta p and P allows for 40 possible Dirac covariants t^^^ constructed from 

S^^^il, i>, ti^f], {P^P^in X {p',P',Y] X {1, i), r,i>r]- (A.94) 

The elements PP, 7'', If and become redundant upon contraction with the Rarita- 
Schwinger projector: pP^P^ = ^ppP'^ = Q, f^A-i- = A-|-. This leaves 8 covariants 
for which a convenient orthogonal set is displayed in Table A. 4. The corresponding 
orthogonality relation reads 

ilV{afO = (-1)^+^4;, (A-95) 



108 



A collection of propagators, vertices and amplitudes 



<^>'"'{p, P) 



Table A. 4: Orthogonal basis for the A quark-diquark amplitudes, with g'^ := ipr^- 



= l/2,P^'q'' 



where a^'^{p,P) = tI^^{p,P)V^^{P) and the conjugated amplitude is defined as (the 
superscript denotes the Dirac transpose): 

= C^rf (A.96) 

The basis of Table A. 4 again corresponds to a partial-wave decomposition in terms of 
quark-diquark total spin and orbital angular momentum in the A rest frame [180,255]. 
Since there is only one spherically symmetric s-wave component {t^^)^ the A's deviation 
from sphericity amounts to an admixture of p (r2'4 g), d (rg'g y) and / {t^'^) waves which 
contribute a significant amount of orbital angular momentum to its amplitude. 

In the same way as for the nucleon, the isospin matrices of the A quark-diquark 
amplitude are constructed via removal of the diquark contributions from the full A 
flavor-amplitude. The equivalent of Eq. (A. 91) reads: 



A++ = (u,0,0), A+=f;)5d, :^u,0 



O'TT^d, T^^u , A- = (0,0,d), 



(A.97) 



where the three entries correspond to the axial-vector diquark's three symmetric isospin- 
1 states. Contraction with the diquark flavor matrices s"^'^''^ of Eq. (A. 35) yields the 
final flavor tensors corresponding to the four A states. 



Normalization. The quark-diquark amplitudes for nucleon and A are normalized 
via the canonical normalization integral of Eq. (2.24). In this context, ^ is the BSE 
solution of the quark-diquark amplitude. Go is the product of a dressed quark and 
diquark propagator, and K is the quark-diquark kernel (5.13). In contrast to the 
meson and diquark case, the kernel depends on the total nucleon momentum P such 
that the derivative dK~^ / dP^ cannot be omitted. The normalization condition for the 
nucleon is equivalent to the normalization of the proton's electric charge: G^{0) = 1. 
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A. 7 Diquark-photon vertex 

General properties. With the notation k = + k-)/2 and Q = — k-, where 
k^ and /c_ are outgoing and incoming diquark momenta, the general form of the scalar 
and axial-vector diquark-photon vertices is (cf. Fig. A. 6) 

2 

rfdq)(^' Q) = Y1 ft'^''\k\ k ■ Q, Q2) rtik, Q) , (A.98) 

i=l 

Q) = fl ft"'''\k\ k . Q, Q2) rr^^k, Q) , (A.99) 

where the set of basis elements is given by 

rf(A;,g) G {A;^Q'^}, 

rf "^(fc, Q) G {P, Q^} X {5"^, k'^k^, Q'^Q^, k^Q^, Q^k^} , (A.lOO) 

S^""" X {A;^,Q^}, (5^^ X {A;",Q"}. 

Both vertices satisfy vector Ward-Takahashi identities which reflect electromagnetic 
current conservation. They constrain the longitudinal contributions ~ by relating 
them to the corresponding diquark propagators: 

Q^rfdq)(^' Q) = D-\k+) - D-\k^) = {2FA(^-i)} , 

Q'^rf';;f(A:,Q) = D~^^{k+) - D-l{k.) = (A.lOl) 

The average S and difference quotient A were defined in Eq. (A.21). D-^{k^) is the 
inverse scalar diquark propagator; D^^{k'^) and cr{k'^) := (^D^^{k'^) — Dj^^{k^)) //c^ are 
the dressing functions of the inverse axial- vector diquark propagator: 

D-^{k) = D^^k") Tf + Dl\k^) Lf = D^\k^) - aik") k'^k^ . (A.102) 

The differential Ward identities read: 

^':,^ik,0) = ^^^^ r^f?ik,0) = ^-^4^. (A.103) 
(dq)^ ' ' dki" ^ ' ' dki' ^ ' 

The final expressions for the vertices upon implementation of the WTIs can be read 
off from Eq. (A.lOl): 

r(dq)(^'Q) = 2fc''A(B-i) +r^"r^(A:,Q) , (A.104) 
rf-f(fc,Q) = 2k^ (A(^_.)5"^ - A(.) {k-k^ + ^)) 

- U^'^kf^ + 6'"^k'') + T^" r^°^(A;, Q) . (A.105) 
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Figure A. 6: Momentum routing in the di quark-photon vertex. 



Similar to the quark-photon vertex of Section A. 2, the physical properties are encoded 
in the transverse projection of the vertices (A.104-A.105). The transverse parts include 
the elements of (A. 100) in which no is involved, implying one further transverse 
structure ~ in the scalar vertex and nine for the axial-vector vertex. In the axial- 
vector case the symmetry requirement 

r5^f(A:,g) = r5^°(fc,-Q) (A.ioe) 

leads to relations between the transverse dressing functions which reduces the number of 
independent components to six [129,256]. The transverse parts must vanish for — >■ 
and do not participate in electromagnetic current conservation. Nevertheless they give 
contributions to magnetic moments and electric quadrupole moments in the case of a 
spin-1 diquark. In this respect it is advantageous to add the transverse contribution 

+ A(,) jfc^'Q^Q'^ - ^ + ^^/^Q") I + (A.107) 

+ A(,) {fe^ (k^Qf^ - kl^Q"'^ + k-Q (j'^^fc^ - J'^^A;") } 
to Eq. (A. 105). The resulting vertex [256] may be more conveniently written as 

/ X ^ (A.108) 

a{kl) 5^1^ kl + aikl) J/^" kt] + T^" r^°^(A;, Q) . 



In the special case of free scalar and axial-vector propagators, i.e. 

D~\k^) = k" + Ml, D^\k^) = e + Ml, Dj}{k^) = Ml, (A.109) 
the corresponding vertices satisfy 

A(^-i) = A(^-i) = 1, a(/c2) = = 1, A(,)=0. (A.llO) 
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On-shell vertices I. To obtain the electromagnetic current matrix of a spin-0 or spin- 
1 particle with mass M, the respective photon vertex must be taken on its mass shell, 
where 

= -M^ k-Q = 0, (A.lll) 

and hence any transverse projector Tq' acting on k"^ can be replaced by 5^'^. This 
yields for a scalar diquark: 



(q2) = rf, Jk, Q) = 2k^Fi(Q^) . (A. 112) 



The electric form factor Fi{Q'^) is the sum of Aj-^-i) = 1 and the transverse component 
~ k^^ in (A. 104). Electromagnetic current conservation entails 

Qi'Ji'{Q'^) = 0, J'^(O) = 2k^', Fi(0) = 1. (A.113) 

The electromagnetic current of a spin-1 particle is found by applying transverse pro- 
jectors with respect to outgoing and incoming diquark momentum (corresponding to 
the transverse pole in the propagator) onto the on-shell vertex: 

Ignoring all terms which become redundant upon contraction with the transverse pro- 
jectors, the on-shell diquark-photon vertex reads [62,258,259]: 

rfdqf Q) = (^1 -J"^ + ^3 - Gu (5^"Q^ - <5^^Q") , (A.115) 

with the three Q^-dependent electromagnetic form factors -Fi, Gu and -F3. Again, 
Fi{Q'^) is the combination of A^-^-i^ = 1 and the transverse component ~ k^ 5°^^ in 

(A. 108); and current conservation guarantees the relations of Eq. (A.113). GMiQ"^) is 
the magnetic dipole form factor with Gm{0) = A* being the magnetic moment of the 
axial-vector diquark. -F3(0) is related to an electric quadrupole moment. Note that the 
on-shell projection of Eq. (A. 108) without its transverse part yields the vertex 2k^ 6"^ 
while that of Eq. (A. 105) already includes a constant contribution to the magnetic 
moment: Gm{Q'^) = a{-M^^) 



Eq. (A.115) is of limited use in the description of a baryon as a bound state of quark 
and diquark since the internal diquark is always off-shell. Minimal ansatze for the full 
scalar and axial- vector diquark-photon vertices in the literature involve Eqs. (A. 104) 
and (A. 108) which ensure the WTI for arbitrary diquark propagators, together with a 
transverse term of Eq. (A.115) including the constant diquark's magnetic moment and 
possibly its electric quadrupole moment [62,63, 184]. To account for the suppression of 
the generalized form factors corresponding to Gm and F3 for non-zero photon momen- 
tum and off-shell momenta k, the magnetic moment must be chosen much smaller 
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Figure A. 7: The fully resolved diquark-photon vertex of Eq. (A. 119). The first row de- 
notes the impulse approximation (A. 120), the second row the seagull coupling (A. 121), 
and the third row the seagull coupling where the inverse rainbow-ladder kernel is in- 
volved, Eq. (A. 123). 



than its point-like value 2 [63]. This issue is not manifest in our calculation where the 
diquark-photon vertex is constructed from a different principle, cf. Eq. (A. 119). 



Scalar-axialvector transition. Electromagnetic transitions between scalar and axial- 
vector diquarks appear naturally in the context of nucleon form factors and mediate 
a spin flip within the diquark correlation that contributes to the Faddeev amplitude. 
These contributions are not constrained by current conservation and hence purely trans- 
verse with respect to the photon momentum: 



(dq) 



0, 



(A.116) 



In analogy to the radiative p — > ir^ decay [87], the transition vertex is described by the 
Lorentz structure 



(dq)' 



IE' 



fllBpcr 



Mn 



(A.117) 



Since in this case the incoming and outgoing diquarks have different mass shells, the 
onshell values of the involved momenta are 



Ml + A£ 



2 

av 



4 ' 



k-Q 



Mi 



Ml 



(A.118) 



Resolving the diquark structure. In the actual calculation we follow an approach 
explored in Ref. [99, 182] and express the diquark-photon vertices in terms of the 
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diquark's substructure, i.e. by the "gauged" diquark propagators, cf. Eqs. (A.59), 
(A. 67) and the discussion in Section 6.1: 

rf4(A:, Q) = T^^Uk, Q) + r^G^'Cfe, Q) + Q) , (A.119) 



where each of the three terms corresponds to a row in Fig. A. 7. This is the most 
general form which is consistent with Dyson's equation for the 2-quark T-matrix. It 
satisfies electromagnetic current conservation and provides expressions for the hitherto 
unknown transverse components of the vertices in Eqs. (A.104-A.105). The indices a, 
h collect scalar (a, 6 = 5) and axial-vector (a,b = 1 ... 4) quantities as well as possible 
scalar- axial vector transitions. 

The impulse-approximation part in (A.119) expresses the coupling of the photon to 
the quark propagators within the vertex. The two diagrams from the coupling to the 
upper and the lower quark line are equal, cf. Eq. (A. 125); their sum yields: 



rZUk,Q) = -l jTr {T''iq-Q/2,k+)S{q+)T\q,k^) x 



siq. + Q) rfq)(g_ + Q/2, Q) 5(g_)] ^1 



(A.120) 



The second row of Fig. A. 7 denotes the part of the vertex which owes to the gauged 
diquark amplitudes, i.e. the seag ulls M'''^ M^'*^ which are discussed in App. A. 8: 

r^G"'(fc,Q)=-J fTv{T-{qM)S{c^+)M^^\q,k^,Q)S^{<^_) + 

271^ (A.121) 

' + M^''^(g, k+, Q) S{q+) T\q, k^) 5^(g„)} . 

The internal quark momenta in both integrals are q± = zizq + k-/2 and q'^ = ibg-|-A;_|_/2. 
In order to satisfy the Ward-Takahashi identities (A. 101) for the full diquark propaga- 
tors (A.59), the diagrams involving the gauged inverse ladder kernel need to be taken 
into account as well. The photon cannot couple to the rainbow-ladder kernel itself 
but interacts with the diquark amplitudes, hence its generic form is T'^R-^M''^^ + 
M^^'"' K~^T^. Implementation of the ansatz (A. 67) for the parts of the diquark propa- 
gators which owe to the kernel, 

D],HP') = Ml {Asc + Psc Fsc(x)}, 
(D^l)'^'' (P) = Mf, {Aav + /3av Fa.(x)}5'^^ 

yields with Eqs. (A.104-A.105) and axik^) = 0: 

Tf'{k,Q) = 2fc^M2 /Jsc A^,^, 
rr^(A;,g) = 2A;^M2,/3avAp,,, 

where A again denotes the difference quotient defined in (A. 21). Here we neglected 
possible further transverse terms and contributions from a scalar-axialvector transition. 
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On-shell vertices II. Only the impulse approximation Tj'j^p contributes on the di- 



quarlts' mass shells and hence to the diquarl^s' form factors Fi and (in the axial- vector 
case) F3, Gm- Schematically, the full vertex (A. 119) is the gauged diquark propagator: 



^ (dq) 



{T{Go-K-')ry = 

f r + f (Go - K'^) M'^ + (Go - R-^) T . 



(A.124) 



The seagull contributions vanish on the mass shell because of the diquark bound- 
state equation: KGqT = T. They are nonetheless necessary to guarantee the Ward- 
Takahashi identities for general off-shell momenta, and they ensure transversality for 
the scalar-axialvector transition vertices which can be shown using the relations (A. 39) 
for the Dirac parts of the diquark amplitudes that follow from the general antisymmetry 
relation (A. 36). 

At the physical u/d mass, the vertex (A. 119) together with the input (CI) of Sec- 
tion 3.2 yields the on-shell values /Uav = 2.7 and Ksa = 2.3. These are in the ballpark of 
values used to calculate nucleon electromagnetic form factors [62, 184] in combination 
with the ansiitze (A.104), (A.108), (A.115) and (A.117). 

Color and flavor. The color and flavor factors in Eq. (A. 119) have already been 
worked out. The color trace for all contributions is 1. The flavor-charge traces yield 
e_ = C-i- = edg/2 = 1/2 in the scalar and axial-vector case and e_ = — e+ = 1/2, 
edq = for the transition vertex (e± and edq appear in the definition of the seagulls, cf. 
App A.8). For instance, using the diquark flavor matrices Sj (i = . . . 3) of Eq. (A. 35) 
entails that the impulse-approximation contribution (A. 120) is given by 



Tt {sf Sj q} A;'''^^ + Tr {sf Q sj} Aj^'"^ = 2 TY {sf sj q} A^''^^ 



(A.125) 



where Q = diag(g„, qd) is the quark electric charge matrix, and A-^j denotes the photon's 
coupling to the upper and lower quark line, respectively. In App. A. 9 the diquark charge 
factors 



fQu + 



2IV{sfsjQ} 



Qd 



qd - 



\ 



qd-qu \ 



2qu 
qu + qd 
2qdJ 



(A.126) 



are explicitly attached to the current matrix diagrams at each occurrence of the diquark- 
photon vertex. 
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A.8 Seagulls 

Seagull contributions represent the photon's coupling to the scalar and axial- vector di- 
quark amplitudes. They are fermion-scalar- vector and fermion-axialvector- vector four- 
point functions and, together with the diquark-photon vertex, reflect the diquark's 
internal substructure. 

Ward-Takahashi identities. The seagull WTIs were derived in [183, 260] and read 

Q^M^'{q,P,Q) = +e_{r(g+,P)-r(g,P)} 
+ e+{r(g_,P)-r(g,P)} 

-ed,{T{q,P+)-T{q,P)] , 

_ _ _ (A.127) 

Q''M^'{q,P,Q) = -e+{T{q+,P)-T{q,P)} 

-e_{r(g„,P)-r(g,P)} 

+ ed,{T{q,P-)-nq,P)} , 

where q is the relative momentum, P is the diquark's total momentum, Q is the photon 
momentum, q± = q±Q/2, and P± = P±Q. The charges e+, e_ and Cdq correspond to 
quark and diquark legs (cf. Fig. A.8). For Q ^ the WTIs reduce to the differential 
Ward identities: 



M''(,.P.0)=(5^^-e„J,)r(,.P). 

Fn,.F,0) = (S^A_,„_^)r(,,p). 



(A.128) 



Dominant diquark amplitudes. In the case where only the dominant diquark 
amplitudes are retained and their dressing functions /i (either scalar or axial-vector) 
depend solely on g^, i.e. 

FseCg, P) = r{q') h'C, r^,(g, P) = r{q^) i^C, (A.129) 

the scalar WTI (and the axial-vector WTI accordingly) reduces to 

Q^M^'iq, P, Q) = e_ {T{q^,P) - T{q, P) } + e+ {r(g_ , P) - T{q, P)} = 
= Q^{e.iq + Q/4r A+ - e+ (q - Q/AY A] } . 

with the difference quotient 



(A.130) 



_ f{qi)-f{q') _ J{ql)-f{q^) 

ql-q' " (g±Q/4)-Q- ^""-'''^ 
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Figure A. 8: Conventions for the seagulls. In accordance with the Ward-Takahashi 
identity (A. 127), outgoing charges are taken to be positive, incoming charges as nega- 
tive. 



The expressions for the seagulls which are consistent with the WTI can be immediately 
read off from Eq. (A.130), cf. Refs. [62,99, 183, 184]: 





{e_ (<7 + g/4rA+ 


-e+ {q 










(A.132) 


M^'"(g,P,Q) = 


{e_ (g + Q/4)^A+ 




-Q/4)^Aj}i7°C. 





The seagulls vanish if the diquark amplitudes are taken to be pointlike, i.e., /i = const. 
and all other /i = 0. 

Full diquark substructure. The situation becomes naturally more complicated if 
the full diquark substructure, i.e. all the remaining Dirac covariants, are taken into 
account. By virtue of the WTI (A. 127), the difference quotients in {g-t, q} and {P±, -P} 
must be evaluated in all variables which appear in the diquark amplitude. The resulting 
expressions are lengthy yet necessary to guarantee the differential Ward identity. A 
violation of the latter inevitably contaminates the small-Q^ structure of the nucleon's 
electromagnetic form factors and can lead to sizeable deviations from the neutron's 
zero electric charge. 

To simplify the discussion, we transform the diquark basis of the amplitudes (A.51) 
with dressing functions fk{q'^,0, — M^), given in Table A. 2, into the general (unnormal- 
ized) form 

Tfc=i...4(g,p) = {i,f, 
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such that the diquark amplitudes read: 

4 



k=l 

r'..{q,P) = Y.fk''il'^9-P,P')ri:{q,P)iC. 



(A.133) 



k=l 



The new coefficients /l'^ and f^^ are the respective linear combinations of the those in 
Eq. (A. 51) and pick up a dependence on the variables q ■ P = \/q^V P^ z and P^. In 
the following we drop the superscripts 'sc' and 'av' for the purpose of brevity. 

In the same way as encountered in previous examples, the strategy to obtain a 
certain vertex (free of kinematic singularities) from its WTI is to write Eq. (A. 127) as 
a sum of difference quotients in one variable and extract a factor from the linear 
terms. For example, the first row in the scalar-diquark WTI yields 

{r(g+,p)-r(g,p)}(-icV) = 

k 
k 

+ J2^A('l+-P)f,{ql,q.p,p^)'^ {q^-P-q.P) n{q+,P) + 
k 

+ M'i^^ ^' (^^^(9+, P) - Tkiq, P)) 
k 

where A^'^ denotes a difference quotient with respect to the variable (•), e.g.: 

^ ^ q±-P -q-P ^ ' 

Expressing the Q-linear terms through scalar products with Q^, 

{ql-q^)=±Q-{q±Q/A), {q^. p - q. p) = ±^ , (A.135) 



allows to read off the respective vertex contributions. The differences in the basis 
elements Tk{q+,P) — Tk{q,P) only contribute for k = 3,4. Proceeding this way yields 
the following form for the scalar and axial- vector seagull vertices: 



M^'{q,P,Q) = 


(Mf (g, P, Q) + M!^{q, P, g)) iT^C , 






(A.136) 


M^^^^{q,P,Q) = 


[U^iq, P, Q) + M2^'"(g, P, Q)) iC , 
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where the M{^^"^ include the difference quotients of the diquarks' dressing functions fk- 



4(8) 



k=l 



The quantities (V±)^ and are free of kinematic singularities and given by 



{V-. 



Their limits for vanishing photon momentum are 



^ Q'^^O 1 pM dfk 



±>k 



dp^ 



d{q-P) 



The M^^"^ contain the differences of the basis elements: 



_i 4^ 

" 2 dq^' 
dfk 



dP^' 



ptj. 



(A.137) 



2 

e_ — e+ 



r{h+hr)-ed,{h+h 



2 {t^'^C/s + /4f ) + 5^"(/5 + /ef ) + {^^''i + q'^inifi + /sf )} 



+ 



e_ + e+ 



Q'^l^ifl + hf) - ed<,{7"(/2 + hi) + q'^ik + M)}l^ ■ 



Together with (A.137) this allows for a simple check of the Ward identities (A. 128). 
The final seagull vertices (A. 136) are identical to those introduced in Ref. [129]; yet 
the form (A. 133) allows for a somewhat simplified notation. The conjugate seagulls are 
obtained from 

M'^Cg, P, Q) = -dMf'fi-q, -P, Q) , 
M'^'^Cg, P, Q) = C{M'''''f{-q, -P, Q) . 



Transverse terms. In the same way as encountered for the quark-photon and diquark- 
photon vertices, further terms transverse to the photon momentum may also contribute 
to the seagull vertices. Such contributions were found to be important for the medium- 
structure of electromagnetic form factors [113]: ignoring them may cause form 
factors to rise with increasing photon momentum which is clearly unphysical. In the 
absence of a constraint on the seagulls' transverse structure, a minimal ansatz is to 
introduce a contribution which is proportional to the seagull itself: 

M^^"^ = M'^^") + m(x) T^" ikf^^") = (l - m(x)) T^" M^(") + L^'' M''^°\ (A.138) 
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where x = Q"^ /m?p (the introduction of the p-meson mass scale will be explained be- 
low). To avoid a kinematic singularity, the dimensionless function m{x) must vanish at 
zero photon momentum but is otherwise arbitrary. Via electromagnetic current con- 
servation only the transverse projection of the vertex contributes to observables, hence 
Eq. (A. 138) is physically equivalent to multiplying the full vertex Af^^") by a function 
1 — m{x) which is modeled on phenomenological assumptions. 

In the approach presented herein, the diquark-photon vertex is obtained by resolving 
the diquark's substructure which exposes its constituents, i.e. the quark propagator, 
quark-photon vertex and seagull amplitudes. Resolving the seagull structure is not 
feasible in our context, yet it is a reasonable strategy to model the transverse seagull 
part after that in the quark-photon vertex which would remain as fundamental quantity 
is such an approach. Specifically, we use the following parametrization: 

m{x) = -—^!—e-P^^^+''\ (A.139) 

which resembles (A. 28) and again includes a transverse vector-meson pole with residue 
l/(7p. The additional factor x in the numerator was implemented to leave the nucleon's 
static properties such as magnetic moments and charge radii unchanged via m'(0) = 0. 
In connection with the parameter p2 which appears in (A. 28), the values 

P2 = 0.001, Pi = 0.075 (A.140) 

optimize agreement with the polarization-transfer data for the proton's form factor 
ratio Ge{Q'^)/Gm{Q'^) and thereby enable a realistic Q^-evolution of the nucleon form 
factors. Nevertheless we emphasize that the implementation of the transverse seagull 
term (A. 138) has no noticeable consequences for the form factors' small-Q^ behavior, 
i.e. for photon momenta ^2 GeV^. 



A. 9 Nucleon electromagnetic current in the 
quark-diquark model 

In this appendix we collect the ingredients of the nucleon's electromagnetic current 
matrix in the quark-diquark model, given by the terms in Eq. (6.17) and depicted in 
Fig. 6.2. They depend on the quark-photon vertex Eq, the scalar and axial-vector 
diquark photon vertices E^J^, and the seagull vertices which have been specified in 
Apps. A.2, A. 7 and A. 8. The explicit form of the current is given by 

JapiQ') = 1 1 [^''{Pf,Pf)X'^^'^\pf,p,,Pf,P,)'^\p,,P,)]^^, (A.141) 

Pf Pi 



where Pi and Pf = Pi + Q are incoming and outgoing on-shell nucleon momenta. The 
loop momenta pi and pj are arbitrary, a, P = 1 ... 4 are quark and a,b = 1 ... 5 are 
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diquark indices. The quark-diquark amplitudes are the solutions of the quark- 
diquark Bethe-Salpeter equation (5.12). The quantity X^''^^ is given by 



,ab 



Xq^''^n2^)V(p/-pi-(l-r?)Q) + 



fi,ab 



(A.142) 



and involves the quark and diquark impulse-approximation diagrams (the left column 
in Fig. 6.2): 



5(p+)r^;(p+,p_)5(p_) D'^'iPd-) 

ti,a'b' 



■\b'b/ 



(A.143) 
(A. 144) 



and a two-loop diagram which represents the gauged quark-diquark kernel (right col- 
umn of Fig. 6.2): 

(XK)^f = D^-'{p,+) \s{p+) K^^^'%f,p„Pf,n) S{p.)] D''\p,_) . (A.145) 

The quark-photon coupling Xq""^ connects scalar (a, b = 5) or axial-vector (a, b = 
1 ... 4) amplitudes whereas X^^^ and X^"'^ additionally allow for scalar-axial-vector 
transitions. The quark and diquark momenta are: 



P- = Pi + 11 Pi, 
P+ = Pf + vPf , 



Pd- = -Pi + (1 - ^/) -Pi , 
Pd+ = -Pf + {l-v)Pf ■ 



The gauged kernel K^^'""^ which appears in (A.145) contains the exchange-quark diagram 
and the seagull contributions: 



SG 



(A.146) 



with 



,ab 



EX 



j^^^b ^ M^^'\h,pd.,Q) S''{q')T^{p2,Pd+) 

dq 



r'dqiPuPd-) S{q')T^!^{q',q)S{q) P(p2,Pd+) , 



(A.147) 



Kl^' = ^dciPuPd-) S^'iq) M^'"(A;2,Pd+, Q) , 



and momenta: 



Pd- - P+ 



Pd+ - P- 



Pi 



P+-Q 



P+-Q 



P2 

k2 



p- -q 

2 

p~ -q 
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Breit frame. For explicit calculations we work in the Breit frame where ingoing and 
outgoing nucleon have opposite 3-momenta and the photon consequentially carries zero 
energy. With r = Q"^ /{AM"^), the momenta read: 



iM 






V VI + T y 



Q 



( \ 



\Q\ 
\ ) 



p 



Pf + Pi 



iM 



Color, flavor and charge coefficients. The current matrix diagrams still have to 
be endued with color and flavor-charge coefficients. The color traces for the impulse 
approximation and exchange/seagull diagrams are given by 

^BA Sab . Sba £aed £ceb 5c d , t \ -x ao\ 

respectively. With the diquark and quark-diquark isospin matrices of Eqs. (A. 35) and 
(A. 90): Si and ti, i = . . . 3, and the quark charge matrix Q = diag(g„, g^), the flavor- 
charge matrices for the quark-photon, diquark photon and exchange diagrams read: 

5]5.,ttQtj, 5;tttj2Tr{stsjQ} , ^ttsjQsttj. (A.149) 

ij ij ij 

The traces for proton and neutron are obtained by sandwiching these expressions be- 
tween the isospinors u = (1,0) or d = (0, 1). The index range in z, j of the sums in 
(A.149) depends on the type of the quark-diquark amplitude in the initial and final 
state: for a scalar quark-diquark amplitude in the final state (i.e., on the left-hand 
side): i = 0, for an axial- vector amplitude: z = 1 . . . 3, and likewise for the index j and 
an incoming amplitude. For instance, with the four contributions: 

(Id 
2 ' 

2^/3 ' 

_2qu+_qd 
2^/3 ' 

Mu - Qd 
6 



S^S: uT (^tisoQsitoj u 

S-A: ut I^J^ttsjQs^tjj u 



A ^ S : 



A-A: ut(5;^ttsjQsttj|u 

i=i j=i 



the full exchange contribution to the proton's electromagnetic current (including the 
color factor) is given by 

J^"" = -f JM + {41 + Jl^) + Jf^ , (A.150) 
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A collection of propagators, vertices and amplitudes 



where Jgg, JsA^ "^AS "^AA denote the Dirac parts of the exchange contributions in 
Eq. (A. 141). For the seagull diagrams, the coupling of the photon to both quark lines 
and the diquark line in the seagull vertex (cf. Fig. A. 8) need to be taken into account; 
hence all occurences of e_, e+ and edq must be replaced by the combined flavor-charge 
factors 



SG: 



e_ ^ J^tfsjQsftj , J^tfQsjsft 



-J ' 




SG: 



e_ ^^tfsjQsftj , e+^^tfsjslQt 



-J ' 




and equipped with an overall color factor — 1. 



Appendix B 

Utilities 



B.l Euclidean conventions 

Throughout this thesis we work in Euchdean momentum space with the conventions 



P-q = ^Pkqk, V 

k=l 



p-p, i>=p-i, {7^7''} = 2 5'^^ i^ = {r)^. (b.i) 



A vector is spacehke if > and timehke if p^ < 0. Moreover: 



a 



flU 



7 



-7S^'7^ 



2 ' ' ' ' ^ ' ' ' ' ' ' 24 

The standard representation for the gamma matrices reads: 



7 



-ia' 
ia^ 



7 



1 

-1 



7 



1 

1 



The charge conjugation matrix is defined as 

C = 7^72, = = = -C. 
We express four-momenta through hyperspherical coordinates: 



y"^ sin( 



2/2 cos ( 



y 



\ 



\ 



I 



^ sin il) sin Q sin <f) \ 
sin sin Q cos (\. 
smi\) cos 6 

y cos ^p j 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



Should p describe the on-shell momentum of a bound state with mass M, then: 
— M^, and in the rest frame: 2; = 1. A four-momentum integration reads: 



/-/ 



d^p 
(2^ 



1 2tx 



1 1 

(2^ 2 



y dp'^p'^ j dz \/l - ^2 j dy j d<j). (B.6) 



-1 



-1 
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B.2 Angular dependence 

The Lorentz-invariant coefficients of bound-state amplitudes depend on Euclidean scalar 
products of the involved momenta. Denoting the relative momenta by pi and the total 
bound-state momentum with P, where = —M"^, the coefficients fk {pj,Zi) carry a 
dependence on radial and angular variables 



Pi, 



Zi 



cos ipi = Pi ■ P , Pi ■ P2, 



(B.7) 



In the rest frame of the respective bound state: pf E (0,oo) and Zi G (—1,1). A 
polynomial expansion in the radial arguments is not practicable as the dependence 
of the amplitudes can differ enormously for different amplitudes fk- individual zero 
crossings might appear, and the specific power laws for pf ^ 0, — > oo depend on 
the choice of Dirac covariants (e.g., normalized vs. unnormalized) and the inherent IR 
and UV properties of the physical problem. 

On the other hand, the dependence on the angular variables Zi is usually weak and 
well-suited for an expansion into orthogonal polynomials. We denote these generically 
by Yn{z), n £ INq, where the (continuous and discretized) orthogonality relations are 
given by 

1 N 

[ dz ^y{z) Y:^{z) Yn{z) = lim V fiy (Zfe) y^(zfc) y„(zfc) = 5mn (B.8) 



-1 



k=l 



and thus allow for an expansion of a function f{z) via 



N 



f{z) = Jim 5^/„K„(^ 



fn = < 



n=0 



dz^Y{z)Y:{z) f{z) 



N 



(B.9) 



lim Y.nY{zk)Y:{zk)f{zu). 



k=l 



Chebyshev expansion. The common use of Chebyshev polynomials in the bound- 
state framework is motivated by an approximate 0(4) symmetry in a certain type of 
ladder-exchange Bethe-Salpeter equations [72] . For / = 0, i.e., for s-wave ground- 
state wave functions in the context of quark models, the hyperspherical harmonics 
ynimi'fp, ^) (p) reduce to Chebyshev polynomials of the second kind in the cosine of the 
azimuthal angle ip. The Chebyshev polynomials of the first (T„) and second kind {Un) 
are given by 

Tn{z) := l 



1 + 



Un{z) :-- 



COS (n arccos z) 
sin [(n -|- 1) arccos z] 



n— 1 



Z — 
n 

n 

k=l 

2-n 



1 



cos k 



k=l 



COS 



1 

~ 2 

k IT 
n + 1 



{n > 1) 



(n > 0) 



(B.IO) 
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with To := to satisfy the orthogonality relation of Eq. (B.9) together with (B.12). 
For n > the Chebyshev polynomials of the second kind are related to the Tn{z) via 
Tn{z) = Un{z) — zUn~i{z). The first few Chebyshev polynomials (n > 0) are 



Tn{z) = {j^, z, 2z2 _ 1, _ 3z, - + 1, . . . I , 
Uniz) = {l, 2z, - 1, 8z^ - 4z, X^z^ - I2z^ + 1, . . . } . 



(B.ll) 



For the explicit expansion in Eq. (B.9) one uses Yn{z) = i"'Tn{z) or Yn{z) = i" Un{z), 
for each of which the continuous and discrete integral measures are given by 



1^*: nT{z) = l-^=L=; nT{zk) = ^, ^,=cos(fc-0^; 
2-'^: nu{z) = ^Vl^; nu{zk) = '^t~''\ z, = cos 



(B.12) 



vr ' ' ' A^ + 1 A^ + 1 



The Zk {k = 1 . . . N) which appear in the discrete versions were chosen to be the roots 
of the respective polynomials Ti^{z) or Un{z) which, for finite N > N, optimizes the 
agreement between the coefficients fn that appear in the expansion and the projection 
of (B.9) [72]. 

A Lorentz boost of the bound-state amplitude leads the domain of the angular 
arguments z = pi ■ P into the complex plane whereas the rest-frame solution of the 
bound-state equation is given only for z £ (—1,1). While a Chebyshev expansion allows 
for a straightforward analytic continuation within the convergence radius \z\ < 1, its 
convergence properties are lost for |z| > 1 with increasing A^. As a possible workaround 
one may drop the ideally weak dependence on z altogether, as it was done for the off- 
shell diquark amplitudes in Eq. (A. 51). In the context of nucleon form factors, where 
the nucleon amplitude is evaluated in the Breit frame, one may alternatively introduce 
a different expansion variable zq with \zq\ < 1 for each value of the photon momentum 
[63]. In general the problem can be overcome by solving the bound-state equation 
in each boosted Lorentz frame anew where the dependence on the complex variable z 
is substituted into two real arguments 2:1^2 £ (~1) 1) 



B.3 Limitations from the singularity structure 

For a bound-state equation which is solved in the rest frame, the radial arguments 
of the off-shell Green functions which appear in the equation's kernel become complex 
as well. In the following discussion we denote the associated momenta pi generically by 
P(x)i where X is e.g. a quark propagator, a diquark propagator or a diquark amplitude. 
For instance, the arguments of the quark propagators S{±q±) in the meson BSE (3.9) 
read 



■.= q+=q + crP, 
P(S2) ■= 9- = -9 + (1 - cr) P. 



(B.13) 
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: — kq - 


P{S2) 


:=q = 


PiD) 


■ = kd = 



While q G R^, the square of the on-shell meson momentum is complex: \/p^ = iM, 
and thus also ^'S- 

p2^^^ = |g|2 _ 0-2^2 + 2i aM\q\ q-P. 

Another example is given by the quark and diquark propagators in the quark-diquark 
BSE (5.12): 

k + rjP, 

-p-k + {l-27])P, (B.14) 
-k + {l-i])P, 

where p, k £ R^. Eqs. (B.13-B.14) can be summarized as 

Pix)=Rx + Px{v)P, (B.15) 

where Rx G R^ denotes the sum of all involved real momenta (e.g. the integration vari- 
ables) and P is the bound-state momentum which satisfies P^ = —M"^. The coefficients 
(3x £ R depend on an arbitrary momentum partitioning parameter r/ G [0, 1]. 
The complex arguments P^^x) sampled on a domain 

p^^^^ = \R\'^ - + 2if3M\R\Z , (B.16) 

where we defined Z := R - P G (—1, 1). This is the interior of a parabola (|i?| it i/3 M)"^ 
whose apex —^^M"^ depends on the bound-state mass. It is constrained by the occur- 
rence of the nearest singularity in the quantity X which defines the limiting parabola 

l'(X),max = ± i^X? , t G R+ , (B.17) 

where mx is the respective pole mass (e.g., nig of Fig. 3.6 for the complex conjugate 
poles in the quark propagator, Msc, Mav for the timelike poles in the diquark propagator 
and so forth). This leads to a restriction M < mx /\Px{v)\j ^-^d in total to an upper 
limit for the bound-state mass: 

M < /(r?) := min | , (B.18) 

Physics must be independent of rj, as explicitly demonstrated in [86, 180]. As a con- 
sequence one can fix the momentum partitioning r] to that value ryo which maximizes 
the upper boundary /(?/) for the calculable mass M. Assuming isospin symmetry with 
equal current-quark masses, Mgc < Mav, ^sc < 2mq and a singularity-free diquark 
amplitude, exemplary values are: 

• Meson/diquark BSE: 770 = 1/2, M < 2mg , 

• Quark-diquark BSE: r/o = mq/{mq + Mgc), M < niq + Msc > 

• Faddeev equation: tjq = ^/3, M < 3mq . 
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If the bound-state momentum P in Eq. (B.15) is replaced by the Breit momentum with 
= — (M^ + the above considerations can be applied to the nucleon form 

factor diagrams as well. The result 



implies an upper limit for the photon momentum if the nucleon mass M is already 
known: 

• Quark-diquark BSE: < 4 [(^^ + M^^f - M^] , 

• Faddeev equation: < 4 [{Zniqf - M^] . 

This is typically in the range of several GeV^ and explains the restrictions encountered 
in Section 6.3. Refined numerical methods, such as an inclusion of pole residues, are 
necessary to circumvent these limitations and establish a connection with the large-Q'^ 
domain in the form factors. 




(B.19) 
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